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Abstract 



I study a model for a massive one-dimensional particle in a singular periodic potential that is 
receiving kicks from a gas. The model is described by a Lindblad equation in which the Hamiltonian 
is a Schrodinger operator with a periodic 5-potential and the noise has a frictionless form arising 
in a Brownian limit. I prove that an emergent Markov process in a semi-classical limit governs the 
momentum distribution in the extended-zone scheme. The main result is a central limit theorem for 
a time integral of the momentum process, which is closely related to the particle's position. When 
normalized by frz , the integral process converges to a time-changed Brownian motion whose rate 
depends on the momentum process. The scaling t* contrasts with which would be expected 
p^j ■ for the case of a smooth periodic potential or for a comparable classical process. The difference is 

a wave effect driven by Bragg reflections that occur when the particle's momentum is kicked near 
. the half-spaced reciprocal lattice. 
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Mathematical models of a quantum particle in a periodic environment have been used to describe an 
electron in a metal and, more recently, an atom in an optical lattice. One topic of mathematical and 
physical interest is the transport behavior for the particle in the periodic environment. The periodic 
situation stands in contrast with the random or quasi-periodic situation, which may exhibit Anderson 
localization [3]. Another important topic is the study of the motion of the particle when acted upon 
by a static force. Zener predicted an electron in a metal would exhibit some periodic motion when a 
constant force was applied |41| . This behavior, called Bloch oscillations, is related to Bragg scattering 
and has been observed experimentally in conductor superlattices [18]. More recently, atoms in optical 
lattices have provided an analogous setting in which it is possible to measure Bloch oscillations with 
fewer noise effects [5] . 

The current article studies a suppressed dispersion effect that, like Bloch oscillations, is generated 
by a combination of outside forcing (in this case from a noise) and Bragg scattering in a periodic 
potential. My model concerns a one-dimensional massive quantum particle in a periodic singular 
potential that receives random momentum kicks (e.g. from a gas of light particles). The massive 
particle effectively does not "feel" the potential except for infrequent instances when its momentum 
is kicked near an element of the half-spaced reciprocal lattice of the potential. Near the lattice values, 
the particle's momentum has a chance of being reflected, and these reflections in momentum occur 
often enough to inhibit the motion of the particle. I imagine the model to describe an atom in a very 
singular one-dimensional optical potential. In the physics literature, the article [21J discusses Bragg 
reflections of atoms from a weak a one-dimensional optical potentials. The articles |35] report 
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the experimental observation of Bragg scattering in atoms with lower kinetic energy through optical 
potentials. 

My mathematical starting point for modeling the particle is a quantum Markovian dynamics 
generated by a Lindblad equation. In the following section, I introduce the dynamics, state the main 
theorems, discuss some background for the model, and make conjectures for a similar model. Section[3] 
contains an outline for the proof of the central limit theorem that is the main mathematical result 
of this article. Section [5] contains a proof that the probability density of the extended-zone scheme 
momentum behaves approximately as an autonomous Markov process when the Hamiltonian dynamics 
operates on a faster scale than the noise. Section [5] connects basic facts from the original quantum 
model to the limiting Markovian dynamics for the momentum process. Section [6] contains the details 
for the proof sketched in Sect. I show that a time integral of the momentum process, when properly 
rescaled, converges in distribution to a variable diffusion process whose diffusion rate depends on the 
absolute value of the momentum. 

2 Results and discussion 

2.1 The model and statement of the main results 

Let ^i(L 2 (IR)) be the space of trace class operators over the Hilbert space L 2 (W). I begin with a 
quantum Markovian dynamics in which the state of the particle, as expressed by a density matrix 
Px,t G Bi(L 2 (^)) , evolves according to a Lindblad equation 

j t Px,t = -{[P 2 + V, Px , t ] + *( PA)i ) - ^*(I),p X) t} (2-1) 

from the initial state pa,o = P- in this equation, P = — i^ is the momentum operator, V is a periodic 
5-potential (i.e. Dirac comb potential) with strength a > and period 2tt, and \£ : /3i(L 2 (K)) is a 
completely positive map describing the noise acting on the system and having the form 

*(/>) = / dvj(v)e ivX pe~ ivX , (2.2) 

where p £ ,6i(L 2 (R)), X is the position operator, and j(v) £ L 1 (M) is the rate-density for momentum 
kicks of size v. I will assume the rates satisfy j(v) = j(—v) and a = J M dv j(v) v 2 < oo. In (12. ip . 
is the adjoint map \P* evaluated for the identity operator I on L 2 (IR), and it happens that = 1ZI 

for 1Z = f R dvj(v) in my case. 

Equation (|2.ip describes a quantum particle in dimension one evolving in a potential V and 
receiving random momentum kicks v with rate-density j(v). The noise is effectively frictionless, since 
intuitively, the rate of momentum kicks does not depend on the current momentum of the particle. 
This excludes the possibility of energy relaxation in the model, and there is a linear rate of growth for 
the mean energy of the particle: 

Ti[ PX}t {P 2 + V)} = Tr[p(P 2 + V)} + a t. (2.3) 

By Bloch theory, the Hamiltonian H = P 2 + V has continuous spectrum and decomposes through 
a fiber decomposition of the Hilbert space over the Brillouin zone <j) E [— |, ^) as 

L 2 (R)= f & dW<f>, H^ = H\ H ^, 

J i-l\) 

where the Hilbert spaces are canonically identified with L 2 ([— 7r,7r)), and the restriction of H to 
the <^-fiber is a self-adjoint operator H^. The operators have a complete set of eigenvectors Vv^' 
n G N with eigenvalues E n ^ satisfying 

< E n ^ < E n+1 ^, E n ^ — > oo as n — > oo. 
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Through the extended-zone scheme, the eigenvectors ip n ,if> can be associated with a collection of eigen- 
kets \k) Q parameterized by k G R such that 

H= I dkE(k)\k) QQ (k\, 

JR 

where the dispersion relation has the form E(k) = q 2 (fc) for the anti-symmetric, increasing function 
q : R — > R satisfying the Kronig- Penney relation 

cos(2vr/c) = cos (2vrq(A;)) + — ^— sin (2irq(k)) (2.4) 

for k G R — ^Z and q(§) = § for n G Z (see (I4.ip for the corresponding Bloch functions in the 
position representation). The Bloch structure and my conventions are discussed in Appendix [Bl The 
dispersion relation essentially has the form E(k) ~ \k\ 2 + for \k\ 3> 1 except for values of k in small 
neighborhoods around the lattice ^Z, where E{k) makes jumps E(^+) — E(^—) ^ 0. The kets \k) Q 
also have discontinuities at values k G ^Z: lim e _>o \ k — e) Q / |& + e) Q . 




-3/2 -1 -1/2 1/2 1 3/2 



Figure 1: Qualitative plot of the dispersion relation in the extended-zone scheme. 

The first result is concerned with the limiting behavior as A — > of the diagonal distributions in 
the extended-zone scheme representation: Q(k\p\ tt \k) Q = D\ jt (k). I show, in a sense defined below, 
that D\ t converges for small A to the solution V t of a classical Markov process 

^Vt(k) = L(V t )(k) = J dk'^J(k,k')V t (k')-J(k',k)V t (k)), V (k) = p(k,k), (2.5) 

where the rates J(k,k') are determined by the rates j(v) and values K v (k,n) (defined below) through 
the formula 

J{k, k') :=J2j(k-k' -n)\K k „ k ,„ n (k',n)\ 2 . (2.6) 

The values n v (k,n) G C arise as coefficients in the formula 

^2 Kv (k,n)\k + v + n) Q = e ivX \k) Q , k,v£R. (2.7) 

The fact that e wX \k) Q is a combination of the \k + v + n) Q , n G Z is a consequence of the fiber 
decomposition. By the unitarity of e lvX , the coefficients satisfy Ylnez \ K v(k,n)\ 2 = 1, and the process 
has a constant escape rate: 1Z = J R dk' J(k', k). 
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Theorem 2.1 (Freidlin-Wentzell/semi-classical limit). Let p\j satisfy \2. p\ t t(k,k) = D\j(k), and 
T>t be the solution of \2.5}) . There exists a C > such that for all t > 

\\Dxj-VtW! < CX(l + t). 

Also, for t < X, the norm of the difference is bounded by a multiple of t. 

Let K r be the Markov process satisfying the master equation (j2.5j) and define the integral func- 
tional Y t = £drK r . My main result concerns the limiting distributional behavior for the processes 

1 5 

(at)~2 \K s t\ and t~iY s t, s G [0, 1] for t > 1. I will make the following technical assumptions on j(v): 
List of rate assumptions 2.2. There is a p > such that 

1. f R dv j(v) e a \ v \ < p for some a > 0, 

2. sup_i<0<! £„, eZ j(# + I) < M» 

3. mf ve[ _ 1A] j(v) > IT 1 . 

The theorem below states that the processes \t~^K st \ converge in law as t — > oo to the absolute 
value of a Brownian motion and t~iY st converges to a time-changed Brownian motion whose rate of 
diffusion emerges as the limit law of ^ 1 1 i a K s t\ 3 for v := TZa. It is clear from the above statement 
that the process t~^K s t itself does not behave as a Brownian motion; otherwise the appropriate scaling 
for Y s t would be £5, and the limiting process would be differentiable rather than diffusive. 

Theorem 2.3 (Main result). Assume j{v) satisfies List \2.2\ Let Kt be the Markov process whose 
probability densities qt satisfy \2. 5|) . qo have finite second moment, and Yt be the time integral of K r 
up to time t. As t — > oo, there is convergence in law with respect to the Skorokhod metric 

((o-t)-3\K st \, o--ivh-iY st ^ =J> (|B S |, J dB' r |B r |§), s G [0, 1], 

where B s and B' s are independent standard Brownian motions. 



2.2 Further discussion and background 

2.2.1 The Lindblad dynamics, the noise, and the Hamiltonian. 

Introductory material on Lindblad equations can be found in [3]. Some basic mathematical questions 
regarding existence and uniqueness of solutions to Lindblad equations with unbounded generators 
are not completely understood except for specific classes such as those for which the generator is 
translation covariant [27]. Because the Hamiltonian part of the Lindblad equation (|2.ip is unbounded, 
the mathematical definition of a solution to the Lindblad equation is less direct than the bounded case. 
I discuss the rigorous definition for the dynamics and related technical issues in Appendix [XJ In the 
case discussed here, these issues are not interesting or challenging, since the noise term ^ : ^i(L 2 (M)) 
is bounded and iff* (I) = 1ZL commutes with the Hamiltonian. 

Consider a Levy process with density j(v), and let ti,...,t n and v\,...,v n be the Poisson times 
and increments of the Levy process up to time t. The state for the particle p\ :t = <l?A,t(p) a t time t is 
equal to 

$ x +(p)=E\u x>t (0 P Ui t (0*1, (2.8) 
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where the expectation is with respect to the law of the Levy process, £ = (vx, ti; Vi-,ti\ ■ ■ ■ ) G (RxR + )°° 
is the full sequence of random events, and the unitary operator U\ : t(£) '■ L 2 (R) is defined by the product 

Ux>t (0 := e-^ ff e w • • • e -^^^ e 4^. (2.9) 

The construction of the maps <I>A,t thus only depends on the existence of the unitary groups e~ irH , 
r£l The equation (|2.8|) implies the trace for the state is preserved such that Tr[px,t] = Tr[p] = 1, 
since the expression on the right is a convex combination of unitary conjugations of p. 

A noise of the type appearing in (|2.ip was originally introduced as a phenomenological model for 
the study of wave collapse in quantum mechanics [23]. It was later derived in [22] starting with a 
heuristic scattering analysis that was meant to model an interaction of a test particle with a gas in 
the limit that the test particle has much greater mass than the gas particles. This scattering analysis 
was clarified in [29], which yielded a minor correction by a unitless multiplicative factor in the final 
expression for the jump rates j(v). Also, the article [26] contains a mathematical derivation for a 
noise of the form (|2,ip through a singular coupling limit of a simple system-reservoir Hamiltonian 
dynamics. The noise model has been discussed in relation to experimental frameworks in matter- 
wave optics [21 [30] and appears in other discussions of decoherence [30]. See [39l Sec. 7.1] for the 
connection of the noise with a quantum linear Boltzmann equation in a large mass limit. A similar 
frictionless noise with some spatial dependence recently appeared in [33] to model the dampening of 
Bloch oscillations for an atom in an optical lattice. 

It is clear from the mean energy growth f|2.3|) that the model for the noise is transient in nature. 
The classical analog of the noise map 

p _> v|/( p ) _ 2 ~ 1 {p, **(!)}, p e B\ (L 2 (R)), 

appearing in the Lindblad equation (|2.ip is given by the map 

T(x,p)— > ( dvj(v)(T(x,p-v)-T(x,p)), (2.10) 

for joint position-momentum densities T S L 1 (M 2 ). The association of (|2.10p with the quantum noise 
can be justified by looking at the Wigner representation W p (x,p) of p. Of course, most quantum 
noises are not so readily identifiable with classical analogs. Equation (|2.10p describes the momentum 
as receiving random kicks v with rate-density j(v). Based on this evidence, the momentum undergoes 
an unbiased random walk, and the ensuing stochastic acceleration explains the mean energy growth 
in (I2T31) . 

The Hamiltonian H = P 2 + ot^2, n&Ij 5{X — 2irn), a > is defined as a particular self-adjoint 
extension of the symmetric operator — with domain consisting of all functions / E H 2,2 (R) n 
{g | g(2im) = 0, n £ Z} having two weak derivatives and taking the value on the lattice 2ttZ. The 
domain of the self-adjoint extension is the space of functions / G H 2,1 (R) n H 2 ' 2 (R — 27rZ) that have 
one weak derivative in the domain R and two weak derivatives in R — 27rZ, and that satisfy 

a/(27rn) = -^(2vm+) - -f-(2vrn-). 

ax dx 

The Dirac comb is a limiting case of the Kronig- Penney model |34j . which is a periodic Schrodinger 
equation in which the potential has the form V(x) = V X^nez Ip^n-a, 2irn+a] f° r a < 2vr. The limit 
connecting them is V — > oo with 2Va = a. The Kronig-Penney model has been used to model the 
transport of electrons through a crystal. One computational advantage of these models is that there 
are closed equations determining the spectral values and the form of the eigenkets. Both the periodic 
5-potential and the Kronig-Penney model are discussed in [1|. Some general theory regarding the 
structure of periodic Schrodinger equations can be found in [TUSSIE]- The articles [23j[9l[l5] contain 
recent results on the dispersion of wave packets evolving according a periodic Schrodinger equation 
with a smooth potential. 
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2.2.2 Bragg reflections 

For a high momentum quantum particle in a periodic potential, the dominant behavior is simply 
transmission over the potential. However, there is a lattice of wave frequencies around which the 
particle is likely to be reflected by the potential. For my purpose, the difference between a smooth 
periodic potential and the Dirac comb is the proximity a high momentum particle must have to a 
lattice momentum to experience reflections. These zones are wider for the Dirac comb, and it is 
possible for the test particle to score accidental reflections in the process of colliding with the gas. A 
more immediate indication of a contrast between the Dirac comb and a smooth period potential is in 
the size of the jumps in the dispersion relation E(k) at momenta k G ^Z (as pictured in Figure 1). 
Unlike smooth potentials, where the energy gaps vanish as k — > oo, the gaps for the periodic 5-potential 
approach the constant value ^. 

Due to the spatial translation symmetry of the potential by 2tt, the kets \k) Q can be written as 
discrete combinations of the momentum kets \k + n) for n G Z: 

\k) Q = Y,v(k,n)\k + n), (2.11) 

where ^nez l 7 ?(^) n )| 2 = 1- For \k\ 3> 1, the n = term dominates the sum except when k is close to 
a lattice point ^ G ^Z, in which case the — n term becomes nonnegligible. When k ~ ^n, the value 
k — n ~ —k is approximately a reflection of the momentum. A high-momentum plane wave tuned near 
a lattice frequency will be driven by the Hamiltonian evolution to oscillate though a cycle of quantum 
superpositions between the original wave and the reflected wave. These are called Pendelldsung 
oscillations, and the corresponding wave velocities E'(k) are small for momenta k close enough to an 
element in |Z to exhibit the Pendellosung oscillations (see also in Figure-1). A plane wave with a high 
momentum sufficiently away from any lattice value will transmit nearly freely through the potential 
with velocity E'(k) ~ 2k. I will refer to the regions near the lattice momenta in which nonnegligible 
reflection occurs as the reflection bands (see |21|). For the Dirac comb potential, the widths of the 
reflection bands scale as oc n - 1 for § G ±Z with \n\ > 1. 

2.2.3 The small A limit 

The regime of A 1 in (|2.ip should be considered as a semi-classical regime in which phase oscillations 
generated by the Hamiltonian term occur on a faster time scale than the mean time between collisions 
with the gas. As it appears in (|2.ip . the parameter A takes the place of 2Mh, where M is the mass 
of the test particle. Since the noise arises in a limit in which M is large, the value h~ l must be even 
"larger" to make A < 1. A more honest comparison of relevant scales requires that I include more 
physical parameters in the model such as the period and strength of the <5-potential, and I do this for 
a slightly richer noise model in Sect. 12.31 

Without the noise, the dynamical evolution generated by the Hamiltonian in the extended-zone 
scheme representation is 

Q (k'\e^Xp)\k) Q = e-^')+W) Q (k'\p\k) Q (2.12) 
for p G Si(L 2 (lR)). The noise map ^ operates in the momentum representation as 

(k'\^(p)\k) = [ dvj(v)(k' -v\p\k-v). (2.13) 

Since the ket \k) Q is a combination (12. lip of the kets \k + n) for n G Z, the equations (I2.12|) and (I2.13P 
imply that only the values Q{k'\p\k) Q with k! — k G Z interact dynamically with the diagonal values 
Q {k\p\k) Q . The relevant phase velocities are \~ 1 ( y E{k') — E{k)^ for k' — k G Z. Due to the non-vanishing 
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energy band gaps for the Dirac comb potential, there is a non-zero minimum A s m ; n for the phase 
speeds, where 

s min := inf \E{k') - E{k)\ > 0. (2.14) 

fc'-fcez 

These phase oscillations thus have period proportional to A, which is small compared to the mean time 
7£ _1 between momentum kicks. This provides the mechanism through which the diagonals Q (k\px,t\k) q 
tend to an autonomous evolution as A — > 0. 

The A — > limit connecting the diagonal of pxt hi the basis of kets diagonalizing P 2 + V to a 
Markovian dynamics is analogous to a Freidlin-Wentzell limit |20t Ch.8] for a classical Hamiltonian 
flow perturbed by a weak noise. In the Freidlin-Wentzell limit, a Markovian dynamics emerges on the 
energy graph in a scaling limit combining a noise of strength A -C 1 and long time intervals oc A -1 . 
The "energy graph" is the collection of connected level curves for the Hamiltonian. My dynamics fits 
the Freidlin-Wentzell limit description with a simple change of time variable s := Xt: 

—p~x,s = -i[H,px, s ] +\(y(px,s) ~ ^{**(I),px, s }) for p x ,s ■= Px,f- 

The connected level curves of the Hamiltonian and the kets \k) Q have parallel roles in imprinting 
the Hamiltonian structure on the limiting Markovian dynamics. A similar study for a time integral 
of a momentum-related quantity after a Freidlin-Wentzell limit has previously been performed for a 
classical model in [25]. This model is a Markovian dynamics for a particle in a periodic potential 
with a white noise, and the authors have shown that the spatial diffusion and Freidlin-Wentzell limits 
commute. 

Theorem 12.31 is a statement only about the classical process Kt and its time integral Yt = J*q drK r . 
This result has no direct consequence for the original model, and in particular, I have not proven 
that the position distribution dx,t(x) ■= (xlpxjlx) itself exhibits subdiffusion for A <C 1. However, I 
believe the suppressed dispersion effect holds in the original model for any hnite A. In other words, 
dxt( x ) has spread oc when the Dirac comb is present rather than the scaling oc that holds other- 
wise (see Appendix [C]). The effect depends on the occasional development of Hamiltonian-generated 
quantum superpositions between momenta in the reflection bands and their reflected counterparts. 
As the particle is kicked out of a reflection band, the quantum superposition collapses into a classical 
superposition. 

2.2.4 Energy submartingales 

I chose the Dirac comb for this study as opposed to some other singular periodic potential because there 
are convenient closed expressions for the eigenkets \ k) Q in the position and momentum representations. 
This is helpful for estimating the coefficients K v (k,n) appearing in the definition (|2.6p of the jump 
rates for the Markov process Kt- Nevertheless, it is easier to perceive certain critical features of the 
classical process Kt by taking a step back into the original quantum framework than through the closed 
formulas. For instance, E^{Kt) is a submartingale. It follows that E(Kt) is also a submartingale, 
although it can be additionally learned from the quantum setting that the increasing part of its Doob- 
Meyer decomposition is at. That in not surprising given the mean energy grown (12. 3j) . The process 
Ez(Kt) is useful for showing Thm. 12.31 since \Kt\ ~ E^{Kt) when \Kt\ ^> 1 by the approximately 
parabolic shape of the dispersion relation, and because the martingale structure allows E%(Kt) to be 
treated (see Thm. [3TT]) . An equally important application of the submartingale structure is to show 
that the process \K r \ spends most of the time interval r £ [0, t] with large values on the order for 
which the Bragg scattering is dominated by only transmitted and reflected waves. 
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The proof that (Kt) is a submartingale is based on the Heisenberg representation of the equa- 
tion (|2,8|) . The stochastic operator- valued process 

1 i 
is a positive submartingale in the sense that for all / G D(H2) C L 2 (M), the process (f\H?J£)f) is a 

positive submartingale. Section [5] contains the analysis of the various relevant operator martingales. 



2.2.5 Heuristics for the limit theorem 



The Kolmogorov equation (|2.5p determines a pseudo-Poisson process K r whose jump times are de- 
termined by a Poisson clock with rate 1Z. The jumps are a sum of two contributions: one coming 
directly from a particle collision, and another being a lattice- valued Bragg scattering. More precisely, 
a jump w from the starting momentum k G R is equal to w = v + n, where the component v G R 
has density Hr^j^v) and the component n G Z has conditional probabilities n)| 2 when given v 

and fc. I refer to these as the Levy and lattice components of the jump. When \k\ 3> 1, most of the 
lattice jump probabilities \n v (k,n)\ 2 are negligible except for n = and a few values corresponding to 
reflections in momentum. The high-energy behavior is dominant, since by earlier discussion, \K r \ will 
typically spend most of a time interval r G [0,t] with \K r \ oc 3> 1. An idealized picture emerges in 
which the process makes a Levy jump i> from k with an additional optional jump -n£Z for 

n = 2(fc + u-6>) with 6 = k + v mod -, 0G [--,-), 

2 L 4 4' 



which occurs with probability 



a 2 



R_(2- 1 6»n) := — =■ (2.15) 

1 ; Svr 2 (2-i0n) 2 + ^ 

If this extra jump occurs, the resulting value k + v — n ~ —k — v is approximately a reflection of the 
value it would have had otherwise. 

The probability (12. 15ft of a momentum reflection decays when \8\ 3> |n| _1 . Hence, when \k+v \ S> 1, 
the value k + v must land near a lattice value in order to have a good chance of reflection. The 
integral 

,i 

* d6 R_ {2~ x 9n) « (2-16) 

i 

serves as an effective reflection probability if the particle is dropped randomly in the cell centered 
around ^. The factor of 2 in (|2.16p normalizes the integration. This suggests the number of reflection 
times N r is approximately a Poisson process with a rate depending on K r as ~ lZa\K r \~ l = v\K r \~ l 
(since the Levy jumps occur with rate 1Z). In other words, the Poisson rate of reflections is inversely 
proportional to the absolute value of the momentum. Over a time interval where \K r \ cx t?, the 
average time between reflections will be oc t2. Hence, if \K r \ behaves as the absolute value of a 
random walk, the number of reflections over a time interval [0, t] will be (at least) on the order of 
tz, and t~2K s t itself can not have a limiting distribution. However, it is reasonable to expect a limit 
theorem for the process Y r , since the sign- flipping of K r is smoothed by the time integration. Roughly 
speaking, Yt can be written 



ft N * fTn + l 

y t « / dr(-l) N "|K r |«^(-ir / dr\K, 



n=0 

/"I n J2{-l) n t-2 / dr|t~2i^ r |J, (2.17) 

n=0 
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where r n are the reflection times. By the considerations above, the integrand \t 2 K r \ and normalized 

interval t~^{r n+ i — r n ) will both be 0(1). Since the number Nj of terms in the sum is on the order 

i I i 

0(t2), a scaling factor of N 4 2 oc is appropriate if I expect central limit theorem- type cancellation 

5 

among the summands in (|2.17p . Hence, it is reasonable to expect t~iY s t, s G [0, 1] to have a nontrivial 

I — - |3 I — — I 

diffusive limit. The diffusion rate will be proportional to \t 2 K s t\ with two powers of t 2 K s t\ 
coming from the integrand (I2.17h . and another factor coming from the less frequent reflections (and 
thus diminished cancellation) that occur when the momentum is large. 



2.2.6 Related limit theorems 



The article [12] is a classical analog of the current work in which the periodic potential is continuous 
and the noise is frictionless. There, the rescaled momentum process t~^K s t, s G [0,1] converges in 
distribution to a Brownian motion, and thus the position process simply converges to the integral of 
a Brownian motion. Also for a classical model, the articles [131 [TO] work to control and characterize 
a periodic potential as a perturbative contribution to a dissipative dynamics driven by a linear Boltz- 
mann equation in the limit of large mass for the test particle. The noise there is analogous to the 
quantum noise discussed in Sect. 12.31 

The study of the limit law of t~± L s drK r for large t G K + fits under the general category of 
central limit theory for integral functionals of Markov processes. It is not covered by previous results 
that I am aware of, since the process K.% is null-recurrent and systematically makes large jumps in 
the form of sign-flips from arbitrarily high values in phase space (see [28] for martingale limit theory 
relevant to a broad class of null-recurrent situations). A simplified version of the problem is given by 
the following: let the Markov process K[ make jumps at times determined by a Poisson clock with 
rate 7Z and transition density T(k' , k) from k to k' given by 



T(k',k) 



lZr x i(k! -k) || ± (k' - k)\ > t^t for all n G Z- {0}, 

(2K)- 1 J2j{± k ' Tk) |f ±{k' -k)\ < A for some n G Z-{0}. 



In words, the process first makes a jump from k to k + v with density 1Z j(v), and if \k + v — § | < A 
for some n G Z — {0}, then the sign either flips or remains the same with equal probability 1/2. 
Note that the process \K[\ has the same law as the absolute value of a Levy process with rates j(v). 
The limit statement of Thm. 12.31 holds with Kt replaced by K' t . If the reflection bands around the 
lattice values 8 G ^Z in the simplified model above are replaced by bands with diameter a|f \~ 2 ^ for 
< i? < 1, then the limit law will be 

(aty^K'J, a-s-fi/ir- 1 -!^) (|B S |, J dB' r \B r \ 1+ ^, s G [0, 1], 
where Y( = J * drK' T . 

The Kronig- Penney model lies at the boundary •d = 1. I conjecture that if the periodic ^-potential 
is replaced by a continuous potential, the limiting behavior will agree with the classical case. The 
Kronig-Penney model will have an intermediary behavior in which t~? K s t converges in law as t — > oo, 
but the limiting process is not a Brownian motion due to a number of random reflections over the 
interval s G [0, 1]. The limit of the time integral process will be differentiable rather than diffusive. 



2.2.7 Further mathematical questions 

The first question is whether the original Lindblad model (|2.ip actually exhibits suppressed spatial 
dispersion for a fixed value of A. Within the current program of passing through a Freidlin-Wentzell 
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limit, there is the mathematical challenge of beginning with a more sophisticated noise that generates 
energy relaxation for the test particle (see Sect. I2.3p . Also, it would be interesting to compare the 
Dirac comb with other singular periodic potentials and to see how the situation changes in higher 
dimensions. 



2.3 Analogous conjectures for a dissipative model 

In this section, I introduce a related quantum Markovian dynamics complex enough to include energy 
relaxation. This material is presented with the intent to broaden the reader's perspective on the 
original model. In other words, this is a continuation of the discussion in the last section and does not 
concern the mathematical results of this paper. The model is a one-dimensional version of the quantum 
linear Boltzmann dynamics discussed in the review [39] , which models a test particle interacting with 
a dilute gas of distinguishable particles. The one-dimensional case certainly can not be derived from 
first principles, since even the classical one-dimensional linear Boltzmann equation does not arise in 
a low density limit from a microscopic Hamiltonian model for a test particle interacting with a gas. 
Analogous mathematical objects in this section to those previously defined will be denoted with a 
tilde, and the meaning of symbols introduced here will reset in future sections. 

Let the state of the particle at time t G M+ be given by a density matrix p t E Si(L 2 (M)) whose 
evolution is determined by the Lindblad equation 

±p t = -~[H,Pt] + Wt) ~ Pt}, (2.18) 

where po = p, and the Hamiltonian H and the completely positive map \I/ are defined below. The 
Hamiltonian H is a Schrodinger operator with Dirac comb potential 



2M 

where M is the mass of the test particle, a is the spatial period, and a > is the strength of the Dirac 
comb. Let 77 be the spatial density of the gas, m be the mass of a single gas particle, and D\(p re i) be 
the reflection coefficient determined by the interaction potential between the test particle and a single 
gas particle. For the noise map I take 

*(P) = / dq\q\ IW&lW L q (P)pL q (P)e--4x, (2.19) 

where L q {P) is the multiplication operator in the momentum representation with function 



1 ■ 

■2 



L q (p) = U(2" 1 (l + —)q + —p) for p(p) 

V M M ' (2-Kmp- 1 ) 

The operator *$>*(I) is also a function of the momentum operator given by 

**« = ™ K>i 2 ^ + > + t,^ - £ < p >- 

The function E (p) serves as an escape rate for getting kicked out of the momentum p. 

I will assume the interaction between the test particle and a reservoir particle is a "hard-point" 
interaction (i.e. an infinite strength 5- interaction) . For the hard-point case, the reflection coefficient 
is £H(p) = 1. When the Dirac comb is set to zero, the distribution in momentum (p\p\j\p) converges 
exponentially fast in L 1 -norm to a Gaussian of width (^-) 5 - This is simple to prove, since it can 
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be reduced to showing exponential ergodicity for a classical Kolmogorov equation. I have discussed 
exponential dissipation in the article [11] for a three-dimensional case (hard-sphere interaction) for 
the purpose of studying diffusion. 

The discussion of the Hamiltonian H is the same as before except for the inclusion of physical 
constants. The Hamiltonian has a basis of kets \p) Q with energies given by E(p) = ^M^ip) f° r the 
anti-symmetric, increasing function q : R — > R determined by the relation 

,a N t a , \\ a M i a \\ ^ 

cos (~p) = cos (^q(p)) + j-tt sin (^q(p)j , P G M - — Z, 

and q(^) = ^ for n G Z. The dispersion relation .E(p) has jumps g n / at the values p = ^ for 
n G Z, which approach — as n goes to infinity. 

I will consider the model (|2,18|) in a limit in which the constants M, h, a, a and a time variable 
s G [0, to] scale with a single parameter < A <C 1 as follows 

M = m\- 1 , t = sX~ 1 , h = h\ 1+s , a = a X 2+g , a = a X 1+s , (2.20) 

for fixed constants m, h, otQ, ao, to, rj, and an exponent g > 0. I will add "A" as a subscript to 
the solution p\ t of the Lindblad equation (|2.18p and the escape rate £\{p) to indicate the parameter 
dependence. 

Let D\ t t(p) = c}{p\p\,t\p) q be the probability density for the momentum given by the diagonal of 
the density matrix p\ t in the extended-zone scheme representation. Define T>\ i G L X (R) to be the 
solution of the master equation 

= J R dp'(jx(p,p'PxAp) - Up\p)VxM) (2-21) 

with V\ )0 {p) = q{p\p\p)q, where the jump rates J\(p,p') are defined below. Let k\ :g (p,m) G C be the 
coefficients in the equation 

«A,ff(p, n)\p + q + n) Q = e [ ^p(2-\l + X)q + XP) * \p) Q . (2.22) 

The rates J\(p,p') have the form 

J\(p,p') = — ^ — y2 \P - P' - n\ \kx,p-p>- n (p ,n)\ 2 . (2.23) 

I will denote the Markovian process whose densities obey the Kolmogorov equation (I2.2ip as Kt, 
although the reader should note that the process has units of momentum rather than wave number. 

The following conjecture is analogous to Thm. 12. 1L The exponent g > from (I2.20P will only 
appear in the error of the following theorem. 

Conjecture 2.4. Freidlin-Wentzell/semi-classical limit. Let Dx^ and T>\ t t be defined as above. For 
A <C 1, then 

sup \\D X ,r-AAl = °( Xe )- 

Next, I state a conjecture analogous to Thm. 12.31 Let Pt be the Ornstein-Uhlenbeck process 
satisfying Po = and the Langevin equation 

dP t = - 7 P t + (^)5dB t , (2.24) 
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where is standard Brownian motion and 7 = 87/(^^)2. Define v = (^?) 2 ^7p- Conjecture 12.51 

states that the processes Aa|iT«,|, s £ [0, to] converge in law for small A to the absolute value of the 

3 — 

Ornstein-Uhlenbeck process above, and the normalized integrals X» J X drK r converge in law to a 
variable-rate diffusion. 

Conjecture 2.5. Let Kt be a Markov process whose probability densities obey the master equa- 
tion \2.21\) for a fixed A > 0. Define the time integral process Yt — ^ J drK r . In the limit A — > 0, 
there is convergence in law with respect to the Skorokhod metric 

(\*\K>\, xlfA (|P S |, -L- f dB' r |P r |§), ae[0,to], 
v A XJ v mvz Jo ' 

where Pf is the Ornstein- Uhlenbeck process {2.21$ , and B£ is standard Brownian motion independent 
ofPf 

If the Dirac comb is set to zero, then there is the standard limit result for Yt given by the conver- 
gence in law 

(**Kt, A§Yt) => (P t , — / drP r ), t G [0,t ]. 

JO 

The Markov process Kt without the Dirac comb has the same jump rates as the classical linear 
Boltzmann equation studied in [13] . The case a > is subdiffusive, since the spread in position is on 

3 1 

the order A~s rather than A~2 for A <C 1. 

2.3.1 Physical characteristics of the scaling regime 

The following lists the essential characteristics for the regime given by (|2.20p . The same mathematical 
mechanisms outlined in Sect. 12.231 should apply for this model, so I focus on a qualitative comparison 
of the relevant physical scales. 

1. The Brownian limit 

The scaling of the mass ratio as A = ^ while considering the dynamics over a time interval 
[0, growing proportionally to A" 1 is the standard regime for a Brownian limit. Since the 
temperature is fixed, the typical speed for a single particle from the reservoir and the test 
particle are (7^)2 and A^^j)^ respectively. Hence, the reservoir particles are moving faster 

than the test particle by a factor A~2 ^> 1. On the other hand, the momentum is (7^)2 for 

a single gas particle and A~ 2(^)2 for the test particle. Individual collisions with gas particles 
impart momenta that are much smaller than the typical momentum of the test particle. 

2. Frequency of phase oscillations versus collisions 

As I described in Section [2231 the autonomous evolution arising for the densities in the extended- 
zone scheme representation depend on the noise operating on a comparatively slow scale to the 
Hamiltonian dynamics. More precisely, certain phase oscillations driven by the Hamiltonian 
occur on a smaller time scale than the mean time between collisions. I characterize the relevant 
Hamiltonian-driven phase cancellations with the frequency 

liminf inf -r\E(p + n) — E(p) 

p-s-oo nGZ H 1 



:-lim inf g n 

h n— >oo 



¥ = » 1 

na nan 
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for A <C 1. I associate the frequency of collisions with the escape rates £\(p) for \p\ on the order 



;«)i = ( ? )ix-i 



Therefore, the phase oscillations occur on a shorter time scale than the collisions for small enough 
A. 

3. Kinetic energy outweighs potential energy 

The kinetic energy of the test particle is typically larger than the momentum stored in the 
potential by a factor of A -1 . For this comparison, I associate the typical potential energy with 
the strength a of the 5-potential divided by the period length: ^ = ^A. This is the potential 

energy, for instance, in a spatial wave (ft(x) = (2L)~z h—L,L]i x ) hi the limit L — > oo. The mean 
kinetic energy is simply 



4. Reciprocal lattice momenta and the reflection bands 



The half-spaced reciprocal lattice of momenta are multiples of ^ = ^2. This is on the same 

order in A as the typical momentum transfers for collisions ~ 2(^)5, and so the test particle's 
momentum has a chance of being kicked out of a given Bloch cell after several collisions. The 
reflection band around a lattice momentum — n for n E Z, where nonnegligible probabilities 
for Bragg reflections may be found, has a width of approximately = for high enough 

n so that 2s£ > ™°s> By a similar idea as in (|2.16p . the probability of a reflection when the 

particle's momentum is randomly dropped in the interval [ ^^a ' ~~^^~^\ ar °und p = ^p- is 
approximately 

2Ma 2ma p 2 a 

H\pJ = lM 2M*a- (2 ' 26) 

5. Bragg reflections are frequent over the time period [0, A _1 io] 

The frequency of Bragg reflections is equal to the frequency of collisions multiplied by a local 
averaged probability for reflection after a collision, which depends on the current momentum p. 
Multiplying (|2,25p with (|2.26[) . the effective frequency of reflections when the test particle has 
momentum \p\ 3> is 

^ } W\ = ¥\' 

Since the momentum will typically be found on the scale \p\ A~2 (7^)2 , a number of reflections 
on the order of A~2 will occur. 



3 Overview for proof of Theorem 12.31 

In this section, I will state the results that enter directly into the proof of Thm. 12.31 and then proceed 
with a presentation of the proof assuming those results. Thus, this section concerns only the classical 
Markovian process K r whose densities obey the Kolmogorov equation (|2.5p . 

Recall that £ r is the square root of the energy at time r £ M + : E^{K r ). In Sect. [5j I show that 
£ r is a submartingale. The following theorem states a central limit theorem for a rescaled version of 
c T . bmce E2(k)-\k\ is a bounded function, the convergence of t 2 £ s t, s £ [0, 1] to the absolute value 
of a Brownian motion as t — > oo is equivalent to the same statement for t _ 2|i^ st |. Thus, the first 
component for the convergence stated in Thm. 12.31 follows directly from Thm. 13.11 
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Theorem 3.1. In the limit t — > oo, the processes t~^£ s t, s G [0, 1] converge in law to the absolute value 
of a Brownian motion with diffusion constant a. Moreover, the martingale and predictable components 
M r , A r in the Doob-Meyer decomposition for E r have convergence in law for large t given by 



((at)-^M st , (aty^Ast) =A (b s , sup -B r Y s G [0,1], 



0<r<s 

where B is a standard Brownian motion. The convergences are with respect to the uniform metric. 

Theorem 13.11 only characterizes the behavior for quantities that depend on the absolute value of 
the momentum process K r , and thus the sign-flipping does not play a role. The concept of "sign- flips" 
is only meaningful when the particle has a high momentum from which a Levy jump to a momentum 
with the opposite sign would be unlikely. It is useful to define a series of stopping times that parse 
the time interval [0, t] into a series of excursions from the low momentum region. The following 
is a list of rough definitions for the notations related to sign-flipping and the time-integral process 
Y r = J'q dvK v . More precise definitions are given below, although the details for the definitions are 
not strictly necessary to understand the structure of the argument in the proof of Thm. 12.31 at the 
end of this section. The technical definition for the sign-flip times will not be quite adapted to the 
original filtration J>. 

rst 

Y s Normalized integral functional: Y® = t~^ / drK r 

Jo 

T m G R+ Time of mth sign-flip 

N r G N Number of r m up to time r G R+ 

3 

(zu n , ? n ) C R+ The nth incursion interval into the low momentum region ~ \K r \ < t» 
(? ra ,n7 n+ i) C R+ The nth excursion interval from the low momentum region 
T r G N Number of excursions to have begun by time r 

T r Information up to time r 

T r Information up to the time of the sign-flip following r 

hg Martingale with respect to J-^ := J- S t approximating t~^M s t for t > 1 

mi^ Martingale with respect to approximating as t S> 1 

Let S : R — > {±1} be the sign function. A sign-flip is said to occur at a Poisson time t n if 
S(Kt n ) = S(Kt n+1 ) and there are an odd number m of sign changes leading up to t n : S(Kt n _ r ) = 
—S(K tn _ r+1 ) for r G [l,m] and S(K tn _ m _ 1 ) = S(K tn _ m ). This definition avoids counting double-flips, 
which occur frequently in the dynamics and would distort the counting. The sign-flips are not hitting 
times, since information from the following Poisson time is required to identify them. Related matters 
are discussed in the beginning of Sect. 16.21 To define the time intervals where sign-flips will be counted, 
let the stopping times q j , Wj be given by so = w \ = and 



3 , , r , \ \ \ -w- r i 3 



'j = min{r G (?j_i,oo) | \K r \ < t& L }, = min{r G (wj, oo) | \K r \ > 2t& 



for some < i <C 1. The intervals [zuj, c^ ) are the incursions into low momentum for the process K r . 
Let T r G N be the number of excursions begun by time r: max{j G N | r > sv}. The times r m and 
increments Ar m are defined inductively for m G N such that 

• n = ft, 

• At to is the waiting time after r m such that either a sign-flip occurs or \K r \ jumps out of 
\MK I ^\K 11 

L2 \ Tm I' 2 \ Tm IJ ) 
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• T m+ i = T m + Ar m when r m + Ar m occurs during an excursion, and 

• T m+ i = $j when r m < Wj and r m + Ar m > ceTj. 

In most cases, r m are sign-flips and AT m = r m+ i — r m . Let N r be the number of T m 's to have occurred 
up to time r. 

I denote the standard filtration generated by the process K r with T r = cr(K s : < s < r). Let 
F r be the cr-algebra given by 



cj(Ar m , if s : < s < r m + Ar m ) r G [r m , r m + Ar m ), 

J" r V(m) : r ^ [r m , x m + Ar m ). 



When r G [T m ,r m + At to ) for some m, the cr-algebra T r includes knowledge of the time r m + Ar m 
and all information about the process K r up to time r m + Ar m . Since the sign-flip times are not 
necessarily adapted by the remark above, T r usually contains some information from the Poisson time 
following r m + Ar m to verify that the momentum does not change sign again. For J 7 ^ := J- S f, define 
the J 7 ]*) -adapted martingales 



m 



(*) 



:=t~i R r m (Ar m — E [Ar m | F Ty 



m=l 



:=t~^ ( £t + At ~£t -E[£ t +At -£ T \F - 



m=l 

The process £ r in the statement of Lem. 13.31 can be equivalently replaced by \K r \. The lemmas 
below all have the purpose of placing my convergence questions within the context of martingale 
theory. In particular, Part (2) of Lem. 13.31 is to establish independence for the copies of Brownian 
motion B and B' in the statement of Thm. 12.31 I will postulate the limiting law in a different way in 
the proof of Thm. 12.31 



Lemma 3.2 (Martingale approximations). In the limit t — > oo, there are the following convergences 
in probability: 



1. supo^i | h^ - t iM st \ 

2. supo^Jm^-yJ*)) = 



0, 



Lemma 3.3 (Quadratic variation approximations). In the limit t 
convergences in probability: 



oo, there are the following 



1- sup < s <i 

2. sup < s <i 

3. sup < s <i 



[ m (*),m(*)] s -I/ s dr(t-^ ri ) 



[mO, h(*)] 



0.. 



0. 



0. 



Lemma 3.4 (Lindberg conditions). The martingale m^*) satisfies the Lindberg condition such that 



as 



.(*)! 



t oo, then E[sup < s <i |mj — m 
integrable. The same statements hold for h^' . 



0. Also, the family mi*) for < s < 1, t G K+ is uniformly 
(t) 
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I will make frequent use of the reference book [31] in the proof below. 



Proof of Thm. \2.3[ All convergence in law will be with respect to the Skorokhod metric. I will show 
that there is convergence in law as t — > oo 

(|HiT st |, Y®) ^ (b s + sup -B r , m s ), s G [0, 1], 

where B is Brownian motion with diffusion rate a and m is a continuous martingale with quadratic 
variation processes satisfying 

1 f' s 

[m,m] s = - dr\B r \ 3 and [m,B] s = (3.1) 
v Jo 

for all s G [0, 1]. Recall that B s + sup 0<r<s — B r is equal in law to the absolute value of a Brown- 
ian motion with rate a. The above description of the limiting law is equivalent to the construction 
given in the statement of Thm. 12.31 By Part (2) of Lemma 13.21 and the fact that \E^{k) — \k\\ is 

bounded, I can approximate the tuple (\t~zK s t\, Y s ®) by (t~^£ s t, m^). Define the process quadru- 
ple Qs = (i~ aha ; raAjt, m s , [mW,mW],), where A r is the increasing part in the Doob-Meyer 
decomposition of £ r . The process Q s is adapted to the filtration F^ l \ and the first and third com- 
ponents are martingales with respect to Fs ■ 

By [SH Cor.VI.3.33], the family of processes t G R + is C-tight if the components are C- 

tight. The first component h" of Q® can be approximated for large t with t~?M s t by Part (1) 
of Lem. 13.21 The families of processes t~^M s t and t~2 A s t indexed by t G R+ are each C-tight, 
since they converge in law to continuous limits by Thm. 13.11 The fourth component of Q® can be 
approximated with ^ dr(t~z£ r t) 3 by Part (1) of Lem 13.31 The processes t~z£ st converge in law as 
t — t- oo to the absolute value of a Brownian motion by Thm, [3~7[] and the functional F : L°°([0, 1]) — > R 
defined by F(f)(s) = Jq drf 3 (r) is continuous with respect to the Skorokhod metric; therefore, the 

processes dr{t^^£ r t) 3 converge in law to J*q dr|B r | 3 for large t. It follows that [m^], t G R+ is a 
C-tight family. The family of martingales m' 4 ' must be tight by the C-tightness of [m^] and |31|, 
Thm. VI. 4. 13]. Finally, is C-tight by |31 1. Prop. VI. 3. 26] and the Lindberg condition in Lemma [3. 41 
Consider a sequence r n G R + such that C/i converges in law as n — > oo to a limit 



(1 f S 3\ 
B s , sup -B r , m s , - / dr(B r + sup -B„) ), sG [0,1]. 
0<r<s v JO 0<v<r ' 



The relationships between the first, second, and fourth components are determined by the consider- 
ations above and Thm. [37Ll By definition of C-tightness for Q^>, t G R+, the third component m 
has continuous trajectories. The second and fourth components are explicitly determined by the copy 
of Brownian motion B, so I will focus on determining the joint law of (B,m). Since 

(r ) (r ) _ 

is a sequence of martingales and the family of random variables , for n G N, s G [0, 1] is 

uniformly integrable by Lem. 13.41 the limit law is a martingale with respect to its own filtration [31, 

Prop. IX. 1.12]. The convergence (h.( Tn \ m^"-*) =^> (B,m) with |31|, Cor. VI. 6. 7] implies joint conver- 
gence with the quadratic variations 



h< r »), m^; 



[h( r »),h( r ™)], [h( r "),m( r ")] 
[m( r '"),h( r ")], [m( r "),m( r ")] 



B, m; 



[B,B], [B,m] 
[m,B], [m,m] 



I thus have that [m, m] s = ~ J* Q S dr (B r + sup 0<1)<r — Bt,) 3 . Finally, by Lem l3.3l the sequence of quadratic 
variation processes [m^ rn \ h( r ")] s converges to zero and hence [m, B] = 0. Therefore, the limiting law 
for subsequences has been determined uniquely as the one given by ()3. 1 j) . and the proof of convergence 
for (\t~2K s t\, Y s ) is complete. 

□ 
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4 The Freidlin-Wentzell/semi-classical limit 



In this section, I prove Thm. 12.11 First, I present Lem. 14.11 which collects some technical estimates 
on the probabilities \K v (k,n)\ 2 appearing in the jump rates (|2.6p for the limiting momentum process. 
The estimates pertain to large and include a technical restriction that the jump v is not too large 
(M < This is an important regime in future sections also, since the particle will stochastically 

accelerate to high momentum values, and the jumps v will be much smaller by the moment assumption 
(1) of List 12.21 on the jump rates j(v). There are only a few possible values of n for which \n v (k, n)\ 2 
is non- vanishing when |fc + u| 3> 1. 

As mentioned in the introduction and is discussed further in Appendix [Bl the Hilbert space L 2 (R) 
for the test particle has a canonical decomposition into a direct integral f r ® i n of copies of L 2 ( [— n, ir)) 
that are invariant under the Hamiltonian dynamics. It is useful to relate the kets \k) Q with their 
associated representations in L 2 ([— 7r,7r)). For k G R — iZ, the kets \k) Q are identified with Bloch 
functions tp^ G L 2 ([— 7r,7r)) given by 



2 



frW)-*'-! -tw(M , i27r(q(fc)-fe) ixq(fc) _ < < n 

1 _ e -i27r(q(k) + fc) C ' — 5 



(4.1) 



where > is a normalization constant, and q : R — )• R is determined by (12. 4j) . The form of the 
Bloch functions (|4.ip can be found in [TJ Sec. III. 2. 3] (under the standard quasimomentum and energy 
band labeling). The Bloch functions formally satisfy {x\k) Q = ipk(x), and tpk is an eigenvector with 
eigenvalue E{k) for the fiber Hamiltonian with k = ^modl. Analogously, the momentum kets \k) 
are identified with Bloch functions ipk G L 2 ([— 7r,7r)) defined by ipk{x) = (2ir)~ ? e lxk . 

The following notations are designed for a description of the reflection bands around the lattice 
momenta iZ. Let G : R — >• [—4, |) and n : R — >• Z be defined through 

8(Jfe) = fcmodi, n(fc) = 2(fc-0(ib)). (4.2) 

Note that the variables @(k) and n(fc) are noi the quasimomentum and energy band, respectively. 
Define the set I(k,v) C Z to include the integers 0, — n(A;), —n(k + v) and — n(k + v ) + n(fc). Also 

define (3(k) = |n(fc)e(fc), 7 (/?) = (£ 2 + - f , 



r_(fc) 



1 + 



|/3(fc)|+ 7 (/3(fc)) |2- 
|/3(fc)|- 7 (/3(fc)) I 



and v+(k) = 1 — r_(fc). The variable /3(fe) is a dilation of Q(k), which is scaled to characterize the 
limiting local profiles of reflection probabilities for high momenta k near lattice values |Z. 

In the following lemma, S : R — > {±1} is the sign function, and dist(n, I(k,v)) is the smallest 
distance between n G Z and an element in the set I(k,v). 

Lemma 4.1. Let \k\ S> 1. There exists a c > suc/i i/iai /or a// |&| large enough and \v\ < i/ie 
following inequalities hold: 

1. 

MM)| 2 < jjL. 
^ dist(n, I(k,v)) \K v (k,n)\ 2 < — . 
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< 



3. When n(fc) / n(fc + v), 

|k„(/c, — n(&)) | 2 — r_(fc)r + (/c + u) 

— n(fc + v)) \" — r + (k)r-(k + v) 

2 



fc|(l + |/3(fc)|) 
c 

< 



K v (k, —n(k + v) + n(k)) | — r_(A;)r_(A; + v) 
and when n(k) = n(k + v), 



\k\{l + \P{k + v)\) 
< 



\k\{l + \P(k)\) [l + \P(k + v)\) 



K v (k,0) 



rl(k)rl(k + v)- S(0{k)0(k + v))rt(k)ri(k + v) 
\ Kv (k, -n(Jfe)) | a - (r£(jfe)ri(Jfe + «) + S(e(Jfe)e(ife + u))ri(fc)r|.(A; + v) 

4- For meZ wii/j |y — A;| < and m ^ n(&), 



< 



< 



\k\ 2 ' 
c 

w 



i de\K^ +g _ k (k,-m)\ < — . 
Proof. 

Part (1): Recall that r](k,m) is defined so that \k) Q = X^mez ^(^' m )l^ + By mapping through 
the fiber decomposition, this is equivalent to the same equality with \k + m) and \k) Q replaced by 
ipk+m an d i>ki respectively. Moreover, the equation (J2.7P defining K v (k,n) is translated to 

^2 K v (k,n)tp k+v+n = e lvX ip k , 

where X is the bounded operator on L 2 ([— it, 7r)) acting as multiplication by the spatial variable 
(Xf)(x) = xf(x) for x G [— vr,7r). 

Below, I will show there exists a C > such that 



|^ fc - S(*0(fc)) rI(A;)Vfc-n(fc) " K(k)i; k \\ 2 < C\k\ 



(4.3) 



Temporarily assuming (|4.3|) , I will continue with the proof. Applying (|4.3p directly and for k replaced 
by k' := k — n(k), it follows that 



\i/j k , - S{k'e(k')) rl(k')^ kf _ n(h/) -4(^)^112 < 3C\k 



h-1 



(4.4) 



where I have used the identity k = k' — n(k'). 
Multiplying ip k by e 1 , then with 



JvX 



ip k = S(k®(k))r 2 _{k)i) k+V _ n{k) + rl(k)?p k+v + 0( 



where the error is with respect to the L 2 ([— 7r,7r)) norm. By the unitarity of e lvX , the norm of the 
error is preserved from (14. 3p . Applying (j4.4j) for k' = k -\-v and k' = k + v — n(k), 



- 1 ) =4(fc)4(fc + ^ fe+1) 

+ S(kQ(k + u))r^(fe)rl(fe + u)^fc+i>-n(jfc+v) 
+ S(ke{k))rl{k)rl(k + v- n(k))A +v _ n(k) 

- S(Q(k)Q(k + «))ri (fc)r'(fc + « - n(fc))^_ n(fc+u)+n(fc) , (4.5) 
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where I have used S(k + v — n(k)) = —S(k) and S(k) = S(k + v) by the restriction \v\ < %\k\. The 
minus sign on the fourth term on the right side appears because k — n(fc) has the opposite sign of k. 
It follows that 



y ] \ K v(k, 

m^I(k,k+v) 



m 



K v (k,n)ipk+ n = 0(\k\~ 

n£l(k,k+v) 



which will complete the proof once (|4.3p is established. 

Now I prove (|4.3p by examining ip^ in the Bloch representation (|4.ip . Note that (|4.3p is a little 
more precise than is strictly required to prove Part (1), but it will be useful later. Computing the 
coefficients r](k,n) yields 



Tj(k,n) = (ifjk+nl'lpk) 



Q 27ri(q(fc)-fc) 



1 



1 



+ 



1 



q(A) + k + n q(&) — k — n 



(4.6) 



It follows from (|2.4p that for large enough k, the inequality holds: \q(k) — k\ < ^m. The normalization 
Ni~ can be written as 



3 27ri(q(fc)-fc) 



4tt 2 1 + 



1 



ii 2 y| i + 

^— ' I q(fc) + A: + n q(fc) — k — n 
q(fc) - fc 2 



q(jfe) - fc + 2G(&) 



+ 0( 



All terms from the sum have been absorbed into the error 0(|fc| 2 ) except for n = 0, — n(k). The 
approximation above is due to the equality q(k) + k — n(k) = q(fe) — A; + 2@(k), and since 



E 

n^0,-n(fc) 



1 



+ 



1 



q(fc) - fc + 29(A) + (n + n(A;)) q(fe) -k-n 



< 2 sup 



(4.7) 



for |A;| large enough so that |q(fe) — A;| < |. When multiplied by |e 27rl ( q ( fc ) fc ) - l| 2 = 0(|/c| 2 ), 
then (|4.7p is 0(|A;|~ 2 ). Thus, the sum Yl n ^o -n(fc) l K f(^> n )P * s 0(|A;|~ 2 ), and the special terms r](k,0) 
and 77 f/c, — n(fc)) have approximations 



r/(fc,0) 
rj(k,-n(k)) 



q(fc)-fc+26(fc) 



q(fc)-fc 



1 + 



q(fc)-fc+26(fc) 
q(fc)-fc 

S(kQ(k)) 



- + o(\k\- 1 ) 



q(fc)-fc+20(fc) 
q(fc)-jfc 



(4.8) 



To complete the proof of (|4.3p . I require that the top and bottom expressions on the right side are 

rj(fc) + Odfcl- 1 ) and fl'(Jfee(Ar))r|(fc) + OG^- 1 ), respectively. 
By the definitions of ©(A;), n(A;), and (3(k), I can write 



q(A) — k 



q(fc) - A; + 2G(A;) _ (q(fc) - ^) + 9(fc) 



(q(Ar) - =£i) - 9(fc) 
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Define the variable 70 (k) = §(q(§ + 



2/3(fc) ■ 



) for n = n(k). The Kronig- Penney relation (j2.4j) to 



second-order for large \k\ yields that 70(A) satisfies 



p\k)= 1 l(k)-^ l0 {k) + O(\k\- 1 ). 

This gives the following asymptotics for 70 (A;): 

%- + 1 (p(k)) + 0(k- 1 ) S(k)Q(k)>0, 

--f(/3(k)) + odki- 1 ) s(k)e(k)<o. 



7o (*0 



I thus obtain 



1 + 



q(fc) - fc + 2G(fc) 



2\ -4 



q(fc) — A; 

where for the case S(k)@(k) > 0, I have used the identity 



r 2 _{k) 



0( 



1-1 



), 



(4.9) 





7(/3) + /3 




7(/3) - (3 



Therefore, combining (|4.9p with (|4.8p . the coefficient r/(fc, — n(fe)) is equal to S[kQ(k))r^(k)+0(\k\ 1 ). 

Part (2): By converting ifi^ to the momentum basis, operating with e mX , and translating back to the 
Bloch basis (ipk+v+n)nez, the coefficients n v (k,n) can be written as 



K v (k, n) = ??(& + v + n, m — n)r](k, m). 

From (|4.6p . I have the inequality 

|r/(p,r)| < 



(4.10) 



|q 2 (p) — p 2 



< 



q(p) + p + r I I q(p) — p — r | ' |q(p)+P + ^|| q(p) — p — r 

2 



(4.11) 



where /i := sup fegR q 2 (p)— p 2 and the first inequality follows by using the lower bound 4q ^^(p)-^^ — ~ 
for N p (i.e. from the single term n = from its sum). Recall that for large enough \p\, then 
|q(p) ~~ p| < ^rb- Roughly speaking, I will find the inequality (|4.11|) useful when \p\ is large enough 
so |q(p) —p\<\ and r is not or —n(p). In that case, the terms in the denominator have the lower 
bounds |q(p) -\-p + r\, |q(p) — p — r\ > 4. When (|4.1ip is not of use, I still have the trivial bound 
\r]ip,r)\ < 1. 

For notational ease, I will restrict to the case n(/c) 7^ n(/c + v) in this proof. Let me first focus on 
the values n such that ^\k + v + n] -1 < j. Define the constants a,j 6 R through 



a>j = dj(k, v, n) 



q(A) - k j = 1, 

-q(A;) - k j = 2, 

q(/c + w + n) — A; — v j = 3, 

— q(fc + w + n) — k — v j = 4. 



I will use the inequality (|4.1ip to bound the sum of the terms m from (|4.10p with the special values 
m = 0, — n(fc), n, —n(k + v) — n removed: 



1 

X 2 



|r/(fc + v + n, m — n)r](k, m)\ < 



m^O,— n(fc), 
n, — n(fc+t>) — n 



m^O,— n(fc), 
n, — n(fc+u)— n 



to — 01 m — 02 m — 83 m — 04 



-. (4.12) 
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Let d > be the largest radius for intervals centered at integer points such that the intervals never 
contain more than two of the a'jS: 

d = d(k, v, n) = sup {s G R + | V(m G Z) : |{di, a2, 03, 04} D [m — s, m + s]| < 2}. 

If no three elements of {ai, 02, 03, 04} are within a radius d > from any single integer, then the 
sum (|4.12p has the following bound 

-1 V" V 1 8/ f°° 1 \ _ 40 

d? ^ ^ \m-aA 2 < d 2 V + i +x )2j " d 2 ' 

j=l,2,3,4 m^O-n(fc), 1 - 71 17 U V4^^ 

n, —n(k+v)—n 

The inequality is a Riemann upper bound using that the distance of the a^-'s to an integer not equal to 
0, — n(/c), n, or — n(k+v ) — n is > \. I claim that the maximal radius d(fc, f , n) increases proportionally 
to \k\ V |n|. Within the nearest integer, I have a\ ~ 0, 02 ~ — 2fc, 02 ~ n, and 04 ~ —2{k + v) — n. For 
77. large enough (i.e. |n| > 6|A;|), then the claim is clearly true, since \v \ < ^\k\. For |n| on the order of 
\k\ or smaller, I already have that and —2k are far apart. In order to have three of the Oj's within a 
small radius, then both n and — 2{k + v) — n must be near or —2k. However, n and — 2{k + v) — n 
can only reach within a distance < j\k\ of each other if both are at least a distance > ||/c| from and 
—2k (again by the constraint \v\ < ^I^D- 

For the finitely many values of n such that ^|/c+w+n| _1 > 4, I can give almost the same treatment. 
For the integers m closest to — q(/c + v + n) — k — v and q(/c + v + n) — k — v, I simply use the bound 
\r](k + v + n, m)\ < 1. The m's closest to those values must be a distance > ^\k\ from either or —2k, 
so \r](k, m)\ < The remainder of the terms can be treated as in the case above. In both cases, the 

sums were bounded by a constant multiple of (\k\ V |n|)~ 2 . 

Now, I deal with the four exceptional terms m = 0, —n(k), n, —n — n{k + v), where exactly one 
the terms rj(k + v + n,m — n) or rj(k, m) is of the form r/(p, r) for r = 0, —n(p). For that term, I use 
the inequality \rj(p, r)\ < 1 rather than (I4.11|) . For the other term, I apply (14. lip and approximate 
2{k + v) ~ n(k + v) and 2k ~ n(fc) (and I double the constant h to cover the error resulting from the 
approximation). The term \rj(k + v + n, m — n)r](k, m)\ is smaller than 

2/i(|n| |n + n(k + 1 m = 0, 

2/i(|n + n(/c)||n + n(/c + — n(/c)|) 1 m = — n(A;), 

m = n, 

2/i(|n + n(/c + u)||n + n(/c + v) — n(/c)|) m = — n — n(k + v). 



All of the four terms are bounded by 2h times 

(\n\ A |n — n(A;)|) _1 (|n + n(&)| A |n + n(k + u)!)" 1 . (4.13) 

For values of n near I(k,v), this will be larger than the expression (\k\ V |?t|)~ 2 that bounds the 
remainder of the terms not in I(k,v), although for larger values of n, they have the same order. Thus, 
the values \n v (k,n)\ are bounded by some multiple h! of f)4. 13|) . 
I now have that 

dist(n, I{k, v)) \K v {k, n)\ 2 

n0(k,v) 

<ti V] dist(n, I(k,v)) (|n| A |n — n(fe) + n(A; + u)|)~ 2 (|n + n(A;)| A |n + n(A; + -u)|) _2 

n0(k,v) 

<ti ^2 (|n + n(A;)| A \n + n(k + v)\)~ 2 + ti ^ (|n| A |n - n(/c) + n(fe + v)\) ~ 2 . 

w£I(k,v) n£I(k,v) 
4S(k)n>-\k\ 4S(k)n<-\k\ 
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Since \v\ < ||&|, the above is bounded by the following 

A, ( i| fc | +m) 2- | fc | | fc | A, (1 + g )2 - | fc |7 x { i + x f- \ k y 

where the inequality is by a Riemann upper bound. 

Part (3): The analysis from Part (1) gives errors for \n v {k,rn)\ 2 , m £ I(k,v) of order 0{\k\~ l ), but 
there is additional decay when (3(k) or /3(k + u) are large. All the cases involve similar reasoning, so 
I focus on the case for \n v (k, —n(k + v)) \ 2 with n(A; + v) ^ n(k). If there is a C" such that 

I I C 

\n v (k, -n(k + «)) - S((k + u)0(Jfe + u))ri(fc + u)r£(fc)| < m (4.14) 



|fc| 



for some C" > 0, then 



2 2C" i i 

\K v (k, — n(k + v))\ — r-(k + v)r + (k) <- — -r^.(fc + v)r? (A;) 

IV |fc| 

^(i + l^fc + f)!)" 1 . 

The second inequality is for some C" > 0, since r_(/c) is bounded by a constant multiple of (1 + 
|/3(/c)|)~ 2 . Thus, establishing (|4.14p is sufficient for the proof. 
By (|4.10p and the triangle inequality, 

\K v (k, —n(k + v)) — rj(k + v — n(k + v),n(k + v))r](k,0)\ 

< \fj\k + v — n(k + v), —n(k) + n(& + v))r)(k, — n(fc)) | 
+ \fj(k + v — n(k + v), 0)77 (A;, — n(& + u)) | 

+ | ry(fc + v — n(A; + w), m + n(A; + v)}r](k, m)\- 

m^I(k,k+v) 

By (|4.3p . Sn^o -n(fc) M^j 71 )! 2 = 0{k~ 2 ) and with Cauchy-Schwarz the last term above is 0(\k\~ 2 ). 
Moreover, by the bounds in Part (2), the first two terms are each bounded by a multiple d > of 
For large enough \k\, I have 

3c' 

\K v (k, -n(k + v)) - rj(k + v — n(k + v), n(k + v))rj(k, 0)| < -rj-r. 
By the analysis in Part (1), there is a C > such that 

C 



|r?(fe + v - n(k + v),n(k + v))rj(k,0) - S((k + v)Q(k + v))r'i(k + v)rl(k)\ 
Putting the above inequalities together, then I obtain (14.141) for C = C + 3c'. 

Part (4): By the bound in Part (3) for I n v (k, — n(fc + v )) | 2 when v = y + — k, 



~ 1*1' 



1 . 1 



y_i 1 2+ 1 -y_i v w v 2 ' ifci(i+i/3(f 



4 

2 r 00 , 1 c 



~ IH /-oo ^ 1 + I M + ?'[ U '! I 2 + 2|fc| 



u)|— 7(ui) I 



where the order equality uses that \m\ > \k\. 



□ 
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Let T : L 1 (]R) be the trace preserving map with integral kernel T(k\,k2) = TlT 1 J{k\, £2)- The 
idealized momentum process (|2.5p has a pseudo Poisson form in which jump times are determined by 
an outside Poisson clock, and the jump transition densities are given by the operator T. A similar 
structure holds for the original Lindblad dynamics by (2) of Lem. IA.H 



Proof of Thm. \2.1\ Define the map p — > [p] D from i3i(L 2 (R)) to L X (]R), which sends a density matrix 
to its diagonal density in the extended-zone scheme representation [/o] D (&) := q(&M&)q. The diagonal 
map is well-defined by the discussion in Appendix lB.21 Let £ = (fi, • • • ; ttf) be the sequence of Poisson 



times less than t. Define the map : B\\L (K)) as 



■Af - i(t ~ tAf } H 



The maps <3?^ are completely positive and preserve trace for all £ and i. Similarly to the construc- 
tion (|2.8p of the solution to the Lindblad dynamics, I have that px t t = E[3>^(p)] by (2) of Lemma lA.H 
where E is the expectation with respect to the Poisson process with rate 1Z for the sequences £. Also 
T>t = E[T^*'^[p] d ] , where J\f(t, £) is the value of the Poisson process at time t. 
I can write 



^A lt -^ = E[[$g(p)] D -^)[4 



E 



[ e ^ ) - n [<^(p)] D -^ ) - n+i [< A L,i(p)] D ; 



n=l 



where I have inserted a telescoping sum. Since T is contractive in the 1-norm, is contractive in 

ir rr 

trace norm, and conjugation by e~~ is contractive in trace norm, I have the second inequality below 



Wxt-Vth <e 



-TZt 



00 A/" 

X>^E 

AT=1 n=l • /n ' 'I'-'.V ' 

00 AA 

^ E^E 

_A/ = 1 ■^0<tl---<i r! ._l<tn+l<---t^/-<t 



i/V-n+1 r^( A ) 



x sup 

l[p]ll=l 



dt n (V 1 [tf 



e a H pe a H ) 



T[ P ]u 



E 

AT=1 



(A/--1) 



sup sup 

>lll=l 



dc (n- 1 [^(e-'f H pe i i H )} u -T[p] D 



where I identify to = and ijv+i = t for the boundary terms. The first inequality above is the triangle 
inequality. With the above and the summation formula 



n + n 2 t = e~ nt 



it follows that 



\D\,t ~ A||i < {Tl + K 2 t) sup sup 

a,b&R+ ||p||i=l 



(4.15) 
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The remainder of the proof is concerned with proving that the supremum in (14.15P is bounded by 
a multiple of A for A <C 1. By a direct calculation, 

"dc {II- 1 [¥(e-x H P ef H )} D (k) - T[p] D (k)) 
r j(v) , 

= — iA / dv —— 2 K v (k — n — v, n)K v (k — m — v, m)p(k — n — v,k — m — v) 

n=£rn 

e -f(E(k-n-v)-E(k-m-v)) _ (E(k-n-v)-E(k-m-v)) 

E(k — n — v) — E(k — m — v) ^ 



Using the triangle inequality and making a change of variables k — n — v — >■ k and M = n — m, I have 



< A yj t dk MM)I M* ± m. n - M)l 

- L L 7e |J5(jfe) - + Af) 



,2vr 



< A (-+^C M ), (4.17) 



where the values Cm are defined as 

n f, j(v) Zn\ K v(k,n)\\ Kv (k + M,n-M)\ 

Cm = sup / dv 



\2k+M\>Jm n \E(k) - E(k + M)| 

The second inequality in (|4.17p can be found by splitting the integration J K dk into the regions 1 2k + 
M\ < 1 and |2fc + M| > 1. This splitting isolates some bad behavior (non-decay for large M) occurring 
in regions of k where \2k + M\ is small. The Cm's in (|4.17p arise by applying Holder's inequality over 
the integration \2k + M\ > 1 and by the inequality \p(ki, kz)\ < \p{k\, k\) + \p{k%, fe) for the integral 
kernel of p. The kernel inequality follows because p is a positive operator. The — term in (|4.17p 
comes from the \2k + M\ < 1 integration for which I apply mainly brute force: 

V f dk t dv 3{v) Z n K(k,n)\\ Kv (k + M,n-M)\ 

^Wi i n \E{k)-E{k + M)\ \P^ k + M )\ 

< ( sup \E(k) - E(k + M)!" 1 ) V f dk (lp(k, k) + \p{k + M,k + M)) 

k,M M _£ Q J\2k+M\<l 2 2 

f 2ir 

< (inf^)-^ / dkp(k,k) = —, 

where g n > are the gaps between the energy bands occurring at momenta k £ — {0}. The 
infemum of the energy gaps g n is — . The sum over n £ Z of |/Cu(fc, n)| |/Cu(A; + M, n — M)\ is bounded 
through the Cauchy-Schwarz inequality and ^ |K„(/c,n)| 2 = 1. The key observation is k and + M 
must lie on different energy bands, and it follows that E(k) and E{k + M) differ by at least the length 
of the smallest energy band gap. 

Next, I need to show that the sum of the Cm's is finite. A single Cm can be bounded by (inf n g n )~ l 
using some of the same reasoning as above. I will show Cm decays on the order of \M\~ 2, and thus 
is a summable series. The difference \E{k) — E(k + M)\ necessarily becomes large for \M\ 3> 1 except 
for cases when k + M is close to —k. However, by my restriction \2k + M\ > 1, the momenta k and 
k + M will not lie on neighboring energy bands, and thus their energies must differ by at least the 
length L\m\ of the |M|th energy band. By [TJ Thm.2.3.3], L n grows with linear order for n>l, Also, 
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if A \k + M\ < l\M\, then \E(k) - E{k + M)| > E -\M\<n<\M\ Ln wm grow on quadratic order in 
0(\M\). 

I have Cm < L M C' M + C M , where C' N and C'^ are defined as 

C' M = sup f dv ^ET y2\K v (k,n)\\K v (k + M,n- M)\, 

\2k+M\>l JR rL n 
|fc|A|fc+Af|>^|M| 

En MM) I \n v (k + M,n- M)\ 



C ' L ~ m *Tm<\\M\k dV: * \E(k)-E(k + M)\ 

By the observation above, C M < sup| fe | A | fc+M | < i M \E(k) — E(k + M)\~ l decays quadratically. The C'^ 
are therefore summable, and I turn to the only somewhat delicate part of the proof, which requires 
isolating the problematic terms in the sum of the \n v {k, n) \ \ n v (k + M, n — M)\ contributing to C' M to 
which I can apply Lem. 14.11 

For fixed k, v, the Cauchy-Schwarz inequality and En I 

)| = 1 yield 

^2\K v (k,n)\\K v (k + M,n- M)\ < ^ \K v (k, n)| \n v {k + M, n - M)\ 

n n£l(k,v) 

n€l(k+M,v)+M 

+ ( MM)| 2 )* + ( Yl K(k + M,n)\ 2 y. (4.18) 

n£I(k,v) n0(k+M,v)+M 

Under the constraint \k\ A \k + M\ > \\M\ and by Part (1) of Lem. 14.11 the terms on the bottom 
line of (14.180 decay on the order \k\~^ < 2|M| _ 2 an d \k + M|~2 < 2\M\~^, respectively. The 
weighted integration J K dv is finite so these terms make contributions to C' M that vanish with 
order 0(|M|~i). 

Controlling the integral J" R dv of the first term on the right side in (|4.18|) will now require 
invoking the decay of j(v) at infinity and its boundedness through Lem. 14. 1[ The constraints n £ 
I(k, v), n£ I(k + M, v) + M, \2k + M\ > 1, and \k\ A |A; + M| > \\M\ leave the following possibilities: 

• The inequality |/c + u + ^M| < | holds and either n = or n = M holds. 

• The inequality \v — \M\ < j holds and either \k + < j or \k + v + < j holds. 
The second case vanishes, since \v\ ~ o|M| » 1 and by Jensen's inequality: 



/ 

Jv 



dt>4rEMM)IM* + ^"-M)|< / dv J -^<^\M\-\ 
„ J\v\>Im 



<\v\>\M '<■ ~ J\v\>\M 

where a = L dvj(v)v 2 . The first case follows by Part (4) of Lem. 14.11 



□ 



5 Submartingales related to energy 

In this section, I discuss certain key submartingales appearing in both the quantum and the limiting 
classical settings. First, let me define an operator-valued submartingale. Consider a probability space 
P) with a filtration Ft, a Hilbert space H, and an operator- valued process Y t : 0, — > B(H) 
adapted to T% and satisfying 

E[|(/|y t /)|] <oo, fen. (5.1) 
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I call Yt a submartingale if E[lj LfJ — Y s is a positive operator for all t > s. Naturally, Y% is a 
martingale if both Yj and —1^ are submartingales. I can extend my definition to the case in which Yt 
may take values as an unbounded operator. In this case, I require that there is a single dense space 
D C % such that domain of Yt almost surely contains D for all t and (|5.1|) holds for all / € D C %. 
In the unbounded case, I refer to the process by the tuple (Yt,D). 

For the following discussion, I will consider a Schrodinger Hamiltonian H = P 2 + V(X) with 
positive potential V and domain T)(H) C L 2 (R). The noise in my model is generated by an underlying 
Levy process L t with jump rate density j(v). As before, let v n £l,t„ £ M + denote the jumps and 
jump-times for the Levy process, and Mt be the Poisson counter for the jump-times. Also let the 
unitaries U\^{C) be defined as in (12. 9p . Define the operator- valued process Gt as 

G t = G x ,t(0 = Kt(OGUx, t (0 (5.2) 

for an observable G acting on L 2 (M). On bounded observables G, the trajectories are right weak*- 
continuous with weak*-limits existing from the left, since the Hamiltonian evolution in the Heisenberg 
representation is weak*-continuous. Analogously to (|2.8p . the Heisenberg evolution for G G S(L 2 (R)) 
can be written as 

$^(G)=E[C Ait (0]. (5.3) 

In Prop. I5.1| the formula (15. 3ft may be interpreted as the definition for the dynamical maps t : 

£>(L 2 (M)). I will suppress the A and £ dependence for the operator processes in the future. 

In the lemma below, I study ()5.2[) for the special cases G = H and G = H?, and prove that 
i j 

(Ht,D(H)) and (-ff t 2 ,D(i7a)) are operator submartingales. I show that each process is a sum of 

a martingale and an increasing part, and I have presented the increasing part of the Doob- Meyer 
i 

decomposition for H t 2 in a form that is not predictable, but which will be useful later. The linear 
spaces D(if) and D(H^) are closed under the operation of the unitaries UxjiQ, since the evolution 
e a clearly leaves the domains invariant, and 

e -ivX H \ e wX _ R \ and e -ivX H( ,ivX _ R 

are relatively bounded to H 2 and H, respectively. These relative bounds can be shown using the 
Wigner- Weyl identity e~ ivX He ivX = (P+v) 2 + V{X) and, in the case of H 2 , a resolvent representation 
for the square root of an operator (see (15. 5|) ). 

Note that I have used the Wigner-Weyl relation e~ ivX H^e ivX = ((P + v) 2 + V{X)Y in the 
statement of Part (2) of the proposition below. 

1 

Proposition 5.1. Let the operator-valued processes Ht and be defined as in A5.2\) for a Schrodinger 

1 x 

operator H = P 2 + V(X) with nonnegative V and domain T)(H). The processes (H t 2 ,D(i^2)) 

1 -l 

and (Ht ,D(ff2)) are submartingales that can be written as a sum of martingale parts (M^, D(i?)) , 
(M' t ,L 2 (R)) and increasing parts (A t ,D(if)), (A' t ,D(H^)) , respectively, with the forms found below. 

1 . The submartingale Ht is equal to 1VL + A^ for 

Aft Aft 
M t := 2y2v n P t - + y^v^ - at and A t := H + at. 

n=l n=l 
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1 

2. The submartingale H£ is equal to + A' t for 

Aft 



M t ■= \Y.^ P t- n + v n? + V{X t -jY-{{P rn -v n f + V(X r J)K 

n=l 

1 Nt i i 

a; := m + -J2((P t - + vn) 2 + v(x t -)Y + {(P t --v n ) 2 + v(x t -)Y 

n=l 

-2{pi+v{x t -)y. 



3. Moreover, I have the operator relations 

®* Xtt (H) =E[H t ] = H + at and $\ t {H 2 ) = E [H 2 ] < H 2 + 3aHt + + ^cr 2 t 2 , 
where ? := j R dv v 4 . 

Proof. 

Part (1): The energy process Ht can be written as 

H t = e^ H e-^ x ■ ■ ■ e -^ x e^ lH He~-^ JllH e iv - x ■ ■ ■ e [ ^ x e-^ H 
= H t - + 2v n P t - + v 2 n , 

for n = Nt- Through iteration of the above calculation, I obtain the relation Ht = + A^. By 
the symmetry of the rates j(v) = j(—v), it is clear that M 4 is a martingale, since Ylm=i v n — ta is & 
martingale. 

For / £ D(H), I will now verify condition ()5.ip for Ht and Mj. For Ht, I have 

E[(f\Ht\f)]=(f\H\f) + ta\\f\g 

since the martingale part vanishes under the expectation. By two applications of Jensen's inequality 
for the first inequality below and using that (f\P 2 \f) < {f\H t \f} for the last inequality, 



E[|</|M t |/>|] 2 <||/||2E[</|M t 2 |/>] 

= 4||/||iE[f;, 2 </|P 2 _|/) 

n=l 

= 4a||/|| 2 , f drE[(f\P 2 \f)} 
Jo 

<4a||/|| 2 (t</| J f/|/> + ^||/|| 2 ). (5.4) 



The domain for At is clear from its form. 

i 

Part (2): The equality Hf = M.' t + A' t can be shown through a telescoping sum and the conservation 
of energy between momentum kicks in a similar way to the argument in Part (1). Also, that M.' t is 
formally a martingale is clear through the symmetry of the rates j(v) = j(—v), however unlike for 
Part (1), it is not immediately clear that A' t is an increasing process. 
I have the equality 



e ivX Hh~ ivX + e ~ ivX Hh ivX -2H s = [" da f db e ibX i [X, i [X, m] ] e~ ibX 

J-\v\ JO 
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by the second-order Taylor expansion 



e lvX m e 



-ivX 



H* +iv[X,H*] + [ da I dbe ibX i[X,i[X,H^]]e 



ibX 



Showing that i[X, i[X, H 2]] is a positive operator will complete the proof that Aj is increasing. 

Using the identity 02 = i f °° dww~%-^-^ and functional calculus, the operator can be repre- 
sented through its resolvents (see [371 Ch.VIII, Ex.50]) as 



7T 



(feu u; 2 . 



H 



H + w 



(5.5) 



Evaluating both sides by the double commutator with X, 



i[X,i[X,H*]] 



X, i 



x, 1 - 

7T 

dw u)2 



aw w 2 . 



H + w 



(H+ W y 

since through canonical commutation relations I obtain 



iJ—pJ—p- 1 



'H + w H + w H + w)' 



(5.6) 



X, i 



X, 



H 



H + w 



2w 



(H + w) 2 H + w H + w H + wJ' 



The following operator inequalities hold: 



1 -p^-p^- < ^-pi P - 1 



H + w H + w H + w H + w P 2 + w H + w 

1 / H x 1 



< 



■(— )■ 



H + w y H + w J H + w' 



where I have used that f(u) = and f(u) = -^-^ are operator monotonically decreasing and 
increasing functions, respectively. Applying this inequality to (|5.6p . 



i[x,i[x,m]] 



> 



7T 



dww^ 



H 



(H + w) 2 [H + wf 



0. 



To see the equality to zero, I compute the two integrals through a change of variables w = /itan 2 (#) 
and find 



1 

U72 



7T Jo ^ (h + W) 2 

aw ■ 



h^TT Jo 
2 



n Jo (h + wy 



h^TT Jo 



d6 sin 2 (6») 



dd sin 2 (2(9) 



1 



, 1 ' 
2/12 



2/tl ' 



Hence, i[X, i[X, -ff 2 ]] is a positive operator, and A£ is increasing. 

Next, I show that satisfies condition (|5.ip . Again by Taylor's formula, 



e wX Hie 
By similar reasoning as above, 



-ivX 



\nX H \ e wX = I dbe ibX i[X,H^]e~ 

J —V 



ibX 



00 , 1 1 „ 1 

dww' 2 — P- 



H + w H + w 
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Using that \P\ < H 2 , I have the inequality 



[f\i[x,m]f)\ < 



7T 



dlV ( / 



(H + w) 2 



f 



where I have used the same change of integration as in the functional calculus above. Since e lbX has 
operator norm one, the above remarks imply 



ivX H2e~ [vX 



< 2\v 



(5.7) 



Finally, for / el 2 

e[|(/|m;/)|] =e 



dr / dvj(v) 



e ivX Hr 2 e~ ivX 



' ivX H r 2 e ivX 



<2t\\ff 2 / dv\v\j(v)<2t\\fgK2a2. 

JR 

1 

The same calculation shows that H t 2 satisfies condition (4.1). That condition also holds for A' t , since 

1 

the process A' t is the difference between and M.' t . 

Part (3): Since the martingale part has expectation zero, the equality 

E[H t ] =H + at (5.8) 
follows trivially from the decomposition in Part (1). For E[i7 t 2 ], the classical theory would have 



H 2 + (M, M) t + 2 / dA r H,, 



E[Hf] =M[(M t + A t ) 2 } =E 



where (M, M)t is the predictable quadratic variation. In my case, dAt = crdt is a multiple of the 
identity operator and therefore commutes with everything. It follows that the equality above holds by 
the same argument as for the classical case. The processes Xm=i Vn -> En^i [ v n ~ > an d ^-M ~~ °^ 
are uncorrelated martingales, and thus 



Mr 



Mr- 



Mr 



<M, M>, =( YJ F C „„, £ F t „„) ( + ( £ R - i) , YJ 



Av 



n=l 



n=l 



71=1 



71=1 



(7 



K) It 



+ (n K - ar >n K - ar n 



(5.9) 



Using P 2 < il r and (JEHJ gives the bound for E[(M, M) t ]. Bounding E[f*drH r ] also follows 
from (1S~%1) . 

□ 

In the proof of Prop. EH I apply Prop.[5j]to gain information about certain martingales related to 
the Markov process K r . Define the energy process E r := E(K r ), where E(k) is the dispersion relation 
determined by (|2.4p . Define £ : R — > R+ as the square root of the energy: £ (k) = E^(k). I also 
use the symbol "£" to refer to the corresponding process £ r = £{K r ). Recall that the kets \k) Q are 
associated though a fiber decomposition of L 2 (R) with normalized Bloch functions ipk 6 L 2 ([— n, 71")) 
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given by (14. ip . The mathematical connection between the results in Prop. ET] and the classical process 
K r is made through formulae such as in the equality 

Q ( k\ e~ ivX H h ivX | k) Q = ^ | *v {K i) 1 2 £{k + v + i) 

=E [Er | K r - = k, dL r = v] , (5.10) 

where L r is the underlying Levy process. The first equality above follows from the definition of 
the coefficients \n v {k, i\\ . The value G [0,1] is the probability for a lattice jump i G Z 

conditioned on a Levy jump v occurring from the momentum k. The rigorous meaning of expressions 
involving bra-ket notation can be traced back to the fiber decomposition such as in 

Q (k\e- ivX Hh ivX \k) Q = (e ivX ^ k \Hl e ivX $ k ), 

where H$ is the fiber Hamiltonian for (f> G [— ^ , ^ ) with A; + v = </>mod 1. 
Define the process m r as 

1 ^ / 2 2 

"V := 2 £ Kn(^in-^)| 2 £(^"tn- + + ») ~ ^ \ K -V n ( K t n -, ^ ^ {K tn - ~ V n + j 

n=l ieZ J 

By the rate symmetry j(-u) = j(—v) of the Levy jumps, m< is a martingale. 
Proposition 5.2. 



1. T/ie process E r is a submartingale, and the predictable increasing part of its Doob-Meyer decom- 
position is or . Moreover, the second moment satisfies 



E[E 2 ] < E[E$] + 3arE[E ] + TZ^r + -a 2 r 2 , 



where ? := L dv 



n 



V 



2. The process £ r is a submartingale, and the martingale part M r of its Doob-Meyer decomposition 

777 { 



has predictable quadratic variation (M,M) r with 4-{M,M) r < a. 



3. The martingales m r and M r — m r are uncorrelated and therefore 4-(M,M) r > 4-(m,m) r . Also, 
the quadratic variation of m r satisfies [m,m] r — [m,m] s < Yl~n>M a v ri f or r — s > w hich means 
that it is dominated by the quadratic variation of the Levy process Lt- 

4- The derivative of the predictable quadratic variation for the martingale m r (i.e. 4p(Tn,m) r ) is a 
function V : M. — > [0,a] of K r . There exists an a > such that for all ligl, then 



Proof. 



Part (1): Let u G L 1 (M) be a probability density with f^dku(k)E 2 (k) < oo. Construct the density 

matrix p = \ f){f \ G £i(L 2 (M)) for f{k) = u3(jfe). Let D x , t , p\,u and V t be defined as in Thm. ED 
For all A > 0, 

/ dkD Xit {k)E(k) =TT[p^ t (P 2 + V(x))] 
=Tr\p$lt{P 2 + V(x))] 

= / dku(k)E(k) +at. (5.12) 
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The third equality is by Part (3) of Prop. 15.11 

By Thm. 12.11 D\j converges to T>t hi the 1-norm as A — > for every hxed t. To guarantee the 
convergence of J R dkD\j(k)E(k) to j R dkV t (k)E(k), it is sufficient to have a uniform bound in A 
on the second moments f R dk D\ :t (k)E 2 (k). Again using Part (3) of Prop. |5TT| there are constants 
ci(t),C2(t) G R + such that 

n,t{( p2 + V ( X )f) = E i( P t + y( x t)) 2 ] <H 2 + c x {t)H + c 2 (t). (5.13) 
Hence, a uniform bound for the second moment of the energy is given by 

/ dkD Xit (k)E 2 (k) < [ dku(k)E 2 (k) + d(t) [ dku(k)E(k) + c 2 (t). (5.14) 
ii Jr Jr 

It follows that E u [i^] = J^dkT>t(k)E(k) is finite and equal to (|5.12p . where E u is the expectation 
beginning from an initial distribution u. The jump rate densities jfc(fc') = J(k',k) are continuous in 
L X (IR) as a function of k over intervals between lattice points in ^Z. I can approximate a 5-distribution 
at k G M — with densities, and I have the result ^E^fEy |( = o = a for k not on the half-spaced lattice. 
Since my initial distribution will always be a density and the jump rates are densities, the behavior 
assigned to the lattice values is irrelevant. It follows that Et is a submartingale with increasing part 
E + at. 

The bound for E[Eyl follows by plugging in the explicit values for c\{r) and C2(r). 

Part (2): Let u G L 1 (IR) be defined as in Part (1). By Part (2) of Prop. EZQ for every t, A > the 
inequality below holds 

/ dkD x . t (k)E 1 2(k) =Tv[p Xjt (P 2 + V)*] 
Jr 

=Tv[p^l t ((P 2 + V)^] 
>Tv[p(P 2 + V)?]. 

i 

A similar argument as in Part (1) shows that is a submartingale. 

By Part (1), the increasing part for the Doob-Meyer decomposition for 6 2 = Et increases with 
linear rate a. I thus have the relation 

a=^(M,M) t + 2£ t ^-A t . 
dt dt 

Since both terms on the right are positive, it follows that 4t(M,M)% < a. 
Part (3): By (|5.10p . the terms of m r can be rewritten 

J2 (Mfc, *) \ 2 £(k + « + »)- \^v{k, i) \ 2 £{k -v + i))= Q (k\e-' wX me lvX - e wX e~' mX \k) Q . 

By (15 .7p . the absolute values for the jumps of mt are bounded by the absolute values for the jumps 
v of the Levy process. Consequently, the increments for the quadratic variation [m, m]t — [m, m} r are 
almost surely smaller than those of the Levy process: [L, L]t — [L, L] r for all t > r > 0. 
The process Et can be written as 

Mt Mt 
£t = £o + XX " £ t- n ) = 5 * + E ( £ t« ~ £ t- ~ E K - £ t~ \^t- . dL t n ] ) , (5-15) 

n=l n=l 
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where St = Eq + Xm=i ^[^t n — £+- \^t~ > dL tn ] . For a single term in the sum and K t - = k, dL x 
then 

E[£ tn - £ t - \T t -, dL tn ] = K(k,j)\ 2 £{k + j + v) - £{k) 

= Q (k\e- ivX Hh ivX \k) Q - Q (k\'H 1 i\k) Q . 

However, I can reorganize St in a way reminiscent of Part (2) in Prop. 15.11 

-. Aft 

S t =£ + - V Q (K t - \ e '^x H \ e iv n x _ e w n x H \ e -w n x\ K \ 

n=l 



n " t 



, Mt 

+ - J2 Q (K t - \e-' lv " x Hh iVnX + ffr» x H*e- iVnX - 2m \K t - ) Q . 



n=l 



The first sum on the right is m t , and I denote the second sum on the right by a t . By the analysis 
in the proof of Part (2) of Prop. 15. 1\ the terms in the sum of at are positive. Also, I note that when 
K,~ = k and \dL+ I = Id, then 

Q (k\e-' lvX Hh ivX + e ivX Hh- ivX -2m\k) Q = 2E[£ tn - £ t - \ T t - , \dL t J], 

since ±v occur with equal probability. 

The process St is the conditional projection of £t on to the set of processes whose value at time 
t depends only on T t - and the Levy process L t . Moreover, the process at is the projection of St 
that depends only on the jump times and the absolute value of the jumps \dLt\. Finally, At is the 
predictable projection of at- It follows that £t — St, St — at = mt, and at — At are uncorrelated 
martingales with the following inequality for their predictable quadratic variations: 

at > (M, M) t = {£ - S,£ - S) t + (m, m) t + {a- A, a- A) t > {m, m) t . 



Part (4): The predictable quadratic variation (m,m)t has the form 

(m,m) t = \ f dr f dvj(v) Q (K r \ e - ivX Hh ivX - e ivX Hh' ivX \K r ) \ (5.16) 

The above gives an expression through which I can examine the dependence of ^(m,m)t on Kt. If 

the expression in the integrand (|5.16|) had \p)q, q(p\, and Hz is replaced respectively by \p), (p\, and 
\P\, then I would have the explicit computation 

(k\e- ivX \P\e ivX - e ivX \P\e- ivX \k) = (k + v\ \P\ \ k + v > - (k - v\ \P\ \k - v) 

= \k + v\-\k-v\=2vS(k-v), (5.17) 

where the last inequality is restricted to \v\ < \k\. In my analysis, I will first work to bound the error 
of substituting Q {p\, \p) Q with (p|, \p), and secondly, I bound the error of substituting Hz with \P\. 

By the proof of Part (2) of Prop. [5.1[ the difference e~ mX H2 e wX — e H% e~~ wX has operator norm 
bounded by 2\v \ . The difference shares the same fiber decomposition of the Hamiltonian. Consequently 



1 1 



for 4> + ,(j)- € [—5,5) with <p + — v = 0_ + v mod 1, then the linear map e wX Hj e lvX — e wX H£ e wX 
on L 2 ([— 7r,7r)) has operator norm < 2\v\. 
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By (I4,3p . there is c > such that the distance between Bloch the vectors ipkji'k £ L 2 [[—tt, 7t)) for 
|fc| 3> 1 is bounded by 



+ rl(jfe) + l-ri(Jfe) < 



l + |/3(fc)| 



where the second inequality holds for some c' > 0, since r^(fc) and 1 — ri(fc) are bounded by a multiple 



of 



Hence, for i^>± 6 [— | , |) with ^± = A: ± v mod 1, then 



fee 



-ivX I 



<3 



;fc e' 



me 



ivX JvX 



me 



-ivX 



k) 



%\e~ ivX Hl + e mX - e wX H\_e-' wX 



%-Me-~ ivX Hl + e ivX 



e wX H\ e 



ivX 



[^ k \e- ivX Hl e lvX - e ivX H}„ e~ ivX 



+ 



<- 



Ad I v 



,,\e ivX Hl + e wX - e ivX Hj e~ ivX 



Wk 
\^k ~ 4>k 



(5.18) 



I now bound the difference Hz — \P\ when evaluated by kets \k') with \k'\ 3> 1. By the formula 
for the square root of an operator in terms of its resolvents, 



m 



H 



P 2 



00 dw 

^\w + H w + P 2 



dww 2 



1 



1 



w + P 2 w + H 



However, the difference between the resolvent of a Laplacian and the resolvent of the Laplacian 
perturbed by a ^-potential has a simple form |lj. To use this, I will focus on a single fiber from the 
decomposition ()B.l|h For w 6 M+ and (j) 6 [—\, \), the Green's function G^ w : [— vr,7r) — >• C for the 
operator to — A^ is given by the form 



) = -£ 



3 f(n+0)a; 



27r Z_. w + (n + ^2 



The operator on L 2 ([— 7r,7r)j determined by the integral kernel G^ w {x — y) is equal to (w — A^) . 
The difference between the resolvents in the d>-Sbei is 



1 



1 



+ P 2 w + H 



a 



l + aGU„(0) 



where the expression on the left in square brackets denotes the operator on L 2 ([— 7r, 7r)j corresponding 
to the ^-fiber, and the operator A^^ w has integral kernel A ( p tW (x,y) = G ( p^ w {x)G t j ) ^ w {—y). 
For k = 6 mod 1, then 



'k'\m -\p\\k')\ = JL 

47T 3 

a 



< 



4-7T 3 



dww^ 



dww? 



a 



1 



l + aG^(0) Vu, + (fc') 2 

1 Q 



(w + (A; / ) 2 ) 2 4|fc'|7r 2 ' 



(5.19) 



33 



Going back to (|5.16j) . I have the relations 

dvj(v) Q (k 

If i i 2 / | 

- / dvj{v) (k + v\ | k + v) — (k — v \ \k — v) + 0( — 

4 J\ v \<l\k\ ^1 + \h 



d 1 

dt {m)t ~ 4 



\ e - ivX Hh ivX - e lvX me~' wX \ 



v\<\\k\ 



\P(k)\ 



-f dvj(v)(\k + v\-\k-v\] +o( r 



U|<||fc| 



dvj(v)(2v) 2 + O 



1 



1 + 



a+ 



1 + 



The first equality follows from (|5.18p . and the second follows by f)5. 19|) for k' = k + v and because 
= (9( iqjg7fcji ) by the restriction |«| < The third inequality is by the definition of cr and 

Chebyshev's inequality through Ji v |>i| fc i cfe < 4|/c| _ \, where ? is the fourth moment of ^r-. 



□ 



Lemma 5.3. Let K t - = k for \k\ S> 1 and \dL t 



v < 



There exists an a > such that 



\^[£{K t )\F t -,dL t ] < - 



+ \m\ 



Proof. Note that 



E[E(K t )-E(K t -)\T t -,\dL t \] = l - Q {k\e-' mP He ivP + e ivP He-' lvP -2H\k) Q <v 2 , 

where the inequality follows from the same computation as in Part (1) of Prop. 15,11 However, the 
quantity Var \£{K t ) | T t - , dL t ~\ appears in the follow equation: 



E[E(K t ) - E(K t -)\T t -,\dL t \] =1^ Var [S(K t )\T t -, dL t = ±v] 



± 



+ -[E[£(K t )\T t -,dL t =v] -E[£{K t )\F t -,dL t 
+ E[£(K t )\F t -,\dL t \ = \v\] 2 -£ 2 {k). 
The bottom line is positive. The second line of (|5.20p can be written as 



(5.20) 



E[£(K t )\T t -,dL t = v]-M[£(K t )\T t -,dL t 



I e~ ivX m e ivX - e ivX m e - ivX 



and the right side is v 2 + Q( 1+ |^ fc ^ ) by the argument in the proof for Part (4) of Prop. 15.21 Since 
all three terms on the right side of (|5.20p are positive, it follows that Var[£ {K t )\Ft-,dL t = v] is 
0{ l+\p(k)\ ) as c l arme d- 



□ 



6 The limiting classical process 

I now shift my focus entirely to the classical stochastic processes K r whose probability density evolves 
according to the equation (|2.5p . In Sect. 16.11 I prove Thm. 13.11 which stated that t~^£ st , s G [0, 1] 
converges in distribution as t — > oo to the absolute value of a Brownian motion. Section 16.21 contains 
various lemmas related to the random variables Tn ,~t}~J n being approximately exponentially distributed 

with expectation v , where r n ,r n+ i are successive reflection times. Finally, Sect. 16.31 contains proofs 
of the lemmas in Sect. [3] and completes the proof of Thm. 12.31 
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6.1 A submartingale central limit theorem 

The difference between the quantities \k\ and £(k) = E^(k) is bounded by a constant, and the 
difference is even for \k\ 3> 1 as a consequence of the Kronig-Penney relation (12. 4h . Working 

with the process £ r is advantageous, since it is a submartingale by Part (2) of Prop. I5.2[ and the 
increasing part of the Doob-Meyer decomposition for the square increases linearly. Thus, the 

processes t~*\K s t\ and t~z£ s t, s E [0, 1] are close, and I will work with the latter. As before, M r and 
A r will denote the martingale and increasing parts the Doob-Meyer decomposition for £ r . The main 
result of this section is proof of Thm. 13.11 One of the key inputs for the proof is Lem. 16.14 which yields 
that £ r typically spends most of the interval r G [0,i] with values £ r S> 1. This is important, since 
some estimates, such as in Part (4) of Prop. [531 only become effective when £ r is large. Lemma 16.21 
offers a lower bound for the predictable quadratic variation of the martingale m% from Prop. [531 This 
lower bound shows that (m, m) r essentially grows linearly with rate a, and when combined with Part 
(2) and (3) of Prop. [531 this implies that t~zM s t and t~^m s i are nearly equal. 

The result of the lemma below holds for a general class of positive processes whose squares are 
submartingales, whose variances increase linearly, and that satisfy some Lindberg condition (to avoid 
the situation where the jumps are very large but very infrequent). The lemma is only stated for £ r 
here. The lemma implies that £ r typically spends most of a time interval r £ [0, t] with values on 
the order of t2~ p , < p < r; for large t, where "most" is all except for a total duration on the order 
0{t p ). Lemma 3.4 in [12] is similar, although the proof here employs different techniques. In the more 
general situation in which the second momentum of £ r grows at varying rates, similar statements can 
be made for a stochastic time change of £ r . 

Lemma 6.1. Let Qi,Q2,Q3 > 0, Q\ + Q2 + £3 = \, and q\ be strictly positive. For e > 0, define 
T e>t = Jo 1 ds X (t- sl £ s t < e). For all e,5 > , 

limsup t e2 ¥\T et > t'^S] < 64cj~^|. 
Proof. I begin with an application of Chebyshev's inequality to obtain 

t e *F[T €>t >t~ e3 8] < t —^E[T €!t \. 
Set £0 = Wi = 0, and define the stopping times q n , zu n < t such that 

w n = min{r G (?n-i, 00) | £ r < et ei }, q n = min{r G (w n , 00) | £ r > 2et Bl }. 
Let St be the number zu ra 's less than t. I clearly have 



T e ,t <t 1 ^ S n - W Tl 



n=l 



Moreover, notice that 



E[T M ] <E[t- x V?n-ron] ^t^E^supE^-ronln^Ht]. (6.1) 



n=l " "™eN 



I hence have an upper bound in terms of the expectation for the number of upcrossings "St and the 
expectation for the duration of a single upcrossing ? n — w n conditioned on the event n < 3*. By the 
submartingale uncrossing inequality |14| . I have the first inequality below 



r 1 E £ t ] E \£?] 2 (E \£n] +ttr)2 , , . 
E \St )< < < — < 2a2t2-^ e - 1 . (6.2 
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The last inequality is for t large enough. 

Next, I focus on the expectation of the incursions q n — w n . Whether or not the event n < Ej 
occurred will be known at time w n , so 

sup E[^ n -w n \n< E t ] < sup E [q n - w n \ F^ n ] . 

n€N neN, uieT mn 

In examining the expression on the right side above, I will set w n = and q = q n — w n . Since at is the 
increasing part of the Doob-Meyer decomposition for £ 2 — £q, the optional sampling theorem gives 

E w [?AT]= a^Eu [£* AT ~ £o] ■ 

If the process St were constrained to make jumps bounded by J > 0, then I could immediately reason 
that 

EcM = lim ct^EJ^V - £ 2 ] < a - 1 {2et^ + J) 2 , 

since £ s is either still in the interval [0, 2et ei ] at time s = q A T or has jumped out by a value of at 
most J. Plugging ()6.2p and the bound <7 _1 (2et ei + J) 2 for sup neN E[? n — w n | n < E 4 ] into (|6.ip leads 
to the desired result. 

For the case in which there is not a cap on the jumps, it is necessary to be more careful. Consider 
a modified process in which jumps £ s —£ s - greater than 2et Ql are removed. The modified process can 
be written as 

Mr 

£' r = ^ n -£ t ,)x(£t n -£ t - n <2et^), 

n=l 

for Poisson times t n < r. The modified process £' r is still positive, and for large enough t, (£' r ) 2 will 
still be a submartingale such that the increasing part of its Doob-Meyer decomposition grows nearly 
at the linear rate a (although growing at a lesser rate, such clS 2" 5 IS sufficient). 

Let q' be analogously defined to q as the hitting time that £' r jumps out of [0, 2et 2ei ]. By the 
definitions of these processes, I always have q < q', 

E w [q] = lim E w [q A T] < lim E w \q' A Tl 

< lim 2a- 1 E w [(^, AT ) 2 -(^) 2 ] < 32a" W", 

since £'^ A j- < 4ei ei . Combining this inequality with (|6.2p in (|6.ip finishes the proof. 

□ 

I will make a few comments about the Markov chain 8 n = Q(Kt n ), which is the sequence of 
momenta Kt n contracted to the torus T = [— j, -A at the Poisson times t n . Since the lattice components 
of the momentum jumps live on Z, they do not influence the jump rates for the contracted process 
9 n . Define the map (•) : L X (R) — > L X (T) by 

(q)(e) = ^(0 + i), ee[-\,\). (6.3) 

If qt n _ 1 , Qt n are the distributions for the momentum at times t n _i,t n , then (qt„) = T(qt„_ 1 ), where 
the operator T has integral kernel 

t(0 2 , e 1 ) = n~ 1 Y J j{o 2 -9 1 + 1 -), q x m € [ - ~ i). 
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As a consequence of (2) of List 12.21 ||(j)||oo = sup_i <6K i 'Yli i ^j{6 + |) is finite, and thus T maps L 1 
functions to L°°: 

\\ThU < n-'WhhJ^j^ + he L 1 (T) < Ill/ill!. (6.4) 

By (16. 4p . the process X r forgets its previous locations on the torus T exponentially fast. In 
particular, the process ®(K r ) will not not be disproportionately concentrated around the value zero. 
This is important for the proof of the following lemma, which relies on Part (4) of Prop. 15.21 since the 
bound in Part (4) of Prop. [O] is only useful for \0(K r )\ = \\Q(K r )\\n.{K r )\ > 1, i.e., when \@(K r )\ is 
not too close to zero. 

Lemma 6.2. The predictable quadratic variation for m r satisfies 

glog(l + t) 



E 

for some constant g > and all t G 



sup |<7S — t 1 (m,m) s t\ 
*e[o,i] 



< 



i 



Proof. By Parts (2) and (3) of Prop. 15.21 I have the inequality a > -^{m,m) r . Thus, the expression 
su Psg[o,i] I " 5 — t~ 1 ( m , m }st\ is equal to a — t~ l (m,m) t . For the mean-7£ -1 exponential waiting times 
e n = t n+ \ — t n between Poisson times, I can write 

Aft Aft 
t- 1 (m, m) t = t- 1 V{K tn )e n - t' 1 ( £ e n - t)v(K Mt ), (6.5) 



n=0 n=l 



where the second term on the right subtracts over-counting in the first term and will be of order 
0(t~ l ), since V is bounded. 

Approximating (m, m)t by the process st := X^=i ^[V(At„) | -^i„-i] e n, yields a difference bounded 

by 

t- 2 E[({m,m) t - s t ) 2 ] < oTLT x + 0(t~ 2 ), 

since sup^, V{k) < a and the e n 's are independent of each other and everything else. The conditional 
expectation E[V(i^t n ) | -7 7 t„_ 1 ] is a function Q : R — > [0, a] of K tn _ x . My expression of choice will be 
L drQ(K r ), which is nearly equal to sj except for two extra boundary terms of order 0{t~ l ). 

The key input for the proof will be an upper bound for a — Q(k) when \k\ > 2is . The probability 
that a Poisson jump v has absolute value \v\ > fs will decay superpolynomially for large t by (1) of 
List 12. 2| and thus 

a-Q(k)= [ , x dvY 3 ^\K v {k,i)\ 2 \o--V{k + i + v)\ + 0{r l s). (6.6) 
■/[-**,*«] ^ 11 

. . 2 l 

By Part (1) of Lem. 14.11 the sum of for i ^ I(k,v) will be < ct~3 for some c > 0. Since 

\v | < ts ; the values \k-\-i + v | must be > ie for i G I(k,v). Applying Part (4) of Prop. 15.21 the integral 
in (|6,6p is bounded by 



, 1 

M<*<i^ 2; M l + IL^J^I 



2a^ log(l + t6) i, 
< 1 1- V\t b , 
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where I have used that the values k + i + v have absolute value > ta and the shifts i G Z do not 
change the values k + v modulo h. The supremum in the second inequality is smaller than \i by (2) 



of List 12.21 Thus, for \k\ > 2t 6, the above gives 

ct-Q(A;) = ^^t ^ + 0(t~e . (6.7) 

Now, I am ready to bound the expectation of f^dr (a — Q{K T )^j. Define T t = t~ l f* drx(\K r \ < 
2te). By considering the events T t > te and T t < t&, I get the bound 



E 



r 1 / dr{a - Q(K r )) <aF[T t > fa] + aE[T tX (T t < t*)] + sup (a - Q(fc)) 



\k\>2t~15 



<a¥[T t > t - s] + <rt~« + log(i). 

For the event Tt < te, the third term on the right corresponds to the part of the trajectory in 
which \K r \ is > 2t6 ; and the second inequality uses (|6.7p . By Lem. 16.11 with S = 1, e = 2, and 
f?i = ^2 = £3 = g, the probability that T t > t~e has order t~6 (since |if r | £ r + 0(|K r | -1 )). 

□ 

The proof that the processes t~2£ s t, s £ [0, 1] converge as t — > oo to the absolute value of a 
Brownian motion relies on the following features of 

1. The process £ r is a positive submartingale with martingale component M r . 

2. The increasing part of the Doob-Meyer decomposition for is £q + or. 

3. When <? r 3> 1, then a — ^(M,M) t goes to zero (at least on average over long periods of time). 

4. The quadratic variation of M r satisfies a Lindberg condition. 

Theorem 13.11 follows from a more general theorem assuming conditions of the type above. In principle, 
if (2) does not strictly hold but the expectation of £t tends to grow with t, then one can obtain 
this assumption by considering a stochastic time change r t : £ Tt = £[. Assumption (3) is required to 
guarantee that the process t~2£ st behaves diffusively away from zero. It forbids, for instance, that 
the process £ r has the deterministic trajectory £ r = a^r?. To make use of (3), I will need to verify 
that £t spends most of the time at "high" values, and this verification is made through Lem. 16.11 
Assumption (4) is to guarantee asymptotic negligibility as is required in the martingale central limit 
theorem. Theorem 13.11 is similar to Thm. 4.1 in 1121 . 



Proof of Thm. E2J 

All "convergence in law" in this proof will refer to the uniform metric. Let M r and A r be the 
components of the Doob-Meyer decomposition for £ r , and define the submartingale G r = M r + 
sup < s < r — M s . The convergence in law of t~2£ st is implied by the following statements: 

(i) . The processes (at)~2G s t, s G [0, 1] converge in law as t — > oo to the absolute value of a standard 

Brownian motion. 

(ii) . The random variables sup < s <i | A st — sup < r < S £ — M r \ converge in probability to zero as 

t — > oo. 
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(i). If (<jt)~z M s t converges in law to a Brownian motion B s as t — > oo, it follows that t~^G s t converges 
in law to L(B S ), since the map F : L°°([0, 1]) defined by F(f)(s) = f(s) + sup r6 r 0)S i —f(r) is bounded 
with respect to the supremum norm. Moreover, by Levy's law, F(B S ) has the same distribution as 
the absolute value of a Brownian motion [62\ Ch.6]. 

Part (3) of Prop. 15.21 states that M r is a sum of the uncorrelated martingales m r and M r — m r . 
By Lem. I6.2| the predictable quadratic variation for t~^m s t, < s < 1 converges in distribution to 
as. By [36, Thm.VIII.2.13], it will follow that t~^m s t converges in law to a Brownian motion with 
diffusion constant a if the following Lindberg condition holds as t — > oo: 



t _1 E 



sup | m r — m r 

0<r<t 



0. 



The quadratic variation of m r is dominated by the quadratic variation of the Levy process by Part 
(3) of Prop. 15.21 and hence I have the first inequality below. 



t _1 E 



sup \m r — m r 

0<r<t 



<t 



<t 



sup I V Tl 

0<n<Nt 



n=l 



i 

=i"2 



dv j(v) v 



i 

4\ 2 



The second inequality is Jensen's, and the equality holds because the jumps v n occur with Levy rate 
j(v). By the assumptions on the jump rates j(v), the fourth moment is finite, and the bottom line 
tends to zero for large t. That proves the convergence in law for t~2m s t- 

Next, I argue that t 2 supQ< s <^ \M S f — rn s t\ converges weakly to zero as t — )■ oo. It will then follow 

that t~?M s t converges to a Brownian motion. Since the martingales m r and M r — m r are uncorrelated, 
the formula holds: 



a > 4-{M,M), 
dr 



d 
dr 



(to, to) 



d 

r + d~r {M 



m,M — to), 



(6.8) 



With Doob's inequality, the inequality (16.81) . and Lem. [6T2| I have 



£ _1 E[ sup \M r 
0<r<t 



TO, 



I 2 ] <4t- 1 E[|Af r . 



=4t _1 E[(M — to, M 
<4E [o--t _1 (TO,TO) t ] 



TO) 



0. 



(ii). Note that G r is the smallest positive submartingale that has Doob- Meyer decomposition with 
M r as the martingale part, so £ r > G r . I begin with a technical point extending the discussion of (i). 
Let A ( s ] be the increasing part in the Doob-Meyer decomposition for t 1 G^ t . I should expect that As 
converges to as, since {at)~^G s t converges to the absolute value of a standard Brownian motion, and 
the increasing part of the Doob-Meyer decomposition for the square of a standard Brownian motion 
B s is s. My first task will be to show there is convergence in distribution 



sup (as 

0<s<l 



4 t} ) 







(6.9) 



for large t. 



The value A s has the form A, 



(0 



when t 2 G s t is likely to jump to zero: Z s := t 1 7Z drF(£ r ,G r ) for 



t- 1 (M,M) st 
(t) 



At) 



where is a term related to the times 



F(G r , K r 



E 



I Mr - M r - + G r \ X ( M r ~ M r - < -Gr) G r , K r , M r - M r - / 
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The dependence on K r enters through the statistics for M r — M r - conditioned on the event that 
r G R + is a Poisson time (i.e. N r - M r - / 0). It follows that t~ 1 (M,M) st - Ay is increasing m s 
and so 

sup (as -A®) = (a-r l (M,M) t ) + zf . 

0<s<\ 

By Part (1), the first term on the right converges in distribution to zero as t — > oo. 

Now to address Z± \ I always have the inequality F(G r ,K r ) < and for the regime \G r \ 3> 1, 
then 



F{G r ,K r ) = E 
< E 
= O 



G r , K r , M r -M r - / 



\M r -M r - +G r \x(Mi 

\M r - M r - \ 2 x{\M r - M r - | > G r ) G r , K r , Mr ~ jV r - ^ 

1 



l + \P(K r )\ 



To see the order equality, I first break M r — M r - into the two parts (M — m) r — (M 



and 



m r - and use the simple inequality 



\Mr-M r - +Gr\ x(\ M r 



M r - | > G r ) 



< 4 ra, 



X[\m r 



> ^G r ) + 4|(M - m) P - (M - m) r -\ 2 



to separate the terms in the expectation. By the reasoning in the proof of Part (3) of Prop. 15.2 



\m 



2 is smaller than the square of the Levy jump \L r — L r -\ 2 . Since the Levy jumps have 
exponential tails, the expectation of \m r — m r - \ 2 x{\ m r — m r -\ > \G r ) must decay exponentially in G r 
(and therefore be negligible). The expectation of | (M — m) r — (M — m) r - | is smaller than ppw^ — jt 
by the proof of Part (4) of Prop. 15.21 Note that £ r > G r and £ r ~ | -KV | to conclude that when G r is 
large, then \K r \ is also large. Doubling the bound to cover the smaller term \m r 



above, I can 



write 



E 



t- l ll / drF(G r ,K r ) 



< E 






Jo 



dr- 



8Ka 



1 + \(3(K r ] 

Define T' et = t~ l J'q dr x(t~^G s t < e). For some c > 0, the following inequality holds: 

limsupP[r e ' t < 5] < c- 



(6.10) 



t— >oo 
1 



This can be either proved using that t~zG s t converges to the absolute value of a Brownian motion 
and a little estimate involving Gaussians or by a small modification of the proof of Lem. 16.11 with £ r 
replaced by G r - Hence, for any 5 > I can pick an e small enough to make the probabilities P[T e ' t < 5} 
uniformly small for large t. By a similar argument as in the proof of Prop. HT2l the right side of (|6,10D 
converges to zero, which finishes the proof of (|6.9p . 

With (j6.9j) as a tool, I go to the core of the proof. Since t" 1 ^ and t^G 2 ^ are both submartingales, 
taking Doob-Meyer decompositions allows me to write their difference as a martingale plus the 
difference of the predictable increasing parts as and A { s ] . I can also write the difference between 
t~ 1 £ 2 t and t~ 1 G 2 t in an expression involving A st — sup < r < st —M r as follows 



2t- l £ st {A, 



St 



sup -M r ) + t-\A st 

0<r<st 



sup 

0<r<st 



-Mr 



1 C 2 

st 



t 1 G 2 st 



as 



(6.11) 
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Since the left side is made up of positive terms, I can immediately conclude that 



sup t 1 (A st - sup -M r ) 2 < sup + sup (as-A^p). 

0<s<l 0<r<st 0<s<l 0<s<l 

Hence, showing that the terms on the right converge in distribution to zero for large t will be sufficient. 
I have already shown sup 0<s<1 (crs — A^) converges to zero, and now I turn to the martingale 
component. 

Since the left side of (|6.1ip is positive, it follows that —AdP < as — AP for all s G [0, 1], and thus 

E[\M®\ X (M®<0)] <E[as-A®]. (6.12) 

Let r be the stopping time when AdP first jumps > e or r = 1 when values above e are not reached 
by time s = 1. By the optional sampling theorem, I have the equality E[.M^] = 0. Since M { r ] has 
mean zero, the equality below holds: 

sup M® > el < e^ElMPxiM^ > e)l < e^RllM® \x(M® > 0)] 

0<s<l - 1 

= e^EllMPlxiM® < 0)] < e^EllM^lxiM^ < 0)] 
Ke^Ela-A?]. (6.13) 

The first inequality is Chebyshev's and the last is (|6.12p . Since the random variables a — A± are 
bounded by a and converge to zero in distribution as t — > oo, the expectation on the right tends to 
zero. 

□ 

6.2 The waiting-times between momentum reflections 

This section focuses on bounding and characterizing the length of the time intervals between "Bragg 
reflections" for the momentum process. The stochastic process K r exhibits a singular rate of sign- 
changing in comparison to an ordinary random walk due to the lattice components of the momentum 
jumps. The lattice jumps yielding reflections are "macroscopic", since they transport the momentum 
to a value that would require many steps through the Levy component to reach. This unusual sign- 
flipping characteristic was not apparent for the study of £ r = E2(K r ) in Sect. 16-H since £ r depends 
only on the absolute value of the momentum. The results below are directed toward showing that 
the waiting-time r m+ i — r m between successive reflections r m , r m +i is approximately an exponential 
distribution with mean u~ 1 \K Tm \ and is independent of the subsequent fluctuations in momentum. 
This idealized behavior emerges for high beginning momentum |i^ Tm | S> 1. 

It will be useful to work with an idealization of the momentum process that makes the same 
Levy jumps with rate j(y) except that the additional lattice jump m G Z takes only the values 
or — n(k + v) with the probabilities defined below. Define the functions H-(k) = -§^i{k 2 + j&£z)~ \ 
R_(A;) = n_ (^0(/c)n(fc)) , and R+ = 1 — R_. If the current state of the momentum is k and the 
Levy increment is v, then the conditional probabilities for the lattice component are defined to be 



m = < 



R + (k + v), 

(6.14) 

-n(k + v) K-(k + v). 



This process describes a particle whose momentum is reflected at the next Poisson time based on 
a weighted coin flip whose bias is determined by the sum of the current position k and next Levy 
increment v. Of course, the simplified statistics is a useful approximation only for time intervals when 
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the momentum is high. Note that the idealized process closely resembles the process proposed in 
Sect. 12.2.61 I will refer to probabilities and expectations in the law above by P and E, respectively. 
For a technical reason explained below, I define the particle to make an additional jump —ii(Kq) with 
probability T-(Kq) at an infinitesimal time after zero . 

Consider a particle beginning with a momentum k with \k\ 3> 1 and making jumps according to 
the simplified law above, and let r be the first reflection time, i.e., the non-trivial lattice jump. I will 
sketch why the random variable 4jj is approximately an exponential with mean v~ , It is equivalent to 

consider a discrete-time random walk X n , n £ N beginning from k with jump increment density 
and a hitting time N marking the first tails outcome for coin tosses with tails probabilities R_(X n ). 
As long as X n does not stray too far from k, the probabilities R_(X n ) are approximately 

R_(x n )«n_(|fc|e(* n )). (6.15) 

The problem is thus further reduced to the contracted Markov chain 9 n := Q(X n ) living on the torus 
T = [—4,3) with the transition operator T : L X (T) (defined above Lem. 16 . 2[) . The probabilities 
H_(|/c|# n ) decay rapidly for \9 n \ ^> so 9 n will typically require many steps to score a reflection. 

Since the chain 9 n is exponentially ergodic to the uniform distribution over T, the probability of a 
reflection on a given time step is approximately 2 J* T d9Tl-(\k\9) « m. The random variable m should 
thus be close in law to a mean-a -1 exponential for \k\ ^> 1. This conclusion is consistent with the 
ansatz (j6. 15|) . since the fluctuations in the random walk X n over [0,iV] will be on the order 0(|fe| a ), 
which is small compared to \k\. Since the Poisson times have rate 1Z and v := IZa, the distribution 
for will be approximately a mean-z/ -1 exponential. 

The connection between the momentum process and the simplified law (|6.14p is not obvious, since 
the approximations for the probabilities \K v (k,n)\ 2 in Prop. [3~T1 reduce for \k\ S> 1 to expressions 
r ei (k)r e2 (k + v), €\,C2 £ {±} when n G I(k,v) and zero for n ^ I(k,v). For a hint of the link, note 
the identity 

R_(fc) = 2r_(fc)r+(A;). (6.16) 

An intermediary law between the original momentum process and (|6.14p is given by the process that 
accompanies a Levy jump v from the momentum k with four possible lattice jumps with values and 
probabilities given by: 

— n(/c) r_(A;)r + (fc + v), 

m - \ + r+(k)r_(k + v), ^ ^ 

n(k) — n(k + v ) r_(/c)r_(fc + v), 
r + (k)r + (k + v). 

Naturally, when n(k) = n(k + v), then the probabilities are summed. I refer to these statistics by 
P',E'. There is an underlying equivalence between the laws (|6.14p and (|6.17p such that they may be 
embedded on a single probability space in which the trajectories match for specific realizations except 
for time periods around reflection times where the values become staggered. To demonstrate this, I 
will consider only skeletal chains, since the waiting-times between jumps are independent of the states. 
As before, let v m , m G N be a sequence of independent random variables with density ^r-. Let the 
(non stationary) Markov chain X n have increments whose law depends on the parity of n > 1 as 



X n — X n _\ 



r + (X n _i) n odd, 

-n(X n _i) r_(A n _i) n odd, 

vt± r + (A„_! + vn ) n even, 

-n(J n _i + U|) r_ (X n _i + vii) n even, 
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where the second column on the right lists the probabilities for the lattice component. Thus, the 
jumps v m occur at even time steps and coins are flipped at each time step to determine whether there 
is an extra lattice jump. The chains X2 m +i and ^2m for m > are each Markovian and have the 
statistics of (|6,14[) and (|6.17|) . respectively. This result depends on the symmetries j(v) = j{—v) and 
r_(fc) = r_(fc — n(/c)). The two idealized laws are interchangeable for the purposes of this section. 

Let S : M — > {±1} be the sign function. I will use the term "sign-flip" in a technical sense to refer 
to a Poisson time t n such that for some odd m G N 

S(K tM ) = S(K tn _J = -S(K tn ) = -S(K tn+1 ) (6.18) 

and S(Kt n _ r ) = —S{Kt n _ r+1 ) for r G [l,m). In other words, there is an odd-numbered sequence of 
sign-changes ending at the time t n . Only sign- flips in which there is a single sign-change are likely to 
occur in practice. The complication in the definition results from the fact that there is a comparatively 
high probability of a sign-change at the Poisson time after a sign-change, and I would like to avoid 
counting these "fake" sign-flips. For example, suppose the beginning momentum is k with \k\ S> 1 and 
the next two Levy jumps are v\ and v%. By the approximate statistics (|6.17p . the probability assigned 
for the lattice jump — n(fc + v\) at the first Poisson time is r + (/c)r_(A; + v±), and the resulting landing 
point is k + x>\ — n(k + vi) ~ —k — v\. The value r_(fc + v\) can not be too small if this jump is 
likely to occur. The probability that the particle makes another lattice jump n(/c + v±) at the second 
Poisson time to land at k + v\ — v 2 is approximately 

r_(fc + vi- n(k + vi))r + (k + v\ + v 2 - n(k + vi)) > ^ r -(& + 

where the inequality is by the lower bound r + > i and the symmetry r_(/c — n(k)) = r_(fc). Thus, 
the second sign-change has a nonnegligible probability > ^r_(A; + vi) regardless of the value v>}. 

A sign-flip time is not a hitting time, since it requires the information from the following Poisson 
verifying that the sign of the momentum does not change again at next the momentum jump. However, 
this will not give much information. Consider Kq = k with \k\ S> 1, and let q',q G L 1 (R) be the 
distribution for the first Poisson time with and without conditioning to not change sign, respectively. 
Using the estimates (1) and (3) from Lem. ETTI it can be proven that the probability of a sign-change 
can be at most ~ ^. Hence, I will have the bound 

q'(p)<2q(p), peR. (6.19) 

In particular, there are no density peaks near lattice points ^Z, where approximations tend to be 
weaker. 

The following proposition is the main result of Sect. 16.21 and the proof is in Sect. 16.2.21 The 
inequalities in Prop. [6731 have the following purposes: Parts (1) and (2) are characterizations for how 
close the random variables Tm \p~~] m are to mean-i^ -1 exponentials for |i^ Tm | 3> 1. Part (3) bounds 
the amount of time that the momentum process spends performing "fake sign-flips" before the actual 
sign-flip occurs in the technical sense (|6.18|) . Part (4) states that the probability of a macroscopic 
fluctuation in the absolute value of the momentum before a sign-flip is small. Part (5) shows that 
the fluctuation in the absolute value of the momentum over the time interval between sign-flips is 
nearly uncorrelated with the length on the time interval. Notice that in Prop. 16.31 the momentum 
process is conditioned not make a sign-flip at the first Poisson time. This is exactly the situation 
that I have for the process K r , r G [r m ,oo) when given the information T- and r m is a sign-flip. 
This does not explicitly cover the case that r m is not a sign-flip as when r m = $j for some j G N or 
\K Tm \ ^ [jjl-Krm-il) ||-Kr TO _i|]> These events do not occur frequently enough to generate nonnegligible 
contributions to the quantities that I study (in fact, Part (4) of Prop. 16.31 bounds the probability of 
the event \K Tm \ £ [^K^^ \, §|l£r m _i |])- Without conditioning the momentum to remain the same 
sign after the first jump, the estimates in Parts (1) and (2) may be distorted if Kq is near a lattice 
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point jZ, since r will have a comparatively high probability of occurring at the first Poisson time. 
However, Prop. 16.31 still offers bounds for the moments of m in that case. There are other natural 
conditions on the beginning momentum that will yield the results of Prop. 16.31 The assumptions that 
K r begins with the value Kq = k and is conditioned to avoid a sign change at the first Poisson time 
can be replaced by the assumptions that K r begins in a distribution given by a density q G 
concentrated around the value k and that the contracted density (q) E L 1 (T) is bounded. 

Proposition 6.3. Let £ > 0, Kq = k for \k\ S> 1, and K r be conditioned not to make a sign change 
at the first Poisson time (i.e. S(Kq) = S(K tl )). Define r to be the first time that either K r has a 
sign-flip or \K r \ jumps out of the set [\\k\, depending on what occurs first. For fixed £, m > 0, 
there exist C and 70 such that the following inequalities hold for all k and < 7 < 70.' 



1. 

2. 



|„-m 



< 



c 



V 



< 



C 



3. E 

4. P 

5. E 



J T dr X {S(K r ) + S(K )) 
\K T \ i [\\k\, §|*|] 
(\K T \ - \K Q \)r\ < Clfcf+t. 



< ClogCfcj) 

- pep • 



Recall by the remarks preceding Lem. 16.21 that the chain 6 n = @(Kt n ) on the torus T 



"4> 4, 



is Markovian with bounded transition operator sup^ g 2e j T(02, 9\) < ^. Since the lattice jumps do 
not appear in T, the contracted process for the idealized momentum process is exactly the same. The 
torus chain has detailed balance, since the kernel is symmetric T(#2,#i) = T(#i,#2) as a consequence 
of the symmetry j(v) = j(—v). The process is thus time-reversible, and its stationary state is the 
uniform distribution. Due to the boundedness of the kernel, the process is exponentially ergodic, and 
even in the supremum norm I have the existence of an e > such that 



sup I (T n h) {9) 



2 < e' 



sup \h(9)-2\. 



(6.20) 



The laws P and P' determine two probability measures on the sequence of momenta Kt . . I denote 
the total variation distance between the measures on sequences of length M by || • ||var,M- Let Am be 
the event \K tj \ G [||*|, ||*|] for < j < M, and x(^4m)P be the positive measure that agrees with P 
for events in Am and assigns events in A C M zero weight. The key inputs for the proof of Lem. 16.41 are 
the estimates (1) and (3) of Lem. 14.11 



Lemma 6.4. Let Kq = k. There is a c> such that for all k, 



\x{A M )F - X (Am) 



log(|fc|) 



Proof. Let refer to the process that obeys the law (|6.17p for the first n jumps and the original 
law for the remaining jumps. The variation norm for the difference x(Am)F — x(Am)^' is bounded 
by 



\\x(a m )f - X (AmW'\ 



M 



Var,M 



<Ell^f)F^-x(^M)P (n) 

n -v - X (A M r {n} 



Var,M 



n=l 



<Msup||x(^Af 



IVar.AP 
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where I have used the triangle inequality with a telescoping sum determined by inserting the measures 
X(A)¥^ n \ Let dL n and d£ n denote the Levy and lattice jumps at the nth Poisson time. The laws 
p( n ) and p( n_1 ) differ only at the rath Poisson time in the probabilities that they assign for the lattice 
jump d£ n = m G Z. I can bound a single term from the sum above through the inequality 



| X (A M )P (n - 1) -x(^M)P (rl) || ValM < / dp\q n (p) ~ Qu(p)\ 

JR 



K 



■P 



(6.21) 



'[d£ n = m | dL n = v, K Tn _ x = p] -P' [di n = m \ dL n = v, K Tn _ 1 

where q n and q' n are the densities 

q n ( J) )=E^[5(K tn -p)x{A n )] and q' n (p) =E'[5(K tn - p) X (A n )}. 
The first inequality (|6,2ip follows since Am C A n for M > n. I will show that the right side of (|6,21D 



is O(iSgp) for \k\ » 1. 

By (1) of List 12,21 and Chebyshev's inequality, I can replace the integration J R dv in (|6.2ip by 
the restriction \v\ < ^\k\ with an error that decays exponentially for \k\ ^> 1 (since the integrand is 
bounded by one). Under P', only the lattice jumps rra G I(p,v) have nonzero probability. For the law 
P, Part (1) of Lem.O yields that for \p\ G [±\k\, §|fc|] and \v\ < \\k\, 



F[dl = m\dL n = v,K Tn _ 1 



P\ 



E 



K v (p,m)\ 2 = 0(\k\~ 



when \k\ 3> 1. Since the density q' n (p) is supported on the interval \p\ G [^\k\, §|fc|], I can restrict the 
summation ^ m in ()6.2ip to rra G I(p,v) with an error 0(|A;|" 2 ). 

Next, I treat the difference (|6.2ip for m G I(p,v). The different cases are handled similarly, and 
I will take m = —n(p) ^ —n(p + v). I make the restriction — n(p) ^ —n(p + v), because I will 
apply Part (3) of Lem. 14.11 to obtain the errors for the approximations of \k v (p, — n(p)) | , and the 
cases n(p) = n(p + v) and n(p) ^ n(p + v) correspond to different approximations in the lemma. By 
definition, 

P' [d£ n = -n(p) | dL n = v, K Tn _ x = p] = r_ (p)r+ (p + v), 



so I am interested in the difference 
dpq'n-iip) 



M<m av n 



(p, -n(p)) | 2 - r_ (p) r + (p + v) 



(6.22) 



I will change summation variables through p = mi + 9\ and p + v = rai2 + #2 for ?rai,rai2 G and 
6*i , ^2 G [—4,2)- The domain \v \ < ^\k\ is roughly the same as | m\— m,2\ < \ I fc I ■ so the expression (|6 . 22 p 
is approximately 



E 



TV 



L 414; 



d9id9 2 q n -i{mi + 6*i)j("i2 - mi + 6 2 — 9\ 



0<|mi-m 2 |<i|fe| 



1112-m 



1+02 _ 01 (rrai + #1, -2mi) I - r_(mi + 9 1 )Y + (m 2 + 9 2 ) 



3.23) 
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By Part (3) of Lem. 14.11 there is a c > such that the sum of terms m\ ^ m 2 in (j6.23|) is smaller than 



0<|mi-m 2 |<2;|fc| 



<^rll(d>lloo||<j>IU / , d0v 



< 32< W Ml II /.-Ml l0 S( 1 + S) 



^ ||\"fn — 1/ || OO || \J I ||UU Ij^lJ 

where the map (•) : — > L X (T) defined in (|6,3p contracts densities on M to densities on the torus 
T = [—4, 4). In the first inequality above, I have bounded \m\ + 6*i | — 1 (1 + |mi#i|) -1 from below by 
+ It an d commuted the sums with the integrals. By the discussion preceding Lem. 16.21 



{q'n)\\oo < [KT^Iloo = llT^Hoo < ^IK^IU, (6.24) 



and where T is the transition kernel for the law (|6.17p . The norm ||(j)||oo is finite by (3) of List 12.21 
Note that q' n {p) < (7~q' n _i)(p), since (Tq' n _i)(p) includes a contribution for \p\ outside the interval 
[g|A;|, 2 \k\] ■ 

□ 

Lemma 6.5. Let tj for j E N be the Poisson times and Kq = k. For \k\ large, there are c, C > such 
that the following inequalities hold: 

1. F[S(K ) = S(K tl ) = -S(K t2 ) = -S(K t3 )] > fa 

2. F[S(K ) = S(K tl ) = -S(K t2 )] < g. 

The same statements hold with the statistics P replaced by P. 
Proof. 

Part (1): The laws P and P' determine measures on the sequences Kt y Let || • ||var,4 an d the set 
be defined as in Lem 16.41 The probability of the event A\ is order 0(\k\~ 2 ) for both P and P', since 
the random variables |-ftTj +1 | — \Kt-\ have uniformly bounded second moments and by Chebyshev's 
inequality. Thus, 

ii^-^iiv ar ,4+°(i fc i") = u^-x{M)n^=o(^-), 

where the second equality above is by Lem 16.41 However, the variational distance ||P — P'||var, 4 will 
bound the difference between the probabilities 

>[S(Ko) = S(K tl ) = -S(K t2 ) = -S(K t3 )] - F'[S(K ) = S(K h ) = -S(K t2 ) = -S(K t3 )] 

Therefore, I can substitute P with P' with an error o ( | | 1 ) . 

One way in which the event S(Kq) = S(K tl ) = —S{K t . 2 ) = —S{K ta ) can occur is if the Levy jumps 
satisfy | | < |fe|2,i = 1, 2, 3 and the corresponding lattice jumps are respectively 0, — rv(k+V\ +v 2 ), 0. 
This possibility gives a lower bound 

F'[S(K ) = S(K tl ) = -S(K t2 ) = -S(K ts )] 
M,M,M<|fc|2 Tl Tl Tl 



> / i —5 7B ^5~ r +( k ) r +( k + v\)r-{k + v x + v 2 ) 

■/l«ll,l« 2 l,l«3l<|fc|2 K ]< - K 

x r + ( k + vi + v 2 - n(k + v\ + v 2 )) r+ (k + v% + v 2 - n(k + v% + v 2 ) + v 3 ) . 
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I will use the identities r+(k') = r+(k' - n(k')) and n_ (2- 1 @(k')n(k')) = 2r + (k')r_(k') for k' = 
k + v\ + v 2 . Since r + is > ^ and j{v2) > — for \v 2 \ < 1 by (3) of List 12.21 I have a bound from below 
given by 

T- I / dv 2 J -^-n_{2- 1 e(k + v 1 +v 2 )n(k + v 1 + v 2 )) 

^ J\vi\<\k\? rC J [-1,1] l< ' 

2/z72. 7[« x |<|jfc|3 ft J[-i,i] 
a 



iin\ky 

where I have used that L d/cn_(A:) = a. 

Part (2): By the argument in Part (1), I can approximate P with P'. There are four possible pairs 
of lattice jumps that will give a sign change for the second jump and not the first: for the first and 
— n(k + v\) or — n(fe + v\ + v 2 ) for the second, or n(k) — n(k + vi) for the first and — 2n(k) + n(k + v\) 
or n(k) — n(k + v\) — n(k + v± + ^2) for the second. The cases require the same analysis (with more 
or less messy notation), so I will analyze the case in which the jumps are and — n(k + v±): 



f 

•/|vi|,|V2|<|*|^ ft ft 



j(vi)j(v 2 ) 

r + (k)r + {k + i>i)r_(fc + vi)r + (k + v\ + v 2 ) 



<\( jMu^(2- 1 e(k + v 1 )n(k + «i)) 



']«i|<|*| 3 
|rra— fc|<-v/jfc| 



7m 



< 



In the first inequality, I have used that r + < 1. The norm ||{j)||oo is smellier than \x by assumption 
(2) of List E2J 

□ 



6.2.1 Results for the idealized law 



Lemma [6.61 concerns the simplified momentum process (|6.14p . or alternatively (|6.17p . and it is applied 
in the proof of Lem. [67 
of the lemma states that random variables 



to obtain the analogous result for the original momentum process. Part (1) 

1 "* are approximately mean-a -1 exponentials for 



J K Tn I 3> 1 . The proof is based on the heuristics in the paragraph ()6. 15|) . 

Lemma 6.6. Let > and Kq = k with \k\ 3> 1 be conditioned not to change signs at the first 
Poisson time. Let r be the first Poisson time that either the process makes a sign flip or its absolute 
value jumps out of the interval [\\k\, For fixed (, there is a C > such that for all fc£K, the 

following inequalities hold: 



1- su PaGK+ 



< C\k\ 



-i+C 



E 



(\K tM \-\K \)x(M T >M) <C\k\? 



-i+C 



where M := [\kf\. 
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Proof. 



Part (1): Since I am considering the probability of the event P[A/" T > a\k\\ that is determined 
when the first reflection occurs (or, less likely, when l-PC/vvl leaves the interval [^\k\, the problem 

can be formulated in terms of a random walk X n = k + v\ + . . . v n , where the increments v n are 
independent and have density ^=r^; and the reflections are decided by flipping coins with heads weight 
IT + (2 n(X n )0(X n n . The problem is thus translated to characterizing the distribution for the first 
tails outcome. This scheme ignores the special rule for the first momentum jump for the law (|6.14j) . 
but that will not generate a difference. I am also ignoring the technical definition (|6.18p for a sign-flip, 
which is only relevant when working with the laws P and P' for which the double flips are likely. Let 
F C N be the set of times n with either a tail outcome or such that \X n \ lands outside [\\k\, §|&|]. 
Also, let N be the smallest element in F. 

First, I will show that P[iV > a|fc|] decays exponentially for o> 1. By a simpler argument than 
in the proof for Part (1) of Lem. 16.51 there is a B > such that 

P[n i F | X n -i] < 1 - ^ (6.25) 

for all |X n _i| G 5 \k\- With an iterated conditional expectation and an application of ()6.25|) . then 

P[iV>n] =I[x(MF)---x(n^F)] 

= E[ X (1 i F) - - - x (n - 1 i F)P[n £ F | X n ^]] 

< (l-^jE[x(HF)... X (n-HF)], 

where B is the constant in the lower bound from Part (1) of Lem. 16.51 By induction, P[iV > n\ < 
(1 — yf|) n - It follows that pry has finite exponential moments E[e 7 ' fc ' lN ] for < 7 < B, which are 
uniformly bounded for all large |fc|. 

By the above P[iV > m] decays exponentially with rate at least (1 — jf|) m - Hence, the values 

P[iV > a\k\^ and e~ aa will be both small for a > \k\**, and I can focus my analysis on a < \k\^. Let 
me set n = [ a l^|J an d "fh = Ll^l J- can rewrite P[iV > n] as 

P[iV>n] = E[x(l£F)---x(ngF)] 

= E [x(l £ F) ■ • ■ x(n - 1 £ F)P[n £ F | F n (n - rh, n) = 0, X h _ rh } 

= F[N > h - 1] / dpr A}A ip)¥[rh £ F | A> > rh, X = p] , (6.26) 
Jr 

where rh )r h{p) is the density for X^-m conditioned on the event N > h — 1. The third equality above 
holds by the strong Markov property. By definition, the density r^ ^(p) will have its support over 
\p\ G [\\k\, \\k\]. I will argue that the integral in the last line of (16. 26ft has the 3> 1 asymptotics 

I dpr h ^(p)¥[m£Y\N>m,X Q =p]=l-^ + 0{\k\^- 2 ). (6.27) 
Jr \k\ 

The result (|6.27|) can be applied inductively in (|6.26|) for P[JV > h — 1], . . . , ¥[N > rh + 1] to obtain 

n rh ft— 1 * • j 

P[iV>n] = (l--^) n " T p[iV>m] + 0( 7 |A ; | 2 ^ 2 ) £ (l - ^ P[N > j] 

' ' j=m+l ' ' 

= e~ aa + 0(\k\*- x )e- Aa , 
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where A = B A a. The second equality follows by a simple argument yielding that P[iV > m] — l| = 
0(|fcp~ 1 ). Since ( > is arbitrary, the result would follow. 

I will separate the analysis showing (16. 27ft into proofs of the statements (i)-(iii) below. For \p— k\ < 
|fc|2 + '' and | A; | ^> 1, the following statements hold: 

(i). The difference between the left side of (|6.27p and 



! dprn,m{p) P[m <£ F | N > rh, X = p] 

\p-k\<\k\? +( 



(6.28) 



decays superpolynomially in 
(ii). P[m £ F | N > rh, X = p] = 



\p-k\<\k\? 

(i) . The value J 



i +< dpr^(p){l -*) = !-& + 0(\k\ 2 <- 2 ) 



i+ij dp Tfa^ (p) is smaller than the probability of the event sup 0<n< , A \X n — k\> 



|p—fc|>|fc|*" 

|fc|2+f. The probability P[sup 0<ri<ft X n — k > \k\5 +< *] will decay superpolynomially fast for \k\ 3> 1, 
since 



;I*I'P 



sup X n - k > |fc|2+£ 

L 0<n<n 



< E 



sup e 

L 0<n<n 



(x„-fe)|fcr^-2 



< 4E 



,(X ft -fc)|fc|-£-£ 



<4(l + 2a'|fc|- 1 -«) n 



< 4e 



2a' 



(6.29) 



where the first inequality is Chebyshev's, the second is Doob's maximal inequality, and a 1 is the second 
moment of ^=A. The last inequality 
bound P[sup < n < ft fc — X n > |&|a+C] 



moment of The last inequality uses the cap h < \k\ 1+( >. The same argument can be applied to 



(ii). I first show that, when beginning from Xq = p for \p\ £ 



the probability of the event 



m ^ F is nearly the same with or without conditioning on the information N > rh. I will bound the 
probability P[m ^ F | Xq = p] in (ii') below. First, I need to show that P[iV > m Xo = p] is nearly 
1. This requires another inductive argument concluding that 

F[N > rh | X = p] = E[x(N > rh - l)P[m $ F | | X = p] 



> E[x{N >rh-l)\X =p][l 
2a\ 



2a 



+ ot 



> i 



\c-i + (\k\t- 2 ). 



The first inequality follows from the analysis below. I can assume that X^ = X^_i + % has absolute 
value |X^j_i + Vm\ > \\k\ by (|6.29|> . so 



1 -P[m £F\X, 



rh—l 



dv ' 



X A _i+i.|>i|fc| 



n 



-n_ (2- 1 n(X ril _ 1 + v)e{x 1h _ l + v)) 



< 



1 

K 



.1 II 

4 >4J 



ft 



d0 (j)(^-e(x A _ 1 ))n_(2- 1 |A;|e) 

i d9U-(2- 1 \k\9) 



< 



2a 



n\k\ 
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Now, I use that P[N > rh | Xq = p] is close to 1. Let h, h! be the distributions for X^-i when 
beginning from Xq = p and conditioned or not conditioned on the event N > rh, respectively. By 
similar reasoning as above 

|P[m i F | N > rh, X = p] - P[m i F | X = p] \ 

dw(h(w)-h'(w)) [ dv J ^-n^(2- 1 n(w + v)Q(w + v)) 



< \\h-ti\\i( 
8a 2 



2a 



n\k\ 



+ o(\k\- 



< 



TZ 2 



™\k\<- 2 + O (\k\< 



(6.30) 



where the last inequality follows from 



\h - h'Wx <2¥[N >rh\X Q =p] < 



4a 



K 



\ k \<-i + (\k\(- 2 ). 



(ii'). Approximating the unconditioned expression P[ra ^ F | Xq = p] is more straightforward. I will 



denote the density P[X^ = w | Xq = p] = T m {6{- — p)) by qm,p- I have the following approximations 
P[m^F| Xq=p] w J &g™ iP (!«)n + (2" 1 n(w)e(w)) 

dw qrh, P {w) n+ (2 _1 n(u;) Q(w)) 



w—p\<rh 



w—p\<rh 



dpq rhjP (w)U + (2- 1 n(p)Q(w)) 



+ 0(\p\<- 1 [ dp g Aj ,(w)n_ (2- 1 n(p) 



By the same argument as (|6.29p . the probability that \K r \ leaves the interval \p — rh, p + rh] after rh 
steps after beginning from p will be superpolynomially small. Hence, the approximation in the second 
line above with the restriction of the integration to [p — rh, p + rh] will have a superpolynomially small 
error. The error in the last equality comes from the inequality 



407T 

n+(x) -n + ( y ) < —^Tl-{x 
1 1 cr 



m - \y\ 



when Ixl 



n 1, 

\v\ \ < ~\x\ 
\v\\ - 2 i 



(6.31) 



for x = 2 1 n(p) @(w) and y = 2 1 n(w) Q(w). Finally, n + (2 1 n(p) @(w)) only depends on w modulo 
^, so I have that 



w— p\<rh 



dp q m , P {w)Yl + (2 1 n(p)Q(w)) 



I a- a! 



^(4 p )(0)n + (2- 1 n(p)0) 



< 



d9(q lh:P )(0)n + (2' 1 n{p) 9) 

_i i\ ' 
414,1 



1 i) 

4 '4 / 



d9U + (2- 1 n(p) 9) 



1-T-r + Ol 



where q' mp {w) := qm,p( w )x(\ w — p\ < to). The second approximation above follows by (|6.20p . since 
(lrh,p) = T m ((<5(- — p)}) must be superpolynomially close for \k\ 3> 1 to the uniform distribution 2. 
The order equality on the bottom line uses that \p\ > 
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(iii). Plugging in the result (i) to (I6,28j) . and expanding \p\ 1 around k results in 



a 



J\p-k\<\k\? +c v \p 



1 -^) + T^I i ^y^,m(y + ^) + 0(|A;| 2 ^ 2 ). (6.32) 
I I I I ^l?/|<|fc|^ + 



It remains to show that f < fc i +c dyyr^^y + k) is 0(|/c| 2< »). This will require an analysis of the 
density r^m and related densities r ft m> which can be written 

E[U + (2- 1 n(X 1 )e(X 1 ))---U + (2- 1 n(X fl _ 1 )Q(X ii _ 1 ))5(X h _ lh -p)} 



>(P) 



E[n+(2-in(X 1 )0(Xi)) • • •n+(2-in(X ft _ 1 )e(X ft _i))] 
_ I [n+ (2- 1 n(fc)9(X 1 )) • • • n+ (2- 1 n(fc)6(X^_ 1 )) ( 5(^_ A - p)] 
E[n + (2- 1 n(fe)G(X 1 )) • • • n + (2-in(fc)9(^_ 1 ))] 

I [n+ (2- 1 n(fc)e(x A+1 )) • • • n+ (2-X*09(-Xft-i))*(-Xft-* - P )] 



E[n+(n(fc)e(x A+1 )) • • • n+(2-in(fc)e(x fl _ 1 ))] 

I claim that the following order equalities hold: 



sup 

\y\<\k\? 



2) and sup 

|y|<|fc|2+< 



(6.33) 
(6.34) 
(6.35) 

(6.36) 



I will return to a discussion of (|6.36p below in (iii'). It follows from (|6.36p that 



,dyy{rw(v + k)-rl A (y + k)) < 0{\k\<-*) / dyWl^y + k) = O {\k\*) . 

In that case, it is sufficient to bound the expression f 1,, dyyr'l ~.{y + k). The density ri' - can 
be written as 

*ft,m(p) 

_ / R du;g A ( U ;)E[n + (2- 1 n(A;)e(X A+1 )) • • • H+ (2-^*0© (*«-i)M(*n-m =p) | ^ 



it; 



E[n + (2-in(fc)9(X^ +1 )) • ■ -n + (2-in(A;)e(X fi _i))] 

where is density for the random variable X^. The density A corresponds to the distribution for 
Xn-ih conditioned on the event Fn[m+1, n— 1] = 0. Moreover, the integral of qm(w)x{\ w ~ k\ > l&F) 
will be superpolynomially small, since the jumps in the random walk Xi, i £ have mean zero 

and exponential tails. Hence, I can split Xn-m, — k in to a sum of components Xfa — k and X^-m — Xm, 
where \Xm — k\ is typically small and Xh-m — has density Um,n,w G L 1 (M) given by 



|(P): 



E[n + (2- 1 n(fc)6(X A+1 )) • • • n + (2- 1 n(fc)9(X A _ 1 )) ( 5(X A _ A - X A = p) | X* 



w 



E[n+(2-in(fc)G(X A+1 ))---n + (2-in(A;)e(X ft _ 1 ))] 

when Xfa = w and conditioned on the event Fn[m + l,n-l] =0. Since the arguments of the 
expectations defining Urh,h,w depend only on the contracted random walk Q(X m ), the density satisfies 
the shift invariance Urh,h,w(p) = ^mhw+^mip) f° r 171 e ^- This shift symmetry gives the inequality in 
the second line below: 



dyyr'l^y + k) =0{\k\<) + 



<0( 



+ 



1 dwqrh{w) / j dyyUrh.h,w(y) 
\w-k\<\k\? +< - J ' \y\<\k\*2 +c - 



dO(qrn)(0) / 1 dy y Urn^oiy) 

1 ^|j/|<|fc|2+ C 



(6.37) 
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The equality above follows by approximating Xh^rh—k using and restricting the integration 

over w 6 K to \w — k\ < \k\i +< *. 

The density {q^) £ L 1 (T) will be superpolynomially close to the uniform distribution over T = 
[—\, \) by (|6.20p . However, if (q^) is replaced by 2, then the resulting expression is zero, since the 
symmetry Urh,h,e{v) = Urh,h,-e(—y) holds as a consequence of the symmetry j(v) = j(—v) for the 
jump rates. Thus, the last term on the right side of ()6.37p is also 0(|A;|^). 



(hi'). I will focus on bounding the difference 



l| over the domain \y\ < |fc|2+£, since the 
order equality on the right side of (16. 36ft follows by simpler analysis. Let G be the event sup <j<ft \Xj — 

k\ < |/c|2 + ^. By (I6.29p . 1 — F[G] decreases superpolynomially with large \k\. The difference between 
the expressions in the denominators of (|6.33p and (|6.34p is smaller than 



E 



|n+ (2- 1 n(fc)e(x 1 )) • • • n+ (2- 1 n(k)e(x il - 1 )) 

- n + (2- 1 n(x 1 )e(x 1 )) • • •n + (2- 1 n(X fi _ 1 )9(X fi _ 1 )) 



G 



h-1 



<-J- E [ n + (2" 1 n(fc)e(X 1 )) • • • n + (2- 1 n(fc)6(X,_ 1 )) 

k |n_(2- 1 n(A : )e(x j ))n + (2- 1 n(A ; )e(x J+1 ))---n + (2- 1 n(x^ 1 )e(x^ 1 )) 



\Xj-k\ 



G 



The above inequality follows from inserting a telescoping sum and using the inequality (|6.3ip . 

Since I have conditioned on the event G, then 1 J fc I < |/c|^ _ 2. Removing the conditioning on G 
will make another superpolynomially small error, and a single term from the sum above is bounded 
by 

|A;|^5irn + (2- 1 n(A;)e(Xi)) •••n+(2- 1 n(A:)e(X j _i)) 

x n_(r 1 n(i)6(i J ))n + (2- 1 n(fe)e(i J+1 )) ■■•n + (r 1 nft_ 1 )e(i A _ 1 )) 

<2/i 4 |A ; |^-ll[n + (2- 1 n(fc)e(X 1 ))---n + (2- 1 n(^_ 1 )e(X A _ 1 )) 

The conditional density for Q(Xj) given 0(X,_i) is bounded by /i as a result of assumption (2) from 
List 12.21 (and the same for 0(X, + i) given Q(Xj)). I can make use of this to obtain bounds in which 
Q(Xj) is decoupled from the other variables, and where Q(Xj) and Q(Xj + i) are integrated with respect 
to the constant density (i over T. Integrating out Xj yields a factor 2 J* T d9H + {2~ 1 n{k)9) m < 

jj^r • Moreover, I can put the variable back in the expression by inserting 2n + (2~ 1 n(£;)0(X,)), since 
2n + > 1. The conditional density for @(Xj) given @(Xj_i) (and 0(X, + i) given Q(Xj)) is bounded 
from below by ^ by (2) of List [2J2J I can use this to recouple Xj with the other variables, so that 
the expression bounding the difference is a constant multiple of the denominator in the definition for 

r h,rh' 

The same argument as above applies in order to bound the difference between the numerators 
of ^M) and (EM}. 

Part (2): I will sketch the proof. Let ( < /3, M = [\kf\, and rh = [|A;| C J. As in Part (1), I can 
phrase the problem in terms of the random walk X n . Since there is no actual sign- flipping for X n 
(only coin tossing), I can remove the absolute values from Xm,Xq in the expression: 



E 



(X M - X ) X (N > M) 



E 



u + (2- 1 n(x 1 )e(x 1 )) •••n + (2- 1 n(x M )e(x M ))(x M - k) 
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where the approximation uses that there is a superpolynomially small probability that the event 
N < M occurs because \X n — k\ > v for some n < M. By the same reasoning as in Part (1), I 
can replace the factors n(X r ) by n(fc) in the arguments of the functions II + with an error of order 
0(|&p _1 ), and I can also remove the factors 11+ (2~ 1 n(Xj)Q(Xj)^ from the expectations for j < m. 
By the triangle inequality, I can bound the right side above by 



E 



U + (2- 1 n(k)Q(X 1%+1 )) ■ ■■U + (2- 1 n(k)e(X M ))(X M - X* 



rh—l 



+ |E[n + (2- 1 n(£ ; )e(X A+1 )) ■■■U + (2- 1 n(k)Q(X M ))(X n+1 -X n )]|. (6.38) 

n=0 

The first term will be superpolynomially close to zero for large \k\, since the distribution (q#i) for 
Q(X 1 f l ) will be superpolynomially close to the uniform distribution and by the symmetry argument 
used in the proof of Part (1). The probability of the event N £ [n + 2, M] when X n+ \ = p for n < N 
and \p — k\ < j\k\ will be Od/cj' 3-1 ), so the nth term on the lower line of (|6.38[) is equal to 

n+(2- 1 n(A:)G(X A+1 )) •••n + (2- 1 n(A;)e(X n+1 ))(X n+1 -X n )]| + OO^ 1 ). 

Finally, I have the equality 

iO) 



E 



n+(2- 1 n(A:)G(X n+1 ))(X n+1 - X n ) 



X r . 



dv v 



K 



■n^(2- 1 n(k)(X n + v)). 



I can use the exponential decay of j(v) to cap the jumps v by |fc|^ with a small error, so the expression 
above is approximately 



dv v 



v\<\k\i 



K 



U^(2- 1 n(k)(X n + v)) < 



K 



den_ (2- 1 n{k)e) = o(|fc| c_1 ), 



where I have followed the usual method for bounding expressions with integrals including j(v) and 



6.2.2 Proof of Proposition HT3] 



□ 



Part (1) of the lemma below makes another step toward showing that the duration between between 
successive sign-flip times r n ,r n +i is approximately an exponential distribution with mean v~ l \K Tn \. 

Lemma 6.7. Let t, be as in the statement of Lemma\6. 7\ There exists a C > such that all k 



1- sup aeI 



W 



> a 



< C|fc|- 



-i+C 



2. 



E 



\K tM \-\K \)x(M T >M) <C\k\^+<, 



where M := [\k\P\ 
Proof. 



Part (1): By Part (1) of Lemma 16.61 and the equivalence between P and P', it is sufficient to show 
that 



sup 

ael + 



> a 



A Mr 
[\k\ 



> a 



01 



ic- 1 



(6.39) 
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for all C > 0. 

First, I will show that the probability ¥[M T > n] decays exponentially for large n. Let F C N be 
the set of times n such that t n is a sign-flip or \Kt n \ jumps out of [^\k\, §|fc|]. I have the following 
relations 



M T > 3n 



3n 

E[n^(^ F ) 



3=1 

n-1 



n 

<e[JJ x (3j-i £f) 



J'=l 



E 



[(n^(^- 1 ^ F )) IP [ 3n - 1 ^ F | K *3n- 



< 1 



n-1 

E[n*(3j-1£F) 



where the second inequality uses the Markov property and Part (1) of Lemma [63] to bound P[3n— 1 ^ 
F | Kt 3n _ 3 ] . My purpose of removing the terms j = 0, 1 mod 3 in the product above was to accommo- 
date the technical definition for a sign-flip, which requires information from the next Poisson time. 
This argument can be applied inductively to obtain a bound for ¥[M T > n] that decays exponentially 



as (1 



< 7 < B and k. The same arguments hold for the law P'. Since the probabilities P[^ > 
and P'f^ > a] decay exponentially in a G M + with a rate that is uniform in k, I can take a cut-off 

\a\ < \k\z as in the proof of Part (1) of Lemma 16.61 

Let the set Am and the variational norm || • ||var,M be defined as in Lem 16.41 Notice that the 
event t^J > a is contained in the set Am for M = [|^| 1+ ^J. The variational norm for the difference 

x(^4m)P — x(^4m)P' bounds the difference in probabilities for any event contained in Am, and thus 
the first inequality below holds: 



for B 



ML 

3 



I remark that the expectations E[e 7 "^] must be uniformly bounded for 

a] 



sup 

a<]fcp 





[Mr I 




[M T I 




p 


iw >a 


-P' 


iw >a 





< \\x(A M )F - x(A M ) 
log(\k\) 



IVar.M 



< cM- 



0( 



The second inequality is by Lem 16.41 

Part (2): If |-K"t M | — \Kq\ were a bounded random variable, then this result would follow immediately 
from the proof of Lem. 16.41 By Part (2) of Lemma 16.61 it is sufficient to prove that 



E 



K tM \-\K \)x{M T >M) 



E' 



K tM \-\K \)x{M T >M) 



O 



log(|fc|) 



(6.40) 



As in the proof of Lem 16.41 I will break down the difference into a sum of differences between the 
expectations E^ -1 ), E( n ) of \K t \ - \K \ conditioned on N T > M: 



M 



n=l 



J2\^ n ~ l) (\K tM \-\K \)x{M T >M) -EW (\K tM \-\K \)x{M T >M 



(6.41) 



The difference K tn — K tn _ 1 is equal to v n + £ n , where v n , l n are the nth Levy and lattice jumps 
respectively. For m G Z, let T(m; p, v) be the closest element to m in I(p, v) (where ties can be assigned 
arbitrarily). It is convenient to split £ n into two parts F(£ n ; K tn _ 1 ,v n ) and £ n — V (£ n \ K tn _ 1 , v n ) that 
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are treated separately. Let K' t n be the momentum at time tu if the nth lattice jump is replaced by 
F(e n ;K tn _ 1 ,v n ): 

K[ M ,n ■= K tM + K tn ^ - K tn + T(£ n ; K tn , v n ) + v n . 
I have the following bound for the nth term of the sum (|6.41j) when Kt M is replaced by K[ n in the 
expectation E^™ -1 ); 



jg(n-l) 



K' tMjn \-\K \) X (N T >M) 



E (n) 



K \ X (N r >M) 



< / dpq[ n _ x (p) / dv 



n 



x \v + i| 

ie/(fc,u) 



|«« (p, m) | 2 - P' [d£ = i | dL n = v, K tn _ x = p] 



m£Z with 
i=T(m; kn—ijV) 



where the density q' n is defined as in the proof of Lem 16.41 Recall that the q' n {p) is supported in the 
set [|H,§M]. By Part (1) of Lem. [4.11 the sum over \k v (v, m)| 2 for m ^ I(p, v) is 0(|A;| -2 ). Thus, 
it is sufficient bound 



\v\<\k\? 



dv J M £ \ v + 



K 



i£l(k,v) 



\k v ( P , [d£ = i\ dL n = v, K tn _, = p] 



.42) 



where the restriction of the integration in v to the set \v\ < \k\ 2 will have a superpolynomially small 
error by the decay of j(v) . I can apply the same estimates as in the proof of Lem l6.4l to show that (|6.42p 
is 0(!2f>). 

To complete the proof, I must bound the error of having substituted |-?Q M | by \K[ M n \ in the 
expectation EO- 1 ) above: 

E^T \K tM \ - \K' tM \ X {K >M)]< [ dpq'^ip) [ dv J -^Y] \m - T(m; p, v)\ \ Kv (p, m)\ 2 . 



Again employing a cutoff \v\ < \k\z and using that the support of q' n (p) is over \p\ G [\\k\, I can 

apply Part (2) of Lem. 14.11 since |m — T(m; p,v)\ = dist(m, I(k,v)). This yields a c > such that 
for all k, 



E («) 



\K, 



\ K L 



< 



\k\- 



Therefore, a single term from the sum (|6.4ip is 0( log Q^ ). Since there are M = [\k\@\ terms, the 
result follows. 

□ 



Proof of Prop, \6.3l 

Part (1): Similar to Part (2). 



Part (2): The Poisson times t n can be thought of as the sum of independent mean-7?. -1 exponentials 
e m , m G N: t n = Ylm=i em - ^ ^ s sufficient to prove the statement of (2) with r replaced by M T and v 
replaced by ^ = a. This can be seen with the short calculation: 



E 



E 



71 \M 



,7? - -2- 

^ |fc| 



E 



E 



+ o(ifcr i ). 
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By the remark in the proof of Part (1) of Lem. 16.71 the expectation of e is uniformly bounded 
for all k E M and 7 in < 70 < 70 for small enough 70 > 0. Using the following expectation formula 
for a random variable X: 



E[. 



,-yXi 



poo 

I d ai e ai ^[X> 
Jo 



I can bound the difference in the case 7 < 70: 



E 


r 

e 7 


a 






a — 7 



< 



7 _1 MW) 
< 2C 7 



2 C -i + 2 



p 




_ e -aa 


+ f 




m- a - 




Jb 



da 'j e 



«7 



\k\ 



> a 



da'j e 



a(7~7o) 



0( 



|2C-1> 



Aft 1 log(|fc|) 

The second inequality uses Part (1) of Lem. 16.71 for the first term. For the second term of the second 
line, the probability Efy^p > a] is smaller than e _a7 ° by the proof of Part (1) of Lem. 16.71 

Part (3): The quantity J Q T drx(S(K r ) 7^ S(i£o)) sums the amount of time that K r spends up to 
time r with the opposite sign of what it began with. By the definition of r, K r cannot spend two 
consecutive Poisson times t n -%,t n < r with S(K tn _ 1 ) = S{K tn ) 7^ S(Kq), and thus S(K tn ) / S(Kq) 
implies S(K tn7l ) ^S(K t J. 

I can rewrite Jq dr\{S(K r ) ^ S(Kq)) in terms of the exponential waiting times e n as follows: 



E 



[ T dr X {S(K r ) ^ S(K )) 
Jo 



.A/V-l 



=E 



[ £ X {S(K tn ) ^ S(K ))e n 

n=2 

j 

<^E[ Y, ^[S(K tn )^S(K tn ^)\K tn _ 
n=2 

E[AT T ] 



A/V— 1 

n=2 

A/V— 1 



< 



■711 



cE 



< 



2c 



where I do not count n = 1, since the sign of the momentum has been conditioned not to change at 
the first Poisson time. The probabilities in the summand on the third line are smaller than a constant 
multiple c > of by the same argument as in the Part (2) of Lem. 16.51 To apply the argument 
from Lem. RT5| I must use that \K r \ > ^\k\ for r < r. The last inequality is for \k\ large enough, since 
E[m] approaches u^ 1 at k — > 00 by Part (1). 

Part (4): In the notation from Part (1), ej^ T+ \ is the waiting time from r to the next Poisson time. 
The random variable r is not a hitting time, since determining r requires information from the time 
t' := T + e/vv+i) but r' is a hitting time. The results for Parts (1) and (2) will also hold for r', since it 
includes only an extra contribution from an additional independent exponential random variable with 
mean Hr 1 . 

I will apply a submartingale argument using the fact that S s ~ \K S \ + O ( | I ~ 1 ) when \K S \ S> 1. 
As before, let M r ,A r be the Doob- Meyer components of the submartingale E r — £q. Recall that M r is 
a sum of two uncorrelated submartingales m r and M r — m r by Part (3) of Prop. 15.21 By an argument 
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similar to Lem. 16.21 I can show (m, m) T is close to err and thus (M — m,M — m) T is small, since 
^{m, m) r + 4p{M — m,M — m) r < a. To see that <tt — (m, m) T is small, notice that 



E[ar - {m,m) T ] = E^ x {K > n) {a - V{K tn ))e n 



n=0 



r r 1 o 

< E[^2x(K > n - l)e n E[(a - V(K tn )) X {\K tn \ G -|fc|]) 



n=0 



where the inequality follows because the event M T > n is contained by the intersection of Af T > n — 1 
and |.KjJ G The exponential random variables e n are independent of everything else, 

so they can be left outside the conditional expectation. By the same argument as in the proof of 
Lem.E21 

log(|fc|) 



E 



{<r-V(K u )) X (\K u \e[±\k\,l\ 



tn-l 



4a 
< — 

- n 



\k\ 



Using the above and the result of Part (1), 

4a 

E[(M -m,M -m) T ] < — 



,!^EH<^ll«IUiog(l*l). 



Now applying standard martingale arguments, I can bound M r — m r and A r 



E 



sup | M r — m r + A r | 

0<r<T 



|2" 


< 2E 


1 |2" 

sup \M r — m r \ 






Lo<r<r' -1 



+ 2E[A^] 



2a 2 



< 8E[(M T , -m T ,) 2 ] + — E[r 2 ] 

= 8crE [( M -m,M -m) T ,]+ — ^ E [r 2 ] 

|k| 

^^ll^llooiogd^D + s^V. 

The second inequality is Doob's for the first term, and for the second term I use that A T < £r, which 
is a consequence of the identity a = 4p(M, M) r +26 r -^A r along with the fact that £ r ~ \K r \ > for 
r G [0, r). By Chebyshev's inequality, it follows that the probability sup < r < T \M r — m T + A r | 2 > 
wiUbeO(^). 

The final part of the story is m T . The m r martingale is well-behaved, since the absolute values 
of its jumps are less than those of the Levy process by Part (3) of Prop. 15.21 and the Levy jumps 
have exponential tails. On the other hand, the moments of r are finite by Part (1). The probability 
\m T \ > tz +s , 5 > will decay superpolynomially, since r has order for \k\ S> 1. 



Part (5): This follows by an extension of the analysis in Part (2) of Lemma 16.71 which I will not 
include. 

□ 



6.3 Proof of Lemmas EOlEOl 

For r,t G R+, s G [0, 1], and n G N, let the processes M r , A r , Y s , mi , and N r ; the filtration J>; and 
the times r n be defined as in Sect. [3j In the proofs of this section, I will treat all the times r m as if 
they occur through sign-flips rather than elaborating on the exceptional occurrences in which T m — 
for some j G N or \K Tm \ £ \h\K Tm _A } ||-K" rm -il]- Although it is slightly incorrect to neglect those 
cases, the omission avoids some messy and unenlightening case considerations, and the estimates in 
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the proofs below imply that those cases have negligible contribution. Whereas the number of sign-flips 
will be on the order 0{t^ ), the expected number of ?j over the time interval [0, t] has the bound a^ts +L 
for < i -C 1 by (|6.2[) . and the event that \K Tm \ £ [| |ifr m _i|j il-^V m _i|] f° r some r m < t is unlikely 
to occur (see (|6.63|) ), 



Proof of Lem. \3.4\ The proof of the Lindberg condition for M r follows from the analysis contained in 
(i) of the proof of Thm. ETTl Also, the random variables t~zM s t for t E R+ and s E [0, 1] are uniformly 
integrable, since 

sup E[(rW si ) 2 ] = t^E[(M t ,M t )} < a, 
se[o,i],teK + 



where the second inequality is by Part (2) of Prop. [5? 
To see the Lindberg condition for mi , notice that 



sup I 

0<s<l 



m 



(t)| 



--t * S up Ar n \K Tn \ 

0<n<N t 



<t * ( sup t 2 \K r 

• 0<r<t 



sup 



l<n<N t \K Tn 



(6.43) 



The random variables t » ( sup < r <t converge to zero as t — > oo, since t i\K s t\ converges 

to a Brownian motion with respect to the uniform metric by Thm. 13.11 Moreover, by Part (2) of 

Prop. 16.31 there is a C > such that for small enough 7 > 0, the expectations E[e |JfTnl | J- T -~\ are 

smaller than C for all n and i > 1. With Chebyshev's inequality, Pl" ,'^ 7 "' 1 > t§ I T -1 < Ce 
Applying this with an inductive argument using conditional expectations, then 



-Ttf 



sup — — - < ts 

L l<n<N t \K Tn \ 



> E 



1 - Ce" 7 '* 



1\N, 



> E 



1 - Ce" 7 " 



The second inequality follows because the number of sign-flips will be less than the number of Poisson 
times, and the number of Poisson times is typically oc t. 

The random variables for s E [0, 1], t E M+ are uniformly integrable, since the variances are 
uniformly bounded: 



N, 



sup E 

se[o,i],teR+ 



m 



--E 



t-2 ^|^ T J 2 |Ar n | 



<^E 
v 



<-E 
v 



n=l 
N t -1 



t~2 \K Tn \ 3 Ar n 



n=l 

_3 

t 2 sup |A r 

0<r<i 



+ -E 
v 



t 2 sup \K r 



0<r<t 



+ 0(t~5) = 0(1). 
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The second inequality uses that ^n^T 1 ^ r «- — ^or ^ ne ^ rs ^ inequality, I use nested conditional 
expectations twice to replace the factor |Ar n | 2 by ^AT n \K Tn \ in the sum. In doing so, I invoke Part 
(1) of Prop. 16.31 to get the approximations 



E\\AT n \ 2 \T 



and 



^ T jE[|Ar n || 7 



For each approximation, I multiply the upper bound by a factor of 2 to cover the error of the approx- 
imation. Uniform bounds of E[sup < r <( t _ ~|-fC r | m ] , m = 3,4 for large t are obtained by using that 
\K r \ < E r and applying Doob's maximal inequality to the submartingale E r . Bounds on the fourth 
moments of £ r are contained in Part (1) of Prop. 15.21 

□ 



5* 



Lemma 6.8. Fix some < i < in the definition for the times r n , n G N. ^4s t 

convergence in probability 



oo, there is 



sup 

0<s<l 



N st -1 



-Ar„ 



= 1 J T n 



Proof. Since my estimates will depend on \K r \ S> 1, I first show that the contribution to Yg that 

3 

accumulates when \K r \ <C \k\ & will be negligible. Let 5 > 0, then 



sup i * 

0<s<l 











< F 











dr\K r \ X ( \K r \ < 2t 

< 128 a' ^S'H" 21 

where the last inequality is for t large enough by Lem. 16.11 with ^ = : 



ds X (t-& +L \K st \ < 2 > 5t-s +t 



(6.45) 
i — t, and 



e = 2. Hence, for any 5, I can pick t large enough to make (|6.45p arbitrarily small. It follows that 

Ns * r(T n +Ar n )Ast 



sup 

0<s<l 



n=l Jr " 



drK r 







as i — > oo, since the contribution during incursions is negligible. The upper bound (r n + Ar ra ) A st for 
the integrals can be replaced by r n + Ar n , since 



sup t 4 
0<n<N t 



t„+At„ 1 , j 2 

c2r|i£" r | < i - 2 f sup t~2|X r 

Tn ^ 0<r<< 



sup 

Kn<N ( 



The right side goes to zero for large t by the discussion following (16. 43ft . For the same reason, I can 
replace the upper bound of the sum by N st — 1. 

□ 

For the remainder of the section, I will set i = in the definition for the excursion intervals 
[<j, t _i,ro n ), and thus indirectly for the definition of the times r n . The rate of decay computed for 
various expressions in the proofs would be a little slower if I kept i > 0, but i can be chosen arbitrarily 
small anyway. The following proof will rely heavily on applications of Lem. 16.71 and Prop. 



Proof of Lem. \ 3.S\ . The proof of Part (1) follows from similar analysis as in (i) of the proof for Lem. [34] 



For Part (2), I can approximate Y} by the expression t~4 Y2n=i~ It" T " drK r for t S> 1 by Lem 
The result follows by showing the convergences in probability (i)-(iii) below. The differences in (i) 
and (ii) involve coarse-graining approximations in which the random time intervals Ar m are parsed 
into shorter intervals A mj „ with duration on the order 0(|-KV m | ) f° r some /3 chosen from the interval 
(0, |). I define the following notations: 



Lm 



\K \P 

A/" Tm +Ar m! 



-Mr 



rn.n 



A, 



oj(m) 

Wrm +nL m 1 



Moreover, I define T m ^ n as the Poisson time preceding r mjn when n > 1 and r mj o = r m . Also, I will 

abuse notation by identifying Jr m0 with T - . 

I will show the following convergences to zero in probability: 
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(i). sup 0<s<1 



5 



v^N st -l / r 
Z^m=l ^ J 7 



+Ar m 



drK r - S(K Tm ) J2to l A m ,n|^r m ,„ 



0, 



(ii). supo^ rf^i 1 (5(K Tm )En=o _1 A m ,„|Kr m ,J-^ m Ar„ 



(iii). supo^^i t-*Y^iKr m ^r m \T 



0. 



The expression i 4 Xmi=l K Tm (/S.T m — E[Ar m | i 7 ,.-] ) is obtained by the right term in (ii) minus the 

expression in (iii) and differs from the expression for by the substitution in the upper summand 
of N S £ with N S £ — 1. By the Lindberg condition in Lemma |3.4| the difference is negligible for t > 1. 

(i). Over an interval r £ [r m , r m + Ar m ), the process K r has the same sign except for isolated Poisson 
times at which it jumps to the opposite sign and back again at the next Poisson time. My first step 
will be to bound the net effect of these rogue sign changes. The first inequality below follows because 
\K r \ e [±\K Tn \, l\K Tn \]: 



E 



sup t 4 ^2 



Nst-1 » Tm +At„ 



0<s<1 



m=l 



N t -1 „ 

<EH V -\K T IE 

m=l 



dr X {K r ^ K Tm )\K T 



m=l 
N t -1 



< 3Ci/t"E 



< 3Cz/t~* 



[E E I 



m=l 

0. 



dr X {K r ^K Tm ) 



.46) 



The second inequality follows by Part (3) of Prop. 16.31 For the third inequality, I have used that 
E[Ax m | J 7 ^] ~ ^ _1 |K Tm | for large t by Part (1) of Prop. 16.31 and I doubled the bound to cover the 
error. The last inequality holds since the Ar m 's sum up to less than t. Hence, I can take the sign of 
K r to be constant over the time intervals [T m ,T m -i- A7 m ). 

It is useful to approximate the expressions in (i) by replacing K r with S(K r )£ r , since £ r is a 
submartingale with convenient analytic properties. By (|6.46j) and because \K r \ = £ r + 0(\K r \ ), I 
have that 



sup 

0<s<l 



N s t-l 

*- f E 

m=l 



T m + AT„ 



drK r -S(K Tm ) 



dr£ 7 



0. 



I can also replace \K r \ with £ r in the expression t 4 Sm=l 1 ^(^r m )Y^n=o 1 ^m,n\Kr m>n \, and the 
difference (i) reduces to 



N st -1 f r m +Ar m N st -1 L m -1 

t-4 ^ S(K Tm ) / dr£r-t-i Yl S ( R rJ E A ™^r m , 

1 J Tm 1 . n 



m=l 



5 

< t"4 



N st -l£ m -l 

E E 

m=l n=0 



m=l 



n=0 



dr(f r -fr m ,J +t"?( sup £ r ) ^ A ™,L m , (6.47) 



0<r<t 



N t -1 



m=l 



where the second term bounds the under-counting by Ar m — (T m ^ m — r m ) < A mj £ m for the length of 
the interval [r m , r m + Ar m ). 
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Next, I will show that the second term on the second line of (16. 47ft goes to zero. The intervals 
between successive Poisson times have mean TZ~ l and are independent of everything else, which gives 
the equality below: 



N t -1 

e[ a ^ 

m=l 



1 r Ni_1 
-E[ u(m) 



< 



t 



-fd- 



n 



I Nt-l 



< 



m=l m=l 
N t -1 



K 



< 



P 



m=l 



K 



(6.48) 



The first inequality follows because cj (m) < |AV m |^, and \K Tm \ > is. By Part (1) of Prop. E3 I have 
that E[Ar m | ET Tm ] > ±\K Tm | for large enough t. The last inequality follows by removing the nested 
conditional expectations and Y^rn=i Ar m < t. Recall that (3 G (0, |), so I have | + |/3 < |. For a 
<J € (0, i - |/3), then (16^81 implies that H+* Em^ 1 A m converges to zero in probability for large 
i, and Thm. [3Tl implies that t~?~ s sup 0<r<< £ r converges to zero in probability. Hence, the product 
goes to zero. 

To bound the first term on the second line of (|6.47p . I use the Doob- Meyer decomposition £ s = 
M s + A s and the triangle inequality to get 



i"4E 



N st -iL m -i ,r m , n+ i 
sup I E £ / dr(£ r -£ r7i 

°- s - 1 m=l n=0 ^ r m.« 

N t -lL m -l 



J 



^1 j. ^ m j. /"l m,n + l 

EE E [ dr(M r -M Tm J 
m=\ n=0 ^r m , n 



+ i _ 2E 



N t -li m -l 

[E E 



m,n + l 



m=l ra=0 



(6.49) 



where I have inserted nested conditional expectations for the martingale term. For the A r term, recall 
that £% is a submartingale with increasing part £q + at by Part (1) of Prop. 15.21 It follows that 
£A r < ^(cj - Jr(M, M) r ) < at~l for the "high energy" part of the trajectory £ r « |if r | > if. The 



last line in (|6,49p is therefore less than at s . 

For a single pair m,n from the sum on the second line of ()6.49p . the following inequalities hold: 



E 



dr(M r -M r „ 



< E 



E 



E 



dr{M r -M Tm>n ) 

Tm : n-\- 1 

dM r (r TOl n+i - r 
r m ,„+i d(M, M) 



Ft 



dr- 



dr 



( r m,n+l - r) Tv m>n 



The first inequality above is Jensen's inequality. The last inequality holds since the predictable 
quadratic variation {M, M) r grows at a rate < a. 
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The second line of (|6.49p is bounded by 2(^)2 multiplied by 



t"4lE 



N t -lL m -l 

[E E 

m=l n=0 



3 

A 2 



<ft"2 ( SUP 

s>i szT(s) 



tv 1 3\ N 4 -l i m -l 



t 4E| > > o,'2(/);) 

m=l n=0 



7?. 2 SU p 



r(s + l)\ _5 



N t -1 



He[ e w ^( m ) E Am ^ 



m=l 



n=0 



<^"2 SUp 



S >i s^T(s) 



t~iM\ sup £ r 2 l, 



(6.50) 



where Y : — >• M_|_ is the gamma function, and it should not to be confused with the times r n m . The 
first inequality and equality above follow because A mj „ is a sum of u){m) independent mean-T^ -1 expo- 
nential random variables. For the second inequality in (|6.50p . I have used that ^2m=i X^w^o 1 ^m,n < t 
and ui m < sup < r < t l-Krl' 3 - Finally, by Jensen's inequality 



i"4E 



SUp £r 

0<r<t 



sup 

L 0<r<t 



4 /3 1 

< 22t"4E 



r- 



where the first inequality is Jensen's and the second is Doob's. The equality is by Part (1) of Prop. [5] 



(ii). The difference is bounded by 

N-t-l Lm-l 



sup 

0<s<l 



i"4 Y, S ( K rJ E A m)ri (|#r m ,J-|# T J) + 2T*( sup |/Y S |) E A - 



m=l 



n=0 



0<s<l 



N t -1 



3.51) 



m=l 



The second term decays to zero by the argument in (i). By partial summation, the left sum in (|6,5ip 
is equal to 



N 3t -1 



L m — 1 



t-4 Y s(K Tm ) E (l^rUH^r m _J)(r m ,z m -r m , n )« 

m=l n=l 
N s t-1 im-l 

rl E E (I^O-l^r^DCr^-r^^H E F ™> 



N st -1 



3.52) 



m=l 



where -ff m is defined as 



n=l 



m=l 



F m := S(K Tm ) Y (l-ffr m , n | - |#r m ,„_il) (r m + Ar m - r mi „)x(r m + Ar m > r m , n ). 

71=1 

For technical reasons involving conditioning, it will be convenient to replace Kr m n by Kr m „ as in the 

first approximation of (|6.52p . The sum t~i J2ni=i Yju=q 1 &m,n\\KT m>n \ - |-Kr m ,„|| is easy to bound 
using \K r \ £ r and techniques used in (i). The second approximation in (|6.52p replaces Y m ^ m by 
r m + Ar m in the expression, and the resulting error decays by the argument in (|6,48p again. 

Since the probability that sup 0<r < t \K r \ > e -1 ^ is small for 1 S> e > 0, I have the equality 

t~4 Y^f^ 1 H m = t~i X^ff 1 HmX{\Kr m | < e^ 1 ^) with probability close to one. Introducing cut- 
off's will be useful to avoid problems with higher moments of K r . So far I have made only minor 
adjustments to the expression. The strategy to show that sup 0<s<1 X^f!^ 1 ^ m converges in 
probability to zero will be to prove: 
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ii'). riE E^iMH m \T-]\ X (\K Tm \<e-Hr 



0, 



ii")- t~*E ^Ll 1 (E[fl? m | jr.-] -E[H m \ T-Y) X {\K Tm \ < e'H 



0. 



By an analogous Lindberg condition as for , the upper summands above can be taken to be either 
— 1 or Nf depending on convenience. Since the sum of {H m — E[i7 m | -T 7 ,--] )x(l-^r m | < e -1 ^) for 

m G [l,N s t] is a J-j*^ -martingale, the convergence (ii") along with Doob's maximal inequality shows 
that t~4 sup <s<i | X^m^r 1 H m — E[-ffm I -T 7 ,.-] I converges to zero. 

(ii'). A single term |E[H m | ^ T -] | is bounded by: 



E 



oo 

)(|^r m ,J-|^r mi „_ 1 |)E[Ar TO + T m - 



n=l 



T- 



oo 

^-^J |E[^ X (Ar m + r m > T m , n ) {\K Tm J - l^r^J) 



n=l 



oo 

+ ^ 1 |E[^x(Ar m + r m > r m , n )(|^r ro ,J - l^-J) (\Kr^J - \K Tm \) 



71=1 



n OO 

+ C(-) C |^ m rE[^ x (Ax m + r m >r min )||Kr^ L |-|^r m , n _ 1 



n=l 



J 7 - 



where I have applied the triangle inequality with 

E[Ar m + r m -r min | Jr m J =^ 1 |^r m | + [v~ l \K T 



-v~ X \K Tm \) 
+ (E[Ar m + r m -r mi „| Jr ra , n ] - f _1 |^r ra ,„|). 



(6.53) 



(6.54) 



For the third term in (|6.53j) . I have applied an analog of Part (1) of Prop. 16.31 to get a C > for a 
given C > such that 



E[Ar m + r m - Y m , n \ Tv m<n \ - ^ifryj < C|i^J c < C(-)^ Tm |\ (6.55) 

where the second inequality uses that |-Kr m ,J < ll-^rvJ f° r ^m,n < AT m + r m . I have said "analog" 
above because the situation is not identical to Prop. 16.31 since I am conditioning with respect to 
J~r m n ra ther than conditioning that the sign of the momentum does not change on the Poisson time 
following r m n . The choice made earlier in the proof to deal with Tt mn rather than Jr m ,„ is useful 
now, since the conditional density for Kr mn and 8(ifr mn ) given Jr m „ will be bounded. It follows 
by the discussion preceding the statement of Prop. [6731 that I can apply the same proof as in Part (1) 
of Prop. 16.31 to this case. 
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The expectation in the second line of (16.53P can be rewritten as 



E 



J- m,n— 1 

n=l 

oo 

- cE [ ^2 x ( Arm + Tm > Tm ,n-i ) \Kr m ,„ 
n=i 



(\ K r m ,„ \ ~ |-^r m ,„_ 1 |)x(Ar m + r m > T m>n J 

-l l 

a; 



> n - 1 



<2c\K Tm f-^m[Y: 

n=l 

< 2c\K T J^ 1+ <-^E[M ATm+Tm -M ATm + 1 1 f. 



V oj(m) 



uj[m 



<-\K K 

a 



(6.56) 



where replacing x(Ar m + r m > r m , n ) by x(Ar m + 

Tm ^ rm,n-i) on the first line makes the expression 
larger. The first inequality is for some c > by Part (2) of Lem. 16.71 For the second inequality, I 
have used that \K r \ > h\K Tfn \ for r G [r m , T m + Ar m ] , and I rewrote the argument of \ in terms of 
the corresponding counts for the Poisson times. The third equality is from the expectation formula 
Y^n=\ > n ] — 1E[X] for a positive random variable X. The fourth inequality in (I6.56P follows from 
the relations 

E[M A r m+ r m -M Tm \T Tm ]=nK[Ar m \T Tm ] < -^J -1 , 



a 



where the inequality is for large enough times by Part (1) of Lem. 16.71 
For the expression on the third line (|6.53j) . I split the factor \K-p r 
terms \Kr . I — \K T I and \Kr I — \Kr , I. For \Kr 

I L m,n—l \ I 'in I I A m , n I I L m,n — l \ I A m , ti 



\K T I into a sum of the 



K Tm | , I introduce a nested conditional 



expectation with respect to Fr m ,n-i as follows: 

oo 



i PE \K T „ 



l-^r m , n _il I Ft 



n=l 



oo 



n=l 



< -\K l^+g- 1 — 
~ u Tml w(m) 



E[(Ar„ 



w(m) 



71- 2 E 



>n-l ) \K T 



m,n— 1 



l-l^rj 



sup 

0<r<Ar m 



I if- 



T m +r| 



T- 



x 



\K l c 

K V T m | ) 



The first inequality is by Part (2) of Lem. 16.71 and | if r m n _ x | > 1 1 if r m I for r G [r m , r m + Ar m ] . The 
first expectation on the third line is ~ -^|ifr m | by Part (1) of Prop. l6T3l The second expectation on the 
third line is also 0(\K Tm \), which can be shown by approximating E T « \K r \, writing £r m , n — £r mn ^ 1 
in terms of the Doob-Meyer components M r , A r , and applying the relation a = ^p(M, M) r + 2£ r -^A r 



in the standard inequalities as before. The term corresponding to |ifr„ 
analysis yielding a bound 0(|if Tn | 1+< *). 



|ifr m , n _il has a simpler 



The fourth line of (|6.53p is 0(|if Tr J^ 2). By the above analysis, the term |E[i? m | -F.-]] is bounded 
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by a constant multiple of l-K^-ml 1 (since £ < \). Thus, (ii') is bounded by a constant multiple of 



Nt-l 



He[£ \K Tm \^ x {\K Tm \<e-H^)_ 



m=l 



<e-<ri^E[Y^\K Tm \ 

m=l 
Nt-1 

m=l 
Nt-1 



m=l 



(6.57) 



where the inequalities follow by the standard method using Part (1) of Prop. [Ol and £m=l A-r m < t. 
Hence, the expression in (ii') goes to zero for large t. 

(ii"). A similar argument as in (ii') yields that 

.Nt-1 



r^[^V[H m \F T -]\{\K Tn \>e- x £\ 



m=l 



is bounded by a constant multiple of e 2 ^t 1+< >. Thus, this term also vanishes as t — > oo. The 
expression (ii") requires an analysis of ~E[H^ m | F T -] • I begin with the inequality 



\E\H? \ < 2 

L r m I r m J I — 



E 



l<n'<ra<L m 



+ 



E 



[ E ^ 



Kn<L„ 



where G n ' in has the form 

G n> := x (Ar m + r m > r m , n ) (|ifr ra>n , | - l^r^.j) (l^rvj - I^.J) 



x E 



(6.58) 



The conditional expectation in (I6.58P can be written as 



E 



(A"7~jti ~i~ Tin Pm,n) {^m,n ^m,n') (^^m ~i~ 1~ m ^m,n) 

={2u- 2 K 2 Tm + ^(rVn - r m , n ,) |i^ Tm |) 

+ (2z,- 2 (^ m n - i^J + u-\T m)n - r m , ft /) (|Kr roi „| - |*T T 



where the error term is bounded by 



.59) 



l^|<c(^) 1+c |K rm |^ + c(^|r min -r min ,|K m |C, 

for C, £ > defined as before. This is analogous to (I6.54D and (16.55H in (i'), and follows by Part (1) 
of Lem. 16.531 and because |-Kr mn | < |l-^r m | for T m ^ n < Ar m + r m . I will not go through the details 
for the bounds of the three terms on the right side of (|6,58p . since they admit the same procedure 
to break them down as applied before (except, for instance, the inequality X]?S=i^[^ > n \ — E[X] 
is replaced by Y,n=i nF i X > n ] < E[X 2 }). The end result is that E[H* m \? T -\ is smaller than a 
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constant multiple of |X Tm | 3 , where the leading term comes from the diagonal sum n = n' . By the 
same argument in (|6.57p . the sum of terms |if Tm | 3 is bounded by 



Nt-1 



[ £ i*v„ 



\(\K Tm \<e-H 



m=l 



which goes to zero for large t. 



(iii). I will first show that there is a vanishing error in replacing the terms E[Ar m | ^F T -] in the 
expression by u \K Tm \. To bound the difference, I apply Part (1) of Prop. [631 to get the first and 
second inequalities below for C > depending on my choice of < ( < %: 



N 3t -1 



E 



sup 

0<s<l 



r% E \ K r m \\^ m \T T -\-v- x \K T „ 



m=l 

< CE 
2Cv 



N t -1 



^ E \ K rn 



|i+C 



m=l 



< 2CVE 



N t -i 



^ E \KrJ C E[Ar m \T Ti 



m=l 



< 



sup \K r \< 

I *-0<r<t 



2Ci/ 




< E 


sup -RT r | 2 


U 


l-0<r<t -1 



* I 2 



4Ct/(E[g 2 ] + 
i 



0. 



The third inequality follows by removing the nested conditional expectations and Y^rn=i Ar m < t. 
The fourth inequality is Jensen's, and the fifth employs \K r \ < £ r and Doob's maximal inequality to 
the positive submartingale £ r . The equality is Part (1) of Prop. 15.21 

I am left to bound the sum Y2 m =i K Tm \K Tm \. Let G C M+ be the set of all times r m such that 
T m S Wj + i) and T m is duii even-nuinbered flip time following (and this includes T m — ). Denote 
the number of sign-flips in the interval (sj,Wj + i) by n^. The sum of terms K Tm \K Tm \ can be written 
as 



£ E 

j— 1 T m e[?j, wj+i) 

r m +i<st 



K Tm \K Tm \ 



m=l 

r m eG 



E ( K Tm\ K r m \ + K Tm+ATm \K Tm+Arm \\ - X(nj even)K Tm+ Ar m \K Tm+A r, 



r s t-i 



(6.60) 



j'=i, 

fj — T TO 



where I have neglected a single extra boundary term ET Tm | i^ Tm | that may occur in the last incomplete 
excursion. 

I can immediately treat the boundary sum on the right side of (|6.60p . Since r m + Ar m occurs 

3 

during an incursion, I have that \K Tm+ A Tm \ < is, and thus the first inequality below holds: 

Tt-l 



t"4 



ME 



I ^ 



r m +Ar„ 



< t-5E[T t ] < 2ah 



1 3 
8 . 



The second inequality is from (|6.2j) in the proof of Lem. IB~T1 and relies on the submartingale upcrossing 
inequality. 
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For the first sum on the right side of (I6,60p . it is convenient to write K Tm \K Tm \ + K Tm+1 \K Tm+l | as 
a sum of the two terms 

(S(K Tm ) + S(K Tm+Ar J)\K Tm+ATm \ 2 and - S(K T J (\K Tm+Arm \ 2 - \K Tm \ 2 ) . (6.61) 
I can bound the sum of terms on the left of ()6.61|) through the inequality 

, N st -1 



sup t-i\ {S(K Tm ) + S(K Tm+ATm ))\K Tm+ATm \ 



0<s<l 



m=l 



Ni-l 



0<r<t 



<2t"*( sup \K r \ 2 ) ^ x {S(K Tm ) = S(K Tm+ATm )). 



(6.62) 



m=l 



For any 5 > 0, the factor t 1 <5 sup 0<r<< \K r \ 2 will go to zero for large t. By definition, S(K Tm ) = 
S(K Tm+Arm ) occurs when \K Tm+Al - m \ jumps out of the interval [||if Tm |i |l^r m |] • By inserting nested 
conditional expectations, I have the equality below 



Nt-l 



Nt-l 



H+ 5 E[ £ x{S(K Tm ) = S(K Tm+AT J)] =t-i +5 E[ P[\K Tm+Arn \ i [\\K Tm \, \\K r J[ 



m=l 



- g t l-S L * 



(6.63) 



where the inequality follows since |-Kr m | > for all times r m and by Part (4) of Prop. 16.31 The 
expectation E[Nt] has a rough bound given by routine arguments: 



Nt-l 



E[N t -l] < r«E[ \^ \K Tr 

m=l 
N t -1 

m=l 



N t -1 



< 2;'/ _ ?e[ ]T E[At,„ ! J 7 . 

m=l 

5 

< «5. 



The first inequality is from \K Tm \ > t&. For the second and third inequalities, I have applied Part (1) 
of Prop. 16.31 and Yjnv=i ^ T m < t, respectively. Hence, the right side of (|6.63|) will go to zero for any 



choice of <5 < |. 

The sum of terms 



N st -1 

Y, -S(K Tm )(\K Tm+AT J 2 -\K Trr 



m=l 

T m €G 

on the right side of (|6.6ip is easier, since I can approximate \K r \ by the submartingale S r = M r + A r 
as before. It is convenient to write \K Tm+Arm \ 2 — \K Tm \ 2 as a sum of (\K Tm+ATm \ — |-^r m |) an d 
2(|-^T m +AT m | — l-^T m |)|-f^T m | and to treat the two terms separately. 

□ 

Proof of Leva. HOI I will show Part (1). Part (2) follows by similar analysis relying on Part (5) of 
Prop. 16.31 Part (3) is a consequence of the convergence in law of t~2M st over s E [0, 1] as t — > oo by 
Thm.EU the Lindberg condition for t~^M st in Lem. E3J and [HI Cor.VI.6.7]. 
The quadratic variation of mW has the following form: 

N st 



[m% =Hj] Kl (Ar m - E[Ar m | F T -]f . 



m=l 
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For technical convenience, I will work with the expression on the right with the upper summand N s t 
replaced by N s < — 1. The difference is negligible by Lem. 13.41 I will show the following convergences 
in probability: 

v2 



W- su Po< s <i 
(ii). sup 0<s<1 



H £^rM <^r m - E [Ar m | )< 



~ X \K I 3 At 



0, 



.-u 



vt 



0. 



(i). For small e, P[sup 0<r<t |t 2 -KV| > e x ] is small. Hence, for K r = K r x{\K r \ < e 1 ), 



N st -1 



E <( A ^- E [ A ^KJ)' 



(6.64) 



m=l 



is typically equal for all s £ [0, 1] to the same expression with K r replaced by K r . The advantage of 
working with K r is that it has arbitrarily many moments. The following is a martingale with respect 
to the filtration J-^p: 



N, 



m=l 



:= H E ^ m (Ar m - E [Ar TO | £ ]) 2 - E (Ar m - E[Ar m | j; ]) 



J" 



(6.65) 



Using Part (1) of Prop. 16.31 for the first two inequalities below, the second moment of the above 
martingale is bounded through the inequalities 



N st -1 



E 



t- 5 E 

32 



771=1 



Ar m - E[Ar m | .F r -]) 2 - E[(Ar m - E[Ar m | F T -]f \ 



N st -1 



< — E 

77^ 



e * 



64 
< -?E 

64 7 i 
/y 3 



m=l 
N st -1 



^ 5 E I^J 7 E[Ar m |^ 



m=l 



I have multiplied the bounds in the first two inequalities by 2 to cover the error terms for the approxi- 
mations of the moments E [(AT m ) m | f r -\ from Prop. 6.3. The last inequality follows by \K Tm I < €~H 1 

and removing the nested conditional expectations and Y^m=i Ar m < t. By Doob's maximal inequal- 
ity, E[sup < s <i |wi^' e | 2 ] tends to zero for any fixed e. By similar applications of Part (1) of Prop 



as above, the difference 
N.t-1 



E 



sup 

0<s<l 



-l 



K Tm \E[Ar m \T- 



t" 2 E Kl{M[(Ar m -E[Ar m \T T -]y 

m=l 

tends to zero. The difference between 

v-h-* i^i 3E [ Ar - 1 Ft-] and v ~ Xt ~^ E i^i 3At » 



0(t ? "5) 



m=l 



m=l 



is a J-s^ -martingale and tends to zero by a similar (but simpler) argument as for (|6.65p . Finally, the 
process v 

Em=i 1 l^r m | 3 Ar m is equal to the same expression with \K Tm \ replaced by \K Tm \ with 
probability arbitrarily close to one as e — > 0. 
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(ii). First, I show the sum over |iC Tm | 3 Ar m can be replaced by a sum over J n 
approximation reduces to some martingale analysis using that \K r \ ~ £ r : 



Tin 
T-m 



dr\K r 



This 



N st -1 



m=l 



Tm ~T ^A Tm 



N si -1 



J Ar n 



m=l 



Nai-1 / . Tm +Ar m 



< sup \K r \ 2 (t~i V / 



dr I if r 



m=l Tm 



(6.66) 



For any 5 > 0, the random variables i -1-5 sup < r < i \K r \ 2 tends to zero, since t~z\K st \ converges to 
the absolute value of a Brownian motion in the uniform metric. The quantity \K r \ in the sum on the 
right side of (16. 66ft can be replaced by £ r with an error 0{\K r \~ l ). Moreover, I have the inequality 



Nft-1 f T m +Ar m 



sup 

L 0<s<l 



E 



m=l Tm 



dr\£ r - £ Tn 



Nt-l „ T » 
m=l 'm 



rfr M r - M T 



J~ T i 



/~r m +A-r m 



[E 

771=1 



(6.67) 



where t^,t^' are the Poisson times following r n and r m + Ar m , respectively. Both t^,t" are hitting 
times with respect to the nitrations J>, J>. For the inequality above, I have used the triangle inequality 
with £ r = M r + A r and introduced nested conditional expectations for the martingale term. Changing 
the lower bounds of the integration from r m to T' m for the martingale term doesn't change the value, 
since M r = M Tm for r £ [r m , r^J. 

For the drift term, I can use that a = -^(M, M) r + 2£ s -^A r to get the bound A r -A Tm < dv^. 

3 

However, since £ v > t~s over the excursion periods, the first inequality below holds: 



^t- 1 /■r m +Ar m 



E 



[E 



777= 1 JTm 



Nt-l 



dr\A r - A T 



< arm 



[ E( A ' 



771=1 



2 5 




< -t?E 


sup |if r | 


V 


-0<r<i - 1 



0(ti 



The second inequality uses that Ylnv=i Ar m < i and Part (1) of Prop. [631 through the standard tricks: 



N t -l 



Nt-l 



E 



771=1 



771=1 

, Nt-l 

-E[ £ Ar m i^ rm 

777=1 



Nt-l 



771=1 



[ J2 (Ar m ) 2 ] = E [ £ E[(Ax m ) 2 | < ±l[ £ K*J < i E [ £ ^r m E[Ar ra | ^ 



Nt-l 



777=1 



8t 
< — E 



sup \K r 



v l 0<r<t 



where the first and second inequalities hold for large enough t. Thus, the second line of (|6.67|) vanishes 
as t — > oo. The martingale term is handled by similar arguments. 

By the above, I can work with the sum of the integrals J Tm+ Tm dr\K r \ 3 . By the definition of the 

m 3 

low energy incursions, the momentum has the upper bound \K r \ <2ts and thus 



sup 

0<s<l 



N st-1 rTm +Ar m 



E 



dr\K r 



i :> 



-1 J Tm, 



dr\K r 



< - f 
v 



The final remainder sup < s <it 2 dr\K r \ 3 for the difference between t 2 Jg TNst (irl^l 3 and the 
final expression t~2 dr\K r \ 3 is smaller than (sup Tm<t Ar m ) sup 0<r<t |iT r | 3 , which vanishes by the 
same argument as for the proof of the Lindberg condition for mi in Lem. 13.41 

□ 
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A Existence and uniqueness of the quantum dynamical semigroup 

Lindblad equations with unbounded generators can pose technical difficulties regarding the existence, 
construction, and uniqueness of their corresponding quantum dynamical semigroups. This situation 
is analogous to that for Kolmogorov equations in classical Markovian dynamics. A quantum problem, 
however, would be considered "solved" if it was somehow reduced to a classical problem, such as in 
the case of models that are fully translation invariant [27j . 

As before let D(if ) C L 2 (M) denote the domain of the Hamiltonian H. The mathematical definition 
of the quantum dynamics will require the form generator C : D(£f) x i3i(L 2 (lR)) x D(iT) — > C given 

by 

£(Vi; P ; ih) := (QfT - 2- 1 **(/))^ 1 | / 9^} + (ih\p(jH - 2- 1 **(/))^ 2 } + (fa; *(/#2>. 

The form generator is designed to draw the operation of the unbounded terms away from p and on 
to the vectors ipi,ip2, which have a restricted domain. In my case, the map $ has finite operator 
norm, since by the complete positivity of \&, its operator norm is equal to the norm of € 
£>(L 2 (IR)). Conveniently, = 1ZI has norm 1Z. If \I/ were an unbounded map, I would have 

to look for a convenient Kraus decomposition ^(p) = ^ ■ A*-pAj, and the last term above would be 

A semigroup of maps $t on /3i(L 2 (M)) is said to be conservative if Tr[$t(p)] = Tr[p] for all 
t E M+ and positive p G £>i(L 2 (M)). Alternatively, this can be expressed by the adjoint semigroup 
as t (J) = /. The representation of the semigroup $A,t in (2) of Lemma lA.ll is a pseudo-Poisson 
property in analogy with classical Markovian semigroups |19t X.l]. The map $\ t is an expectation 
over the Hamiltonian flow interrupted at discrete random times by the operation of a transition map 
T := 7£ _1 \£', where the times occur according to a Poisson clock with rate 1Z. 

Lemma A.l. There is a unique strongly continuous, conservative semigroup of maps $>\ >t on £>i(L 2 (IR)) 
satisfying 

(</>i|$A,tO)V>i> = (HP^i)+ [ drC^QxApy^z) ( A -!) 



o 



for all ipi,ip2 G D(i/) and p G Bi{l/ 

1. The semigroup can be written in the form &\ : t(p) = E[f7\,t(£)P^A where f7\,t(£) = 

e a g^nA . . . e a e 1 " 1 , the expectation is with respect to a Levy process with rate 
density j(v), and £ = (ti,ui; . . . ;t n ,v n ) is the realization of the process over the time interval 
[0,t]. 

,2. Alternatively, $>\.t(p) = E ^ (p)] , where the expectation is with respect to a Poisson clock with 

rate 1Z, the Poisson times over the interval [0,t] are £ = (ti, . . . ,t n ), and : £>i(L 2 (IR)) is 
defined by 



:= ft~ n e ~ • • e — H ^(e~— H pe— H )e^^ H • • • )e~^~ H 
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Proof. By [17] ? solving the integral equation (|A.1|) is equivalent to solving 



In most unbounded cases, there is a technical issue in checking that A x : -j^ff - 2 _1 **(I) defines 
an m-accretive operator over an appropriate domain, and therefore generates a contractive semigroup 
on L 2 (R). However, since = 1ZI, this is trivial given that H is self-adjoint, and I can factor 

I a i 7Z. it tt it jj 

e tA\ = g-'ig j Since * and are bounded, a solution to the above integral equation can be 
constructed by the Dyson series 

oo . 

^x,t = e- tn Y j n n (A.2) 

n=0 J 0<h... t n <t ' 

where <£>j ^ and £ are defined as in statement (2) of the lemma. The semigroup <J?A,t is strongly 

— — H 

continuous, since \I> is bounded and the group e >~ is strongly continuous. To check conservativity, 
it can be computed that <3? Ai (I) = I using = 1ZI and that e~^ H is unitary. 

The summation on the right side of (|A.2p is equal to E (p)] , where the expectation is with 
respect to the a Poisson process with rate 7Z. The other stochastic representation for <&\ t is obtained 
by expanding ^ in its integral form. The resulting integrals can be commuted, because the integrands 
are completely positive maps. 

□ 



B The fiber decomposition for one-dimensional periodic Schrodinger 
equations 

The reader is directed to [38] for a more detailed discussion of the structure of periodic Schrodinger 
equations. A Schrodinger Hamiltonian P 2 + V(X) with a period-27r potential satisfies 

e i27rP {P 2 + V{X))e~ i2nP = P 2 + V(X + 2vr) = P 2 + V(X). 

Commuting with e l2nP implies that P 2 + V(X) must have invariant spaces corresponding to the 
spectral values for e l2nP (i.e. the unit circle in C). The Hilbert space H = L 2 (R) admits a fiber 
decomposition 

% = r i fyUfr -H^ = L 2 {[-tt^)), (B.l) 

with fiber-maps sending / G L 2 (R) to [f]^ G L 2 ([— 7r,7r)) and formally defined according to the partial 
Fourier transform 

[ftW^^e-^/t. + a™), s€[-7r,7r). 
( 27r ) 2 nez 

From the Fourier transform formula above, the eigen space property for e l2nP is clear: [e l2nP f\Ax) = 
e l2v ^[f]^. The parameter <p is called the quasimomentum or the crystal momentum, and its domain 
[—5, 5) is the first Brillouin zone. The operation of the Hamiltonian H on the </>-fiber is [Hf]<f, = H^[f\^ 
for 

H * = + 
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d 2 



where V(x) is interpreted as a multiplication operator on L 2 ([— it, 7r)) , and ( da 
with boundary conditions 



is the Laplacian 



e i2 ^ 5 (7r) 



and 



d<7 ^ 



-7T 



e i2 ^^(vr). 
ax 



For each G [- 
non-degenerate for 
frand index. When 



2' 2 



i), the self-adjoint operator has compact resolvent. The eigenvalues are 



k G R — iZ through the relations 



/ -|,0 and are labeled progressively as E n & by a parameter n G N called the 
^ -±,0, the pair (n, 0) is related in the extended-zone scheme to a parameter 



>mod 1, 



and 



n 



2\k-<j>\ S{k) = S{<f>), 
2|fc-0|-l S(k) = -S {</>), 



(B.2) 



where S : K — )■ {±1} is the sign function. The assignment of for 6 { — ^,0} is a matter 

of convention that I am not concerned with, since the set {— ^,0} has measure zero in the direct 
integral (jB.ip . The dispersion relation E : M — > M + is defined as E{k) = E n ^ for (j) ^ — 1,0 and 
related to (n, </>) as above, and I hold the convention that E{k) is symmetric and left-continuous for 
k > 0. 

For / -i,0, let ^ 

nj be a normalized eigenvector for with eigenvalue E n ^. I can pick 
the eigenvectors ip n ,<j> to vary continuously (and, in fact, smoothly [381 Thm.XIII.90]) as elements 
in L 2 ([— 7r,7r)) for <fi G (— 7r,0) and G (0,vr). Given f £ W, I can assign an extend-zone scheme 
representation / G L 2 (R) through 



/(*0 = W>n,*|[/]*>, 



fc G 



where n,^ are determined by k as above. Again, the assignment of f(k) for k G ^Z is arbitrary. In 
analogy with the position and momentum operators, I can define a self-adjoint operator P Q that acts on 
the domain {/ G Ti \ J^dk\f(k)\ 2 < oo} as multiplication in the extended-zone scheme representation 

(P~f)(k) = kf(k). 

By standard operator calculus, I can define functions of P Q , and the Hamiltonian is given by H = 
E(P Q ). For k G R — ^Z, I define a tempered distribution \k) Q such that for an element / G S(R) in 
Schwartz space, 



(f\k) Q ■= £ ^ [ dx f( X + 2 ™)^nA 



— 7T,7r) 



dx [f]<j,{ x )^n,4>{ X ) '■= /( fc )> 



where (n,(f)) is determined by k. In the usual senses, I have the formal relations Q{k'\k) Q = 5(k' — k) 
and P Q \k) Q = k\k) Q . 



B.l Eigenket conventions for the Dirac comb 

There remains much choice in the phase for the eigenvectors Vn,</> G L 2 ([— -zr, 7r)) determining the 
extend-zone scheme of the last section. For the Dirac comb, I will fix the definition with a specific 
closed expression for the eigenvectors. With the correspondence of (n, <fi) G Nx (— \, ^) with k G M— ^Z 
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in (|B.2h . I will replace the subscript by k: ipj. '■= 4> n ,<t>- The eigenfunctions tpk £ L 2 {\— tt, tt)) for 
k G R — ^Z are given by 



' eU^W-V-l r -ixQ(k) | c i27rfqffc)-fc) r kcq(fc) 



^fc(s) = iV fc ~ 



i27r(q(fc)-fc)_i p -ig;q(fc) + e ixq(fc) 



l_ e -i2ir(fc + q(fc)) 



-7T < X < 0, 



< X < 7T, 



where iVjfc > is a normalization, and q : R — > R is defined as in (|2.4p . For |fc| S> 1, I approximately 
have that q(fe) ~ A;, and the Bloch function ^ is approximately the plane wave (2tt)~2 e lxk , except 
when k + q{k) ~ 2fc is near an integer. 



B.2 The diagonal in the extended-zone scheme representation 

A density matrix p G £>i(%) determines a probability density [p] D G L X (R) corresponding to the 
distribution in the extended-zone scheme variable. Intuitively, this is given by the diagonal of the 
integral kernel p(k±,k2) := q(^i|/o| ^2)^, although kernels are only defined a.e. R x R for Hilbert- 
Schmidt operators, so this does not offer a rigorous definition without some additional condition on 
the kernel such as continuity. For a rigorous definition, notice that there is a unique probability 
measure p p on R such that for all g € L°°(R) 

Tr[g(P Q )p] = [ dp p (k)g(k). 

This follows by the Riesz representation theorem, since the left side is positive for g > 0, bounded 
in absolute value by ||<?||oo) and equal to one for g = 1r. By the continuity of the spectrum of P Q , 
the measure must be continuous with respect to Lebesgue measure, and I denote the Radon-Nikodym 
derivative of p p by [p] D G L 1 (R). 



B.3 Invariant fibers for the Lindblad dynamical semigroup 

The rate of kicks from the gas is invariant of the spatial location of the particle. The total dynamics 
is thus invariant under spatial shifts by 2ir. For the dynamical maps $>\j : /3i(L 2 (R)), this feature is 
expressed as the covariance 

^, t (^ p pe-^ p ) = e^ p ^ x ,(p)e~^ P . (B.3) 

Not surprisingly, this implies that the Banach space i3i(L 2 (R)) decomposes into invariant fibers in- 
dexed by the Brillouin zone. This can be understood formally by the statement that the kernel values 
Q (k\p\k + n + 4>)q for k G R, n G Z do not interact dynamically for different 4> G [— ^, This holds 
also with the kets \k) Q replaced by the standard momentum kets \k). Mathematically, it is easiest 
to discuss the invariant spaces for the adjoint semigroup §* xt . Let At C B(L 2 {R)) be the algebra of 
all bounded operators commuting with e l2nP . By the analogous covariance property ()B.3|) for &X t , 

it follows that ^(Aj) C Aj. More generally, the Banach spaces e^ x Aj for <t> G [-|, |) will also be 
invariant by the Weyl commutation relation 

g i27rP e icj,X e -i27rP = ^vr^X _ 



C Dispersion without the Dirac comb 

The elementary lemma below implies that the spatial dispersion for the particle scales as 1 2 for times 
t > 1 when the Dirac comb is not present. This scaling agrees with the classical case. 
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Lemma C.l. Let p\t G £>i(L 2 (R)) satisfy the Lindblad equation i2.1\) with a = 0. Also let Q\j{x) :- 



3 3 

(x\p\ t\x) be the position distribution for the particle. As t — > oo, the renormalized density t^Q\t(t^x) 
converges in law to a Gaussian with variance j^j. 

3 3 

Proof. I will show pointwise convergence as t — > oo for the characteristic functions (pxj of t^Q\ t t{t^x). 
The function (px,t can be written in terms of a trace formula involving p\ t by the following: 

VX ,M := / dxtlQ x , t (tlx)e ivx = [ dx Qx,t(x)^ vx = Tr [^ hx p X ,t] , 

JR JR 

where X is the position operator. The right side is a special case of the quantum characteristic function 
for the matrix p\ t . The quantum characteristic function has a closed factored form given by 



Tr[e ivX+iqP p x ,t) = Tr[$* Xt (e ivX+iqP )p) = e ^ dr (^ q+ ^ {t ~ r ^~ m ) Tr[e ivX+iqP 



P 



where v,q G R and 4>{q) := J R dvj(v)e lvq . The second equality above can be seen through the 

expressions for Q\ t in Lem. |A.1| Weyl's intertwining relations, and the relation e~ p f(X)e~~ p = 
f(X + %rP) for bounded functions / : R — > C. With the above formula, 



e 3A2 , 



since ^'(0) = and <p"(0) = a. Hence, t-2Qx,t{t 2 x) converges in distribution to a Gaussian with 

4(T 



variance 



□ 
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